Journal 


of 
ROYAL AERONAUTICAL SOCIETY 


MARCH _ 1953 


CONTENTS 
SECRETARY’S NEWS LETTER 
NOTICES 


THE SCOPE AND LIMITATIONS OF THE 
PHOTOELASTIC METHOD OF STRESS ANALYSIS H. T. JESSOP 


THE RESISTANCE TO AIR FLOW OF 
WIRE GAUZES W. J. D. ANNAND 


EXPERIMENTS WITH FLAT PLANES (1866-1890) J. L. PRITCHARD 


STRESS ANALYSIS OF CIRCULAR FRAMES IN A 
NON-TAPERING FUSELAGE K. J. DALLISON 


TECHNICAL NOTES 
C. H. McKOEN: J. LOCKWOOD TAYLOR: C. F. TOMS 


THE LIBRARY : REVIEWS REPORTS 
TEN SHILLINGS MONTHLY 


LONDON 
4 HAMILTON PLACE 


wy 
| (= mar 251953 
THE 
| 
EES 
BS 
= 
\ 

WI 


Control Panel 


for Gas turbine 
Starters 


This panel automatically controls the 
whole gas turbine starting cycle, once the starter 

button is depressed. Protective devices ensure that a 
second start cannot be made untii sufficient time has 
elapsed to allow for the possibility of delayed cartridge 
ignition. The panel also protects the turbine starter 
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faster than the safe starting speed. 


For future developments in Turbine Starting — look to Ple S S ey 
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THE ROYAL AERONAUTICAL SOCIETY 


Scholarships and Awards 


The awards administered by the Society fall into three main groups :— 


(a) Scholarships and Educationai Grants, (b) Medals awarded for achievement in the 
profession of aeronautics, (c) Medals and prizes awarded for papers read before. 


or published by, the Society. 


In addition to the named awards there are special 


premiums attached to the Wilbur Wright Memorial Lecture and the British Com- 

monwealth and Empire Lecture. Also, there is an annual sum of £250 from which 

premium awards are made to authors of those papers published in the JOURNAL 
which are considered to be of particular merit. 


SCHOLARSHIPS 


The Charter Scholarship—Awarded annually, or at such 
other times as the Council may decide, to assist a student 
wishing to undertake advanced work or study in aero- 
nautics. An extension for a second year is looked upon 
favourably by the Council if it is clear that the student and 
the work will benefit. The scholarships are awarded only 
to those who can satisfy the Council that they have the 
necessary qualifications to benefit and whose proposals for 
work or further studies are acceptable to the Council. 
Normally the work should lead to a higher degree or post- 
graduate diploma. 


The Geoffrey de Havilland Memorial Scholarship— 
Awarded annually, at the discretion of the Council, to assist 
those engaged in aeronautical research, either theoretical 
or practical, to undertake advanced study in the field in 
which they are engaged. The scholarship may be extended 
for a second year at the discretion of the Council. Scholar- 
ships are awarded only to those who can satisfy the Council 


that they have the necessary qualifications to benefit and | 
whose proposals for work or further studies are acceptable | 
Applications should either be graduates | 


to the Council. 
or have reached graduate standard. 


The Busk Studentship in Aeronautics Awarded annuall) 
at the discretion of the Council to enable the holder to 


engage in research or preparation for research in aero- | 


nautics of the kind which Edward Tashmaker Busk made 
his own, namely Stability, Meteorology and Gust Research, 
Aeronautical Instruments, Sighting and Appliances for 
Bomb Throwing, the calculations and 
problems relating to these subjects and generally, an) 
scientific work in connection with any branch of art or 


trade which conduces to progress in, or the use of flying | 


The Studentship is awarded for one year but may be 
extended at the discretion of the Council. 
British Subject less than 25 years of age. 


ACHIEVEMENT IN AERONAUTICS 


Society's Gold Medal—The highest honour which the 
Society can confer for work of an outstanding nature in 
Aeronautics. 

Society's Silver Medal—Awarded for work of an outstand- 
ing nature in Aeronautics. 

Society's Bronze Medal—Awarded for work leading to an 
advance in Aeronautics. 

British Gold Medal—Awarded for outstanding practical 
achievement leading to advancement in Aeronautics. 


AWARDS FOR 


The George Taylor (of Australia) Gold Medal—Awarded 
annually, at the discretion of the Council, for the most 
valuable contribution read before, or received by, the 
Society on Aircraft Design, Manufacture or Operation. 


Simms Gold Medal—Awarded annually, at the discretion 
of the Council, for the most valuable contribution read 
before, or received by, the Society on any subject allied to 
aeronautics, such as, structures, meteorology, metrology. 


Herbert Akroyd Stuart Memorial Prize—Awarded, at the 
discretion of the Council, for the most valuable contribution 
read before, or received by, the Society on Applied Thermo- 
dynamics. 


Edward Busk Memorial Prize—Awarded annually, at the 
discretion of the Council, for the most valuable contribution 
read before, or received by, the Society on Applied Aero- 
dynamics. 


British Silver Medal—Awarded for practical achievement | 


leading to advancement in Aeronautics. 


Wakefield Gold Medal—Awarded annually, at the discre- 
tion of the Council. to the designer or inventor of an) 
apparatus tending towards safety in flying. 


R. P. Alston Memorial Prize —Awarded for practical 
achievement associated with the flight testing of aircraft. 


PAPERS 


Pilcher Memorial Prize—Awarded annually, at the disere- © 


tion of the Council, for the most valuable paper read b) 
a Graduate or Student during the previous year at an) 
meeting of the Society or its Branches. 


Ushorne Memorial Prize 


Awarded annually, at the discre- 


tion of the Council, for the best contribution to the Societys) 
publications written by a Graduate or Student on some} 


subject of a technical nature in connection with Aeronautics 
Orville Wright Prize—This prize is offered annually for the 
best contribution received for publication in the AERO 
NAUTICAL QUARTERLY of the Society on some subject of 4 
technical nature in connection with Aeronautics. 


Branch Prize—The Council offer an annual prize for the 
best paper read before the Branches during the previou' 
lecture session. The Prize is open to any member of the 
Society or of any Branch. 
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AST year two new branches were formed, Merthyr 
Tydfil and Henlow, and already this year | am 
able to announce the formation of the 24th Branch of 
the Society at Singapore. Squadron Leader A. D. 
Miller has acted as the interim secretary, and in this 
office he has been confirmed. Air Marshal R. O. Jones 
has accepted the invitation to give the inaugural address. 
which is to be delivered early in March. They are to be 
congratulated on the speed at which they have begun to 
function. I know that every Member of the Society will 
wish them well. 


In this JOURNAL there is printed a letter which I 
received from the Secretary of the New Zealand Branch. 
It was so interesting to me that I have had it printed for 
your interest. It seems to me significant that Mr. 
Coleridge has mentioned that the Society helps to find 
and cement friendships as well as fulfilling its primary 
functions of disseminating technical knowledge. It is a 
music hall joke that where two old boys of the same 
school foregather inevitably there is the foundation of 
an Old Boys’ Association. It may be in the future 


_ when you have aeronautical engineers meeting in like 


circumstances it will result in the formation of Branches 
and Divisions of the Royal Aeronautical Society. 
Shakespeare sums up—‘a consummation devoutly to 
be wished.” 


Council have agreed that exemptions to the 
Associate Fellowship Examination be given for 
academic qualifications obtained at the College of 
Aeronautics, Cranfield, The Technological Institute. 


_ Loughborough College, and the Royal Naval Engineer- 


ing College, Manadon. These exemptions are for those 


obtaining the qualifications after Ist January 1953. 
_ Certain conditions have been laid down by the Council 


and accepted by the establishments. 


Secretary's News Letter 


March 1953 


The British Commonwealth and Empire Lecture will 
be given on Thursday 11th June by Sir Hubert Walker, 
C.B.E.. Chairman of West African Airways Corpora- 
tion. It will be followed by a reception at 4 Hamilton 
Place. This will be open to all Members of the Society. 
Particulars will be given later. 


The Garden Party will be held this year at Hatfield 
on Sunday 14th June. The theme will be “history.” I 
hope that there will be on show, if not flying, some of 
the aircraft which are important in the long roll of 
aeroplanes. I shall be very pleased to learn of the 
whereabouts of any old, or even not-so-old, aircraft. 
which might otherwise be overlooked. 


The Society will celebrate the SOth anniversary of 
the flight by the Wright Brothers by a Dinner which will 
be held in the Dorchester Hotel on Thursday 17th 
December 1953. I hope that the many who have con- 
tributed and are contributing to the 50 years of aviation 
history will be able to be present. Divisions and 
Brariches might also consider the possibility of holding 
a Jubilee Dinner. perhaps on a Regional basis. 


The Annual Dinner of the Portsmouth Branch was 
held on 27th February. This Branch has been some- 
what handicapped by the departure of the design staff of 
de Havillands (Airspeed) from Portsmouth to Christ- 
church. Despite this their dinner was very successful, 
and both the Chairman, Mr. G. H. M. Twyman, and 
the Secretary. Mr. E. M. Bellamy, are to be congratu- 
lated. There may be some. members of the Society. 
living locally, who are not members of the Branch. 
There are certainly some aeronautical engineers in 
Portsmouth who are neither members of the Society nor 
the Branch, and these ought to give their support. 


Secretary 
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NOTICES 


An annual sum of £250 is available for premium awards, usually 15 guineas 


each, for papers published in the Journal. 


Members and non-members of the 


Society are invited to submit Papers on any aspect of aeronautics 


ANNUAL GENERAL MEETING 7TH May 1953 


Notice is hereby given that the Annual General Meeting 
of the Royal Aeronautical Society, with which is incor- 
porated the Institution of Aeronautical Engineers, will be 
held on Thursday, 7th May 1953, at 5.30 p.m. in the offices 
of the Society. 4 Hamilton Place, London, W.1. 


AGENDA 


1. To read the Notice convening the Meeting. 

To receive and deliberate upon the Report of the 

Council on the state of the Society and the Balance 

Sheets and Income and Expenditure Accounts of The 

Royal Aeronautical Society and Aeronautical Trusts 

Limited for the year ended 31st December 1952. 

3. To receive the names of those elected to Council for 
the years 1953-56. 

4. To announce the names of Fellows elected by the 
Council in accordance with By-Law 4. 

5. To elect the Auditors for the year 1953. 

6. Any other business. 


to 


By Order of the Council. 


A. M. BALLANTYNE, 
Secretary. 


NOTE—In accordance with the By-Laws any member 
whose subscription ‘has not been paid before the first day 
of April is not entitled to vote. 

Light refreshments will be served after the meeting. 


EASTER HOLIDAYS 


The offices and Library of the Society will be closed from 
5 p.m. on Thursday 2nd April, until 9 a.m. on Tuesday 
7th April 1953. 


NEW BRANCH 
A Branch of the Society has been formed in Singapore. 
The address of the Honorary Secretary is : — 
Squadron Leader A. D. MILLER, R.A.F., A.F.R.Ae.S.., 
Honorary Secretary, 
Singapore Branch of the Royal Aeronautical Society. 
c/o Officers’ Mess, Royal Air Force, 
Seletar. Singapore 28. 


THE Forty-First WILBUR WRIGHT MEMORIAL LECTURE 


The Forty-First Wilbur Wright Memorial Lecture will 
be given by Professor N. J. Hoff, F.R.Ae.S., F.I.A.S., 
Professor and Head of the Department of Aeronautical 
Engineering and Applied Mechanics, of the Polytechnic 
Institute of Brooklyn, New York, on “ Structures.” The 
Lecture will be given on 14th September 1953, on the eve 
of the Anglo-American Aeronautical Conference (15th- 
17th September 1953). 


1953 GARDEN PARTY 


This year the Society’s Aeronautical Garden Party will 
be held at Hatfield Aerodrome (by kind permission of the 
Directors of de Havilland Aircraft Ltd.) on Sunday 14th 
June 1953 in the afternoon. 


Further particulars will be published in the next Journal. 


NOMINATION OF CANDIDATES FOR COUNCIL 
The following is an extract from the By-Laws:— 


“The Twenty-one ordinary members (of the Council) 
shall be nominated and elected from among the members 
of the Society. At the date of their election at least ten 
shall be Fellows, and one at least shall be in each of 
the following classes: Associate Fellow, Associate and 
Graduate. 

* Of the ordinary members of the Council, that number 
necessary to create seven vacancies shall retire annually. 
The retiring members shall be those with the longest 
service since their last election but they shall be eligible 
for re-election. 


* Nominations of candidates for election to the Council 
must be received by the Secretary not later than /0th April 
in each year and shall include statements in writing by the 
candidates that they are willing to serve. The nomination 
forms shall be signed by one proposer and two seconders, 
all of whom shall be voters.” 


Nomination forms may be obtained on application to 
the Secretary. 


S1xtH Louis BLERIOT LECTURE 


The Sixth Louis Bleriot Lecture will be given by General 
Guy du Merle, F.R.Ae.S., on “Making Commercial Aircraft 
Pay.” The Lecture will be held at the Institution of 
Mechanical Engineers, Storey’s Gate, St. James’s Park on 
Thursday 19th March 1953 at 6 p.m. 


JOINT CONFERENCE OF ASLIB AND NATO 


A Joint Conference of the Aslib Aeronautical Group and 
the Advisory Group for Aeronautical Research and 
Development (A.G.A.R.D.) of the North Atlantic Treaty 
Organisation will be held from Friday 27th March to 
Tuesday 3lst March 1953 at the College of Aeronautics, 
Cranfield. Among the subjects to be discussed are the 
Cataloguing of Aerodynamic Data (with special reference 
to the Netherlands plan set out in N.L.L. Report No. 116) 
and a common international system for indexing, 
abstracting and cataloguing all aeronautical reports. 


The Conference fee will be £1 5s. Od. (including preprints 
memberships of the College club and transport from 
Bedford railway station) and the charge for accommodation 
will be £1 Ss. Od. a day. Full particulars may be obtained 
from Mr. C. W. Cleverdon, College of Aeronautics, 
Cranfield, Bletchley, Bucks. 
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MEMBERS’ NEW APPOINTMENTS 


H. D. Burke (Associate), formerly Technical Sales 
Representative of the Palmer Tyre Ltd., has recently been 
appointed Manager, General Aircraft Division, Palmer 
Aero Sales Department. 

Fur./Lr. W. R. Carus (Associate) has recently been 
appointed a member of the flying instructional staff at the 
Royal Air Force College, Cranwell. He was formerly 
Personal Staff Officer to the Commandant of the Royal Air 
Force Flying College, Manby. 

DouGLas S. CAUSER (Associate Fellow), formerly with 
British Overseas Airways, has been appointed a partner of 
Causer and Company, Drawing Office Specialists, London, 
W.12. 

B. F. DocKkREE (Associate Fellow) has been appointed 
Chief Designer of Delaney Gallay Ltd. 

G. R. Epwarpbs (Fellow) has recently been appointed a 
Director of Vickers-Armstrongs Ltd.,and General Manager 
and Chief Engineer, Aircraft Division. 

H. H. GARDNER (Fellow) has been appointed Chief 
Designer, Guided Weapons, of Vickers-Armstrongs Ltd., 
Weybridge Works. 

V. A. HicGs (Associate Fellow) has recently been 
appointed Manager of the Aero and Auto Equipment Sales 
Department of the British Thomson-Houston Company. 

D. KeitH Lucas (Fellow), Chief Designer, Short Bros. 
and Harland Ltd., has been made a Director of the 
Company. 

T. F. C. LAWRENCE (Associate Fellow) has recently been 
appointed to the High Speed Aerodynamics Laboratory of 
the Department of Supply at Adelaide. 

F. K. Main (Graduate) has been appointed Chief 
Draughtsman with the British Ermeto Corporation Ltd., of 
Maidenhead. 

P. R. OWEN (Associate Fellow), formerly at R.A.E., has 
recently taken up the appointment of Director of the Fluid 
Motion Laboratory, Manchester University. 

A. A. ROSE (Associate Fellow), formerly with Saunders- 
Roe Ltd., has taken up an appointment on the design staff 
of A. V. Roe (Canada) Ltd. 

B. STEPHENSON (Associate Fellow) has been appointed 
Chief Designer Aircraft of Vickers-Armstrongs Ltd., 
Weybridge Works. 

PROFESSOR PHRIXOS J. THEODORIDES (Associate Fellow) 
has taken up a new appointment as Research Professor at 
the Institute for Fluid Dynamics and Applied Mathe- 
matics, University of Maryland. 

J. A. C. Wititiams (Associate Fellow), formerly with 
Guthrie and Co. Ltd., of Singapore, has _ recently 
been appointed Senior Lecturer in Management at the 
Wolverhampton and Staffordshire Technical College. 


ELLiotTt MEMORIAL PRIZE 
The Elliott Memorial Prize has been awarded to 
Corporal Aircraft Apprentice P. E. Newton, of the May 
1950 Entry at R.A.F., Halton, who has obtained the 
highest marks in the General Studies Examination. The 
prize will be presented at the Graduation Prize Giving on 
March 1953. 


ACKNOWLEDGMENTS 


The Council acknowledge with thanks the return of back 
numbers of THE JouRNAL from J. R. Bradford, Esq., 
Associate; N. A. Champion, Esq., Associate Fellow: and 
A. J. L. Mitchell, Esq., Associate Fellow. 


SOCIETY FOR EXPERIMENTAL STRESS ANALYSIS 


The Spring meeting of the Society for Experimental 
Stress Analysis will be held at the Hotel Schroeder, 
Milwaukee, Wisconsin, from 20th to 22nd May 1953. 
Particulars may be obtained from Mr. W. F. Hofmeister, 
Metallurgical Laboratory, Chain Belt Company, 4501 West 
Greenfield Avenue, Milwaukee, Wisconsin. 


ROYAL AERONAUTICAL “SOCIETY — -NOTICES 


ASSOCIATE FELLOWSHIP EXAMINATION RESULTS 
DECEMBER 1952 


The following were successful in the Associate Fellowship 
Examination held in December 1952 :— 


Part 


D. G. A. Bridgnell, Properties of Matter, Heat, Light and 
Sound. 
D. R. Brodie, Aerodynamics I. 
K. G. Draper, Properties of Matter, Heat, Light and Sound. 
J. W. Drife. Pure Mathematics; Properties of Matter, Heat, 
Light and Sound; Strength of Materials and Theory of 
Structures I. 
D. J. Frater, Properties of Matter, Heat, Light and Sound. 
D. R. Henbest, Properties of Matter, Heat, Light and Sound. 
J. L. Johnston, Properties of Matter, Heat, Light and Sound; 
Aerodynamics 1; Strength of Materials and Theory of 
Structures I. 
T. C. Lee, Applied Mathematics; Properties of Matter, Heat, 
Light and Sound; Theory of Machines. 
J. MacNaughton, Properties of Matter, Heat, Light and 
Sound; Strength of Materials and Theory of Structures I; 
Aerodynamics I. 
K. J. Neeves, Applied Mathematics; Magnetism and 
Electricity; Strength of Materials and Theory of Structures I; 
Theory of Machines. 
K. G. Page, Properties of Matter, Heat, Light and Sound. 
N. R. Reddy. Strength of Materials and Theory of 
Structures I. 
D. J. Richardson, Aerodynamics I. 
R. D. Saha, Strength of Materials and Theory of Structures I; 
Theory of Machines. 
B. Singh, Strength of Materials and Theory of Structures I. 
R. P. Stanford, Pure Mathematics; Aerodynamics I; Strength 
of Materials and Theory of Structures I. 
I. K. Trundley, Pure Mathematics; Properties of Matter, 
Heat, Light and Sound; Aerodynamics I; Strength of Materials 
and Theory of Structures I. 
P. Shu-Yan Wong, Applied Mathematics; Properties of 
Matter, Heat, Light and Sound; Strength of Materials and 
Theory of Structures I; Theory of Machines. 


Part I—ABROAD 


W. K. Aiyadurai, Melbourne, Australia, Pure Mathematics: 
Properties of Matter, Heat, Light and Sound; Strength of 
Materials and Theory of Structures I; Theory of Machines. 

K. S. Chandrashekaran, Bangalore 4, India, Aerodynamics 1; 
Strength of Materials and Theory of Structures I. 

S. Singh, Delhi, India, Strength of Materials and Theory of 
Structures 1; Theory of Prime Movers and Fuels. 

A. Subrahamanyan, Bangalore, India, Properties of Matter. 
Heat, Light and Sound; Theory of Prime Movers and Fuels. 


Part II—Home 


P. J. Farmer, Aircraft Materials. 
D. W. Pickston, Aircraft Materials; Engineering Economics. 
C. R. Raynham, Strength of Materials and Theory of 
Structures II. 
D. Rowe, Strength of Materials and Theory of Structures II. 
P. Simpson, Aircraft Design. 
I. H. Smart, Strength of Materials and Theory of Structures II. 
F/Lt. G. Wood. Air Transport: Aircraft Materials. 


ParT II—PROVINCIAL 


A. G. Creak, Cirencester, Glos., Aerodynamics II. 
W. J. Hanlon, Belfast, N. Ireland, Strength of Materials and 
Theory of Structures II. 
W. F. Meichan, R.A.F. Aldergrove, N. Ireland, Meteorology; 
Navigation. 


Part II—ABROAD 


A. F. El-Kerdani, Heliopolis, Egypt, Aerodynamics II. 
P. Krishna Iyer, Chromepet, Madras, India, Engine Design. 
R. B. Lal, S.P.O.. Madras, India, Engine Design. 
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Diary 


LONDON 


March 
GRADUATES’ AND STUDENTS’ SECTION, Annual General 
Meeting. 4 Hamilton Place, W.1. 7.30 p.m. Followed 
by an aeronautical film show. including the ~S.B.A.C. 
Show 1952." Coffee 7-7.30 p.m. 


March 12th 
MAIN LECTURE at DERBY. Some Aspects of High Per- 
formance Jet Aircraft. Sqdn. Ldr. W. A. Waterton, G.M.. 
A.F.C. Rolls-Royce Welfare Hall. Derby. 6.15 p.m. 


March 17th 
GRADUATES’ AND STUDENTS’ Section. Visits to the “ Daily 
Mail” and Daily Herald printing works. 


March 19th 
SixtH Louis”) BLERIOT LECTURE. Making Commercial 
Aircraft Pay. General Guy du Merle. Institute of 
Mechanical Engineers. Storey’s Gate, S.W.1. 6 p.m. 
(Tea 5.30 p.m.) 


March 25th 
GRADUATES’ AND STUDENTS’ SECTION. Convertible Aircraft. 
J. Shapiro. 4 Hamilton Place, W.1. 7.30 p.m. 


March 27th 
FuLL Day Discussion ON FATIGUE. At the Chemistry 
Lecture Theatre, University College. Gower Street, W.C.1. 
10 a.m, to 6.30 p.m. 


March 3lst 
SECTION LECTURE. Experimental Work on Boundary Layer 
Flow. W. E. Gray. In the Library. 4 Hamilton Place, 
pam. 


April 22nd 
GRADUATES’ AND STUDENTS’ SECTION. Aircraft Accident 
Investigation. H. Caplan. In the Library, 4 Hamilton 
Place, W.1. 7.30 p.m. 


April 23rd 
MalINn LECTURE at GLASGOW. Civil Jet Operations. Captain 
A. M. Majendie. The Royal Technical College. Glasgow. 
7.30 p.m. 


BRANCHES 


March 9th 


Halton—Films. 
6.45 p.m. 


Branch Hut, R.A.F. Station, Halton. 


March 11th 
Brough—Airships. Lord Ventry. Electricity Showrooms, 
Ferensway. Hull. 7.30 p.m. 
Halton—Visit to Vickers-Supermarine, Winchester. 
Hatfield—Navigational Problems Associated with the 
Introduction of Comet Aircraft. Captain A. M. Majendie. 
Senior Staff Mess, de Havilland Aircraft Co. Ltd.. Hatfield. 
6.15 p.m. 
Southampton—Can We Reach Mars? E. G. Burgess. 
Institute of Education. University of Southampton. 7 p.m. 


March 12th 
Derby—Main Society Lecture. Some Aspects of High 
Performance Jet Aircraft. Sqdn. Ldr. W. A. Waterton, 
G.M., A.F.C. Rolls-Royce Welfare Hall. 6.15 p.m. 
Isle of Wight—New Materials and Types of Aircraft 
Structure. H. J. Pollard. Club House. Saunders-Roe 
Sports and Social Club. Church Path, E. Cowes. 6 p.m. 
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GIETY 


March 12th 


Portsmouth—Films. Technical School of the de Havilland 
Aircraft Co. (Portsmouth) Ltd.. The Airport. Portsmouth, 
6 p.m. 


March 16th 
Halton—Branch Night. 
Halton. 6.45 p.m. 
Henlow—Helicopter Design Problems. A. H. Yates. 163 
Building, R.A.F. Technical College. 7.30 p.m. 


Branch Hut. R.A.F. Station, 


March 18th 
Coventry—Some Aspects of Airline Flying from the 
Pilot's Point of View. Captain R. Rymer, B.E.A. Wine 
Lodge. 7.30 p.m. 


Leicester—Helicopter Design Problems. A. H. Yates, 
Lecture Theatre. College of Technology, Leicester. 
TAS 


Reading and District—Flight of Birds. Captain J. L. 
Pritchard, C.B.E. Abbey Gateway Rooms, Abotts Walk, 
Reading. 7.45 p.m. 


Weybridge—Making Aircraft Films. P. D. Normanville. 
Vickers-Armstrongs Ltd.. Weybridge Works. 6 p.m. 


March 19th 


Glasgow—Engine Development Flight Testing. Wing 
Commander T. H. Heyworth. Rolls-Royce, Hillington. 
7.30 p.m. 


Southampton—Supper Dance. 
March 2Ist 


Birmingham— Annual 
Birmingham. 


Dinner. White Horse Hotel, 


March 23rd 
Bristol—Junior Members’ Papers Competition. Conference 
Room, Filton House. Bristol Aeroplane Co. Ltd. 6 p.m. 


Halton—Junior Members’ Night. Branch Hut, R.A.F. 
Station, Halton. 6.45 p.m. 


March 25th 
Hatfield Discussion evening for Members’ Papers. Senior : 
Staff Mess, de Havilland Aircraft Co. Ltd., Hatfield. 
6.15 p.m. 

Preston—-Some Modern Aspects of Production. Professor 
J. V. Connolly. The Assembly Hall, Technical College, 
Corporation Street. 7 p.m. 


March 26th 


Portsmouth—Annual General Meeting. Technical School § 
of the de Havilland Aircraft Co. (Portsmouth) Ltd., The 
Airport. 6 p.m. 


March 3lst 


Belfast—Rocket Propulsion and Interplanetary Flight. 
A. V. Cleaver. New Hall. Kensington Hotel, College 
Square East, Belfast. 7 p.m. 

Henlow—Aerodynamics of Wind Tunnel Design. J. C. 
Gibbings, 163 Building, R.A.F. Technical College, Hen- 
low. 7.30 p.m. 


April Ist 
Chester—Fire Prevention in Aircraft. J. White. 


Grosvenor Hotel, Chester. 7.30 p.m. 
April 8th 
Weybridge—The Work of the Estimating, Planning and 


Progress Departments. 
and P. D. Imlach. 
Works. 6 p.m. 


J. F. Hawkins. W. R. Coomber 
Vickers-Armstrongs Ltd., Weybridge 


Om! MARCH 1953 
| 
agar 
for 
rap 
| 
cou 
: 
but 
(fo 
int 
am 
: 
ali 
do 
clo 
col 
air 
| 
IS 


illand 
10uth. 


lation, 


163 


n the 
Wine 


Y ates. 
cester, 


Wing 
ngton. 


Hotel, 


erence 
6 p.m. 
R.A.F. 


»fessor 
ollege, 


School | 


The 


Flight. 


“ollege 


White. 


ng and 
yomber 
y bridge 


April 13th 
Derby—Electrics. H. West. Rolls-Royce Welfare Hall, 
Nightingale Road, Derby. 6.15 p.m. 


Glasgow—Annual General Meeting. St. Enoch Hotel. 
7.30 p.m. 


April 14th 
Bristol—Annual General Meeting and Film Show. Con- 


ference Room, Filton House, Bristol Aeroplane Co. Ltd. 
6 p.m. 


April 15th 
Hatfield—Fatigue: What It Is and Ways of Reducing its 
Incidence. Major P. L. Teed. Senior Staff Mess. de 
Havilland Aircraft Co. Ltd.. Hatfield. 
Preston—Annual Genera! Meeting. followed by a Brains 
Trust. Assembly Hall of the Technical College, Corpora- 
tion Street, Preston. 7 p.m. 
Reading and District—Meteorological Service for B.O.A.C. 
Comet Operation. E. Chambers. The Abbey Gateway 
Room, Abotts Walk, Reading. 7.45 p.m. 
Southampton—Aicrcraft Maintenance. B. S. Shenstone. 
Institute of Education, University of Southampton. 7 p.m. 


ROYAL AERONAUTICAL SOCIETY--NOTICES XIX 


April 20th 
Halton—Branch Night. 
Halton. 6.45 p.m. 


Branch Hut, R.A.F. Station, 


April 
Belfast—Annual General Meeting. New Hall, Kensing- 
ton Hotel, College Square East, Belfast. 7 p.m. 


April 22nd 
Coventry—Annual General Meeting and Films on Lubri- 
cation, by courtesy of Esso. Wine Lodge. 7.30 p.m. 


April 23rd 
Glasgow——-Main Lecture—Civil Jet Operations. Captain 
A. M. Majendie. Royal Technical College, Glasgow. 
7.30 p.m. 


April 27th 
Halton—Bird Flight. Captain J. L. Pritchard, C.B.E. 
Branch Hut. R.A.F. Station, Halton. 6.45 p.m. 


Henlow—Student Members’ Premium Lectures. 
Technical College. Henlow. 7.30 p.m. 


R.A.F. 


THE NEW ZEALAND DIVISION 


Members may be interested in the following letter to the Secretary from the Honorary 
Secretary of the New Zealand Division of the Society. 


Dear Dr. Ballantyne :— 

Our report and financial statement will be coming 
forward early in the New Year and that will detail the 
facts of our activities for the year. Although the growth 
of the Division is not shown by the figures as being very 
rapid. there is an awakening interest in the doings of the 
Society which seems to be fulfilling a real need in the 


country. Aviation activities, though not very profound in 


this country, are on a surprisingly wide basis. We cannot 
muster large numbers of experts in aeronautical matters, 
but we do manage to gather together very large numbers 
(for the size of the country) of people who are extremely 
interested in aeronautical matters in a general or even an 


Senior } amateur way. A large part of this class of person is young 


atfield. | 


and this is an extremely promising sign to find that the 
activities of the Society are catering for the rising gener- 
ation. 

The problems of Aviation in this country are without 
doubt very'different from the United Kingdom, being more 
closely connected with the operating side than with the 
manufacturing. The development of aerial topdressing in 
this country has followed a line more or less unique and 
is certainly becoming a really significant feature in the 
economy of the country. In actual fact New Zealand 
people are extremely air minded. Our internal airways 
communications are becoming more and more popular. 
The awkward geographical obstacles, such as the division 
of the Islands and the mountainous country, are such that 
air travel brings conveniences and advantages which do 
not apply to other forms of transport. These advantages 
have meant that the ordinary member of the population 
is far more likely to be directly interested in Aviation 
matters than would be the case in many other countries. 
Aerial topdressing again has brought directly to the eye 
of the farmer something of the potentialities of air trans- 
port and possibly farmers as a class are far less likely to 
be enamoured with air transport than almost any other 
class of person. Since New Zealand is so primarily a 
farming country it is extremely significant that air trans- 
port has in this way come to the attention of the farmer. 


16th December 1952. 


The expansion of aeronautical activities in this country 
at the present time is rapid and shows every sign of 
continuing to be rapid. It is satisfactory then to find that 
the Royal Aeronautical Society is taking an active and in 
many respects a leading part in this rapid development. 

Another interesting item in N.Z. Aviation is the award 
to Tasman Empire Airways of the Cumberbatch Trophy 
for safety record. Although neither the Society nor the 
Division have any direct connection with this Trophy, 
nevertheless it is indicative of the representation of the 
Division that T.E.A.L. has many of its staff active members 
of the Division. These include the T.E.A.L. General 
Manager as Vice-President, T.E.A.L. Chief Engineer as a 
past Chairman of Auckland Branch and other members 
of T.E.A.L. take an active part as office bearers. 


In the coming year New Zealand will be definitely on 
the map with the England to Christchurch Air Race. 
Although the N.Z. Division is not directly concerned in 
this, nevertheless again it is indicative of the Division’s 
ubiquity that several members are directly connected with 
the organisation of the Race at this end. 

Quite apart from these direct activities, there is no 
doubt that the regular meeting of members of the Division 
one with another is of inestimable benefit. A community 
of interest and an active understanding of the other person’s 
problems are not the least of these. The opportunity for 
making and maintaining friendships is also something 
which the Division supplies of the very greatest value. 
The active interest of yourself and. in fact of the whole 
body of the parent Society, is just another aspect of this 
and this in itself is ample reason enough to send you our 
very sincere and best wishes for Christmas and the New 
Year. 


Yours faithfully, 
T. M. COLERIDGE, 
Honorary Secretary, 
New Zealand Division of the Royal Aeronautical Society. 
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ELECTIONS 


The following is a list of new members and transfers 
of membership of the Society :— 


Associate Fellows 
John Alfred Beadon 
Terence Beville Adair 
Boughton 
Ronald Claude Graham 
Boxer (from Graduate) 
Ralph Brown 
(from Associate) 
Eric Leslie Carlisle 
(from Associate) 
Eleanor Lettice Curtis 
(from Associate) 
Stanley Hargreaves Ellis 
(from Graduate) 
Geoffrey Fawn 
Alan James Greenhalgh 
(from Graduate) 
Stephen John Hancock 
(from Graduate) 
Wm. Marshall Hargreaves 
Ronald Frederick Hill 
(from Associate) 
Edward Henry Holder 
(from Associate) 


Associates 

Cecil Frank Allen 

Dennis Arthur Thomas Ball 

Gysbert Bloemendaal 

Roy Milner Bradley 

Leslie Douglas Chapman 

James Henry Cruddace 

Albert Ernest Dennard 

Walter Miller Fitch 

Eric Arthur Franks 

Neil McLeod Gibson 
(from Companion) 


Graduates 
Richard Cameron Aitken 
John Chas, Arthur Baldock 
(from Student) 
Robin John Balmer 
John Robert Burton 
(from Student) 
Antony Shui-Lau Chan 
(from Student) 
Francis Charles Drinkwater 
Ronald Arthur Dutton 
Edward Vernon Griffith 
Jack Womersley Headley 
(from Student) 
Norman Jackson 
(from Student) 
Morris Vivian Jenkins 
Herbert Judson 
John Irvine Lattey 
Harold Moore Lofthouse 
(from Student) 
William Robert McLean 


Students 

Mond Khan Amir 

Ronald Frederick Ayers 

James Sidney Baker 

Martin Charles Bennett 

John Arthur Boddy 

Thomas’ Patrick Joseph 
Walter d’Alton Bourke 

Derek William Brown 

Derek George Burton 

Frank Richard Dugay 

George Kasin Ferris 

Edgar Ronald Fraser 

Ronald Anderson Gellatly 

Dennis Gortmans 

Geoffrey Leonard Hearne 

Reginald Arthur Jackson 

Edward Henry James 

Brian Greig Jolly 


Thomas George Kent 
(from Graduate) 
Bernard Lovegrove 
(from Graduate) 
William Muirhead 
Peter Jack Palmer 
(from Graduate) 
Bertram Anthony Peaster 
Theordore Stanley Revell 
(from Associate) 
Colin Arthur Richardson 
(from Associate) 
Alan Scrimshaw 
Harold Warren Shattock 
(from Graduate) 
Richard Herbert Taylor 
(from Graduate) 
Walter Edward West 
Mark Vincent Wignall 
Geoffrey Edward Williams 
(from Graduate) 


Derek Norman James 
(from Companion) 

Bernard John Naish 

John Leslie Nunn 

James Ronald Plenderleith 

Clifford Edwin Refoy 

David Royle 

Alexander Rodger Sturgeon 

Thomas Benjamin Taylor 

Geoffrey Thomas White 

Arthur Wilby 


David Needham 
(from Student) 

Denys Randolph 

Peter Marcus Ray 

Robert Chas. Leslie Sample 
(from Student) 

Harold George Sevier 
(from Student) 

Harold Alexander Simon 
{from Student) 

Vivek Ranjan Sinha 

Frederick Thomas Smith 

Christopher John Hamshaw 
Thomas (from Student) 

Ernest Cyril Thomasson 

Denis Edward Weeding 
(from Student) 

Rex Maddock Whalley 
(from Student) 

Frank Gerrie Willox 

Norman James Wilson 


Dennis King 

Brian Richard Luesley 

John David Francis Mac- 
Naughton 

Afzal Haq Malik 

William Colin Frank Moxon 

Gordon Oulsnam 

Anthony John Pape 

Brian William Payne 

John Lloyd Hendy Phillips 

Frederick Bowen Powell 

John Graham Price 

John Holdom Stock 

Michael Edwin Thomas 

William Malcolm Thomas 

Dennis Albert Vaughan 

Wallace Barrie Dannatt 
Wardle 


Companions 
Biaal Chandra Dutta 
Paul Leon Albert Raes 
(from Student) 


Alan Charles Keith 
Robertson 


ANNUAL SUBSCRIPTIONS 


Members are reminded that their annual subscriptions 
became due on Ist January 1953. The rates are:— 


HOME ABROAD 

Associate Fellows .. 44 0 
* Associates 3 0 -3 10 
Graduates (aged under 26) 2; 210 2 2 6 
Graduates (aged 26 and over 2 26 22 
Students (aged under 21) .. i £40 1 1 0 
Students (aged 21 and over) 36 
Companions 3 3 10 
Founder Members .. 2 2 219 


*Any Associate elected before Ist October 1947 may, 
if he wishes, elect not to receive the JOURNAL, and in this 
case his subscription will be reduced by £1 Is. Od. to 
£2 2s. Od. 

It will avoid delay and confusion if members, when 
sending remittances for subscriptions, will state their names 
clearly and give their addresses and grades of membership. 
Remittances should be made payable to the Royal 
Aeronautical Society. 


CHANGES OF ADDRESS 
To assist in keeping the records of members correct and 
up to date the Secretary will be glad if all members will 
notify him as soon as possible of changes of address. 
When notifying changes please give the following 
particulars : — 
Name (in block letters). 
Grade of membership. 
New address (in block letters). 
Old address. 
Changes of address must be received before the [5th 
of the month in order to be effective for the JOURNAL 
for the following month. 


JOURNAL BINDING 
Self-Binder Cases 

Self-Binder cases of the “ Easibind ” type are available 
from the offices of the Society. These binders are for 
members who do not have their Journals permanently 
bound, or who wish to keep their Journals together during 
the year for binding later. 

These cases will hold 12 Journals which are kept in 
place by means of flexible steel wires. Journals can be 
inserted or withdrawn easily without damage, so preserving 
the contents for permanent binding later. 
will open flat at any page. 


The binder is strongly made in durable dark blue leather | 
cloth on stiff board covers and has gold lettering on the | 


spine. The year is not blocked on the spine but there is 
a panel on which members who wish to use the binder 
as a permanent case can put the date. 

The cost is lls. 6d. each including postage and packing 
for either the size to fit 1952 and previous Journals, or 


for the size to fit the Journal from January 1953, which | 
Orders and remittances should ; 


has been increased in size. 
be sent direct to the Secretary at the Offices of the Society 
and it is important to state whether the old size or new 
size is required. 


Permanent Binding 
There is no increase in the price of permanent binding 
of Journals. The prices are :— 
1952 Volume (including packing and postage 16s. Od. 
Previous Volumes (including packing and postage) 18s. 0d. 
Journals, with a note of the name and address of the 
sender, should be sent direct to the Lewes Press, Friars 
Walk, Lewes, Sussex, and the remittance to the Secretary 
at the Offices of the Society. 


The Journals | 
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The Scope and Limitations of the 
Photoelastic Method of Stress Analysis 


H. T. JESSOP, T.D., M.Sc., F.lnst.P. 


The 86lst Lecture to be given before the Royal Aeronautical Society was held at the 
Institution of Mechanical Engineers, Storey’s Gate, London, S.W.1, on 27th November 
1952. Mr. G. R. Edwards, C.B.E., F.R.Ae.S., Vice-President of the Society, presided and 
introduced the lecturer, Colonel H. T. Jessop, T.D., M.Sc., F.Inst.P., Reader in Photo- 
elasticity at University College, London, who was a pupil of Professor Filon and later 
worked as Professor Filon’s Research Assistant Mr. Edwards said that his company 
had been fortunate in making use of Colonel Jessop’s facilities for some photoelastic 
analyses of their own pieces of machinery: the photoelastic method of analysis was one 
which the Aircraft Industry was finding of increasing value. 

During the lecture Colonel Jessop demonstrated the stress lines in various shapes of 

materials with a polariscope. 


Introduction 


The problem of determining the distribution of stress 
over the various components of a structure such as an 
aircraft frame is one which has always presented 
considerable difficulty. For a large proportion of the 
components calculations of the stresses can only be 
carried out on the basis of assumptions and approxima- 
tions, the effects of which cannot be estimated with any 
degree of accuracy and, since considerations of size and 
weight impose limitations upon the employment of a 
large safety factor, there has always existed a need for 
some practical means of checking the results of calcula- 
tions or of obtaining more reliable estimates of the 
stresses. 

One such means is to be found in the photoelastic 
method of stress-exploration. In comparison with other 
methods it is rapid and inexpensive, in many cases it 
gives much fuller and more exact information, and it is 
moreover the only method by which the magnitudes of 
the stresses in highly localised stress-concentrations can 
be determined with any degree of accuracy. This last 
point is of particular importance in modern aircraft 
design, for it is becoming increasingly urgent to have a 
more exact knowledge of the magnitudes of these local 
Stresses in order to avoid the possibility of their leading 
to fatigue failures. The photoelastic method is also the 
only one which gives any information as to the stress- 
distribution at internal points, and while it is usually the 
boundary stresses which are of chief interest to the 
engineer, there are cases in which the rate at which 
the stresses fall off below the surface of the material is 
likely to have a significant effect upon the possibility of 
failure. 


by 


Outlines of the photoelastic method have appeared in 
various publications''’, and the general principles of the 
method are too well known to need recapitulation here, 
but many developments in materials and techniques have 
taken place in recent years and these have increased the 
power and range of the method to an extent which is 
not yet generally appreciated by the engineers for whose 
use it has been designed. The new developments have 
taken place largely in research laboratories which have 
not the resources to make extensive investigations, but 
the methods have been applied to enough practical 
examples to demonstrate their value in a wide variety of 
stress-problems. In this paper an attempt is made to 
give an account of the present scope of the method and 
of the techniques employed, and to indicate some of the 
types of problem in which it could be used with 
advantage. 


Two-Dimensional Techniques 


The two-dimensional method is substantially that 
developed by Coker and Filon™? in the years 1912-1930, 
although the introduction of new materials and new 
techniques has led to simplified procedure and improved 
accuracy. This method is strictly applicable only to 
components whose form is that of a flat plate of uniform 
small thickness to which the loads are all applied at the 
edges and all act in the plane of the plate. For such 
components the photoelastic model is generally 
examined under load in the ordinary polariscope, and 
the complete stress-distribution is obtained in most cases 
rapidly and with considerable accuracy. Measurements 
of the fringe-pattern seen in the polariscope yield 
directly (1) the magnitude of the stress at all points on 
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Ficure 1. Fringe-pattern showing the stress-distribution around 

an array of holes in a flat plate under tension, and demonstrat- 

ing that the peak tensions are approximately the same in the 
wide and narrow ligaments. 


the free boundary, and (2) the magnitude and direction 
of the greatest shear-stress at all internal points. In 
many practical problems this information is all that is 
required to check the suitability of a design or to deter- 
mine the safe load for a component. and often a mere 
inspection of the fringe-pattern will suffice. If the 
magnitudes of the separate stresses at any point are 
required, these may be found by simple calculations 
based on measurements of the fringe-pattern. 

Figure 1 is a photograph of the fringe-pattern in a 
plane model which was one of those used to investigate 
the stresses round the steam-pipe holes in a boiler drum. 
The diameter of the drum was very large compared with 
the thickness of the shell and a flat plate under uniform 
tensile load reproduced, with good accuracy, the 
conditions of hoop tension in the shell. The investigation 
was undertaken to find the effect of uneven spacing of the 
holes, and the main interest was in the strains around 
the holes and the possibility of these being so large in the 
narrower ligaments as to cause leakage under high 
pressure. The fringe-patterns in loaded models showed 
at once that there was little difference between the peak 
tensions in the wide and narrow ligaments: the fringe- 
orders at opposite sides of the holes in Fig. | are 
obviously approximately the same. Measurements and 
calculation of the actual stresses showed that for a wide 
range of differences in spacing the difference in peak 
stresses did not exceed 5 per cent. 

Figures 2 and 3 illustrate a bolted or riveted cleat 
joint, a type of problem which might well be investigated 
by photoelastic methods. In the model shown, the 
cleat is made of Perspex, the photoelastic sensitivity of 
which is very small, so that the pattern observed through 
the whole model is sensibly that due to the stresses in the 
middle bar. This gives a qualitative comparison of the 
stresses in such a bar with the two different bolt spacings 
shown, but even this is sufficient to show that with the 
closer spacing the bolt loads are more equally 
distributed, while the peak shear-stress in the bar is 
reduced. 

A full and accurate investigation of any specified 
pattern of such a joint could be made by using either a 
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reflection method on loaded models, or the frozen stress 
method. The stresses in the cleat and the bar could 
then be completely separated and the effects of varying 
both spacing and size of bolts explored. This preliminary 
test suggests that there will be in any particular design 
an optimum spacing, and that this spacing will be closer 
than that normally employed in such joints. 

The two-dimensional technique may also be 
employed in many cases in which the stress-distribution 
is not strictly two-dimensional, such as an I-section spar 
in which the loading is in the mid-plane of the spar and 
parallel to the plane of the web. Here the stresses in 
the web will be sensibly two-dimensional and may be 
deduced with good accuracy from the fringe-pattern, 


= 
(1) 
FIGURE 3. 
FiGuRe 2. 


Figure 2. Simple model representing a bar attached to a cleat 
by a row of pins. 


FIGURE 3, Fringe-patterns in the bars of two models similar to 
that in Fig. 2 under tension. The patterns show on inspection 
(a) that the mean stress across each of the corresponding 
sections AB, CD is greater in the model with more closely 
spaced rivets, and therefore that the load is more equally 
shared by these rivets, and (4) that the greatest shear-stress in 
each bar occurs at the boundary of the first hole, and is less in 
the model with the closer spacing. The patterns around the 
lower holes are not significant, for here the higher stresses in 
the Perspex cleat begin to interfere with the smaller ones in 
the more optically sensitive bar. 
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while the pattern in the flanges will allow the estimation 
of the mean stresses through the thickness of the flange. 
The stresses in the fillets between flange and web would 
need a full three-dimensional exploration. In tests made 
recently on such a spar model by Vickers-Armstrongs 
Lid.. the photoelastic investigation revealed high stresses 
which their calculations had failed to indicate, and 
sugvested the type of modification which should be made 
in the design. A second test, on a model built to the 
modified design, showed a satisfactory reduction in these 
stresses. In this case photoelastic tests occupying three 
weeks yielded better information than was obtained 
from calculations which had taken a much longer time. 


Three-Dimensional Techniques 
TH! FROZEN STRESS METHOD 

The general characteristics of the frozen-stress effect 
are Well known to all interested in practical stress 
analysis. The features of the effect which are relied 
upon in stress investigations are first, that the “ frozen ~ 
photoelastic fringe-pattern shall correspond with the 
stress distribution in the model under direct load, and 
second, that no change in the “ frozen” stress shall be 
produced by cutting the model. That these two 
conditions are, in fact, fulfilled in the two-dimensional 
case and in simple three-dimensional cases, has been 
demonstrated on a variety of models. The photographs 
of Figs. 4 and 5 illustrate these two features of the 
“frozen” stress. The method renders it possible to 
freeze the stresses into a model of any shape and, by 
cutting thin slices and examining them photoelastically, 
to explore the stresses at all points of the model. 


Three main methods of examination are available : — 


(1) A thin slice of the model may be examined in the 
polariscope using the ordinary two-dimensional tech- 
nique, with the light passing through the slice either 
normally or obliquely (Figs. 6 and 7). 


(2) A narrow pencil of light may be passed through 
the uncut model (or any part of it) and the stresses at 
points on the path of the pencil explored by measure- 
ment of the fringes which are observed by scattered 
light’ (Fig. 8). 

(Since this method does not involve slicing the model 
it may equally well be employed on a model under load, 
and in certain cases it will give all the information 
required without necessitating the freezing of the model.) 


(3) A thin slice of the model may be examined in 
polarised light on a universal tilting stage’. By making 
a systematic series of small displacements about two 


FIGURE 8. 
along the line A B of the slice from the grooved cylinder shown 


Fringes seen by scattered light in a beam passing 


in Fig. 6, The spacing of the fringes gives a measure of the 
difference between the axial and circumferential stresses. 


PHOTOELASTIC METHODS OF STRESS ANALYSIS. 


FicurE 4, Fringe-patterns in two similar models, one (B) under 

direct load, the other (A) having the stresses frozen into it 

under a similar load. The stresses were frozen into model A 
eighteen months before this photograph was taken. 


FiGuRE 5. FIGURE 6. 
FiGurRE 5. Frozen fringe-pattern in a model, showing that the 
pattern is not disturbed by cutting the model. 
FIGURE 6. Fringe-pattern in a slice cut in a diametral plane 
from a grooved cylinder in axial compression. The fringes 
give the difference of the principal stresses in the plane of 
the slice. 
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FIGURE 7. Fringe-pattern in a slice cut from the web of a crank- 
shaft. The shaded portion in the diagram indicates the slice. 
The fringe-pattern shows clearly the distribution of shear 

stress around the journal. 
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Ficure 9. Fringe-patterns showing the stress concentrations in 


fillets of different radii. In this example doubling the fillet 
radius has reduced the peak stress from 6 units to 5. 


perpendicular axes the directions in space of the three 
principal axes at any chosen point of the slice may be 
determined. 

In most cases the quickest and most accurate 
approach will be by a combination of two, or all, of these 
methods of examination and will lead to the determina- 
tion of the principal stress differences. As in the two- 
dimensional case the separation of the stresses will, in 
general, require an integration process. 

The time and labour involved in a three-dimensional 
investigation vary enormously with the type of problem. 
A complete exploration of the stresses at all points of 
even a simple model would be a very long process, 
probably involving the freezing and slicing of two or 
three models, but for practical purposes the engineer is 
interested mainly in the stresses in certain specified 
regions, and in most cases these happen to be the regions 
in which the solution is comparatively simple. In order 
of simplicity regions may be classified as :— 


(a) Free boundaries in a plane of symmetry. Here 
only the magnitudes of two stresses have to be found, 
and this may be done most easily by a combination of 
methods | and 2, or in some cases by two observations 
by method 2 without slicing the model. ° 


(b) Free boundary not in a plane of symmetry. This 
involves finding the directions of the two stresses, as 
well as their magnitudes, but the solution is very little 
more complex than for case (a). 


(c) Interior regions in a plane of symmetry. Here 
the direction of one stress is known. The directions of 
the stresses in the plane of symmetry and the magnitude 
of the stress difference are obtained by method |. The 
second stress difference is obtained by combining 
observations by method 2 with those of method 1. A 
comparatively simple integration leads to the values of 
the separate stresses’. 


(d) Interior regions not in a plane of symmetry. In 
this case the directions of the stresses at any point have 
first to be found by method 3. The slice used for this 
purpose will probably be found to have been cut in a 
plane which is not suitable for accurate measurement of 
both stress differences, and this will necessitate the use 
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of a second model. Another slice may be cut from this, 
and the stress differences found. The separation of 
the stresses in this case involves the integration of 
measurements taken in two perpendicular planes, and, 
in general, will involve the use of two more models'”, 

It may be noted that in cases (c) and (d) the most 
troublesome part of the solution is in the separation of 
the stresses. In problems where a knowledge of the 
principal shear stresses at a point is all that is required, 
these cases present little more difficulty than do the 
previous ones. 

There are available, therefore, the means of measur- 
ing the stresses at practically any point of a “frozen” 
model. The conditions necessary in order to apply these 
means to the determination of the stresses in an engineer- 
ing component are, in general, three : — 


(7) It must be possible to construct a scale model of 
the component in a homogeneous photoelastic material. 


(ii) It must be possible to reproduce the loading 
conditions of the prototype in an annealing oven and to 
keep these conditions constant during the “ freezing” 
period (16 hours or more). 


(iii) It must be possible to cut from the model slices 
sufficiently thin for the stresses in them not to vary 
appreciably through their thickness. 

These conditions being satisfied, the solution to the 
problem is possible and its difficulty depends only on 
the complexity of the model and on the amount of 
information which is required. 


Stresses in Fillets 


An important and frequently occurring problem in 
which the photoelastic method gives the required 
information very easily is that of the _ stress 
concentrations in fillets. The available data regarding 
stress concentration factors in fillets are all based on 
experiments or calculations in which the fillet is an 
isolated “ stress-raiser” in a test-bar in which the stress 
distribution in the absence of the fillet is known. In 
practice, the critical cases generally occur in regions neat 
to points of load application where no basic stress is 
known to which a concentration factor can be applied 


Ficure 10. Diagram of rail anchor. 
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An illustration of such a problem was given in a recent 
paper in The Aeronautical Quarterly” in which a 
photoelastic investigation revealed a very high stress 
concentration in the sharp corner of a spot-facing on an 
elevator hinge-bracket. In most such components the 
sharp corners are replaced by fillets, and this will 
undoubtedly lead to a reduction in the peak stresses, but 
the cxact amount of the reduction is still unknown and 
is likely to vary widely with the design and loading of 
the particular component. The choice of a radius for the 
fillet, therefore, is a matter of guess work backed by 
experience, and the designer is still in ignorance of the 
actual peak stresses involved. 

Figure 9 shows a comparison of the stresses in fillets 
of different radii in a loaded bracket. Two-dimensional 
models have been used in this illustration in order to 
obtain fringe-patterns from which the stresses may be 
compared by inspection, but equally good information 
may be obtained for three-dimensional models by taking 
accurate Measurements of the smaller optical effects 
which appear in thin slices cut from the models. 


Transference of Results from Model to 
Prototype 


One of the first questions raised by the engineer in 
relation to photoelasticity is that of the effect of the 
material of which the model is made upon the stress 
distribution. In this connection the common description 
of photoelastic materials as “ plastics ” tends to be rather 
misleading. In fact the stress-strain relation in these 
materials is strictly linear for a wide range of loads, 
extending in most materials right up to the point of 
fracture. and the fractures, when they occur, are brittle 
ones. This perfectly elastic behaviour has _ been 
confirmed in many two-dimensional cases by a com- 
parison of the stress distribution obtained from the 
photoelastic pattern with that given by the exact 
theoretical solution for a perfectly elastic material. The 
only respect in which two-dimensional photoelastic 
materials show any plastic behaviour is in their being 
subject to a certain amount of creep with time, both in 
strain and in the stress-optical effect. In the materials 
chosen for two-dimensional stress investigations this 
creep is small and the greater part of it takes place 
during a comparatively short period after applying the 
loads. In most materials, also, the magnitude of the 
creep is proportional to the stress and so there is no 
change produced in the relative optical effects at different 
points of the model. The only effect of the creep. 
therefore, is to produce a small change in the stress- 
optical coefficient of the material. and this may be 
allowed for by delaying measurements on both model 
and calibrating test-piece until the rate of creep has 
become negligibly small. The frozen-stress materials 
exhibit no creep at their critical “ freezing ” temperature. 
and the recovery creep after removing the loads is 
negligible. 

The stress distribution given by photoelastic methods 
therefore in both two- and three-dimensional investiga- 
tions, is that which would occur in the prototype if the 
material of which the latter is made retained a linear 
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FiGure 11. Fringe-pattern in curved part of anchor. 


stress-strain relation throughout the range of the stresses 
operating. If the peak stresses, as indicated by the 
photoelastic method, are higher than the yield stress of 
the prototye material, then the latter will undergo plastic 
yield in the most highly stressed regions with a conse- 
quent redistribution of stress which can only be 
estimated. The magnitude and extent of the high 
stresses in the elastic distribution, however, give some 
measure of the amount of plastic yield which will occur 
in the prototype. - 

The rail anchor illustrated in Fig. 10 provides a good 
example of such a case. The anchor, as indicated in the 
diagram, clips on to the base of the flat-bottomed rail 
where it is held by the spring grip of its jaws. The face 
of the lower part of the anchor bears against a sleeper, 
preventing longitudinal creep of the rail. It was found 
necessary to modify the shape of the anchor and a photo- 
elastic test was undertaken to determine the type of 
modification which would produce the best “ grip” on 
the rail. Fig. 11 is a photograph of the fringe-pattern in 
the curved portion of one of the models examined. On 
scaling up the stresses observed in the model it was 
found that the peak stress in the anchor, on the basis of 
an elastic distribution, would be very much greater than 
the yield-stress of the steel used. The distribution across 
the critical section XY of the anchor would therefore be 
something of the type indicated by the continuous graph 
in Fig. 12, in which only the stresses in the portion EF 
of the section would contribute to the grip. 

The crucial factor, therefore, was the amount of 
plastic yield which would occur and, since the criterion 
of yield might be taken to be the magnitude of the shear 
stress, it was not necessary to carry out any separation of 
the stresses, but only to compare the shear stress 
distributions in different models. Fig. 13 shows the 
graphs of the shear-stress across the critical section XY 
for three different models tested. 


EFFECT OF POISSON’S RATIO 


One other factor may affect the transference of 
results to the prototype and that is the value of Poisson’s 
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ratio for the materials. In two-dimensional problems 
the stress distribution is a function of this elastic constant 
only in the case in which an unbalanced force acts on 
the boundary of a closed hole in the model. Even in this 
case the effect is not large, for the value of Poisson's 
ratio for the materials used in two-dimensional work 
does not differ much from that for most structural 
materials. Results in typical cases for which the 
stresses can be calculated indicate that errors from this 
cause are not likely to exceed 7 per cent'”’. 

The value of Poisson’s ratio under frozen stress 
conditions, however, is very nearly 0.5, and in general 
this must have some effect upon the distribution of the 
stresses. Exactly how large this effect will be in any 
particular case it has been impossible, up to now, to 
estimate, for there are no theoretical solutions to 
compare with practical observations. The development 
of the scattered light method of observation does now. 
however, offer a means of assessing this effect. for it will 
now be possible to use this method on a directly loaded 
model made of the cold loading type of materia! and 
compare the results with those obtained from a similar 
frozen model. A cold loading material with the requisite 
light-scattering properties is now being produced, and 
experiments will be started shortly which, it is hoped. 
will result in a reduction in the possible errors due to 
this factor. 


[he accuracy of the Photoelastic Method 


The accuracy obtainable in a photoelastic investiga- 
tion necessarily varies enormously with the nature of the 
problem, and it is impossible to lay down any general 
limits. All that can be done is to examine the sources 
of error and to estimate in each individual problem the 
possible effect of each on the determination of the 
stresses in which we are interested. 

The accuracy of the actual optical observations is 
generally extremely high. Measurements by direct 
visual observation may be made to an accuracy of better 
than 0.01 of a “fringe” and the use of photometric 
methods reduces this possible error to something of the 
order of 0.002 fringe''’’. The peak stresses even in a 
thin slice from a frozen model. will generally produce an 
optical effect greater than half a fringe, and consequently 
errors from this source may be considered negligible. In 
the case of a thin slice, for example, the error in 
measurement of the thickness of the slice is likely to be 
greater than that in the optical observations. 

The main sources of error appear in the photoelastic 
material itself, in the construction of the model. and in 
the reproduction of the loading conditions which apply 
in the prototype. The material may be subject to lack 
of homogeneity and to initial stresses which are bound 
into it by chemical action during its manufacture. These 
defects can only be avoided by the exercise of great care 
in the casting'''’. The models used for three-dimensional 
work are either machined from a solid block of material. 
or cast in a mould. The heat generated during polymeri- 
sation of the resins employed makes the casting of large 
blocks which are homogeneous and free from initial 
stresses a matter of some difficulty. The casting of a 
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FIGURE 12. Comparison of stress distribution in the photo- 
elastic model of the rail anchor and in the anchor itself. 


model in or near to its finished shape, by reducing the 
mass of reacting material, helps to minimise the 
temperature effects. The technique of casting models in 
this way is still in the experimental stage. but enough 
progress has been made to show that in many cases it is 
quicker and more satisfactory than direct machining. 

Figure 14 shows a photograph of a model of a 
“tucked-in ” porcelain insulator cap, one quarter full 
size, which was cast directly in a mould. The investiga- 
tion was to determine the peak tensions in the cap which 
would be produced by the mechanical forces needed to 
hold it in position and produce an air-tight joint. The 
model was frozen under load, and measurements on a 
slice cut in a plane of symmetry yielded the necessary 
information. 

It may be mentioned here that one of the most urgent 
needs in photoelasticity today is for some systematic 
work by skilled “ plastics *” chemists on the development 
of photoelastic materials. The properties required for 
this application are different from those required for 
industrial uses, and although the firms manufacturing 
the resins have given much useful information and help, 


most of the research up to the present time has been ; 


done by physicists and engineers having very little know- 
ledge of the chemistry of the subject. 
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FiGurE 13. Graphs showing stresses across section XY of three 
different models of anchor. 
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n both directly loaded two-dimensional models and 
in siices Cut from a frozen model it is possible to apply 
corrections for initial stresses, but since such stresses 
ma) indicate a lack of homogeneity in both elastic and 
pho elastic properties of the material, these corrections 
cannot be relied upon unless the initial stresses are small. 
It is accordingly frequently necessary to discard models 
and blocks of material, and this still leads to a great 
wastage of time and material. If, however, care is taken 
in the selection of the materials, experiments have shown 
that excellent agreement is obtained between measure- 
ments taken on several similar models. 

The reproduction of loading conditions is, in many 
cases. the most serious source of error. The relatively 
large strains which occur in the model, especially in 
frozen stress conditions, may in some cases tend to 
produce serious differences in the distribution of applied 
loads. No general rules for avoiding such errors can be 
laid down. The conditions of each individual problem 
must be studied and some means devised of applying the 
loads in such a way as to minimise the differences in 
loading. In some cases, for example, very small loads 
must be used for investigation of stresses in the vicinity 
of load points and larger loads for the more remote 
regions. In others it may be possible, by studying the 
effect of varying the loading on the model, to obtain a 
good estimate of the stress distribution under the 
conditions existing in the prototype. 

As an estimate of the order of accuracy to be 
expected from photoelastic investigations, one might 
say that the peak stresses in a straight forward two- 
dimensional model would be liable to an error of from 
one to 5 per cent., rising to a possible 10 or 12 per cent. 
in cases Where unbalanced forces act on the boundaries 
of holes, and to an error of from 5 to 30 per cent. in a 
three-dimensional model. While the possible error in 
the last case may appear to be high, it must be borne in 
mind that these highest possible errors will occur chiefly 
in problems for which no other method offers a solution 
of any kind, and that in such cases even a qualitative 
comparison of the stresses in two different designs may 
be of considerable value. 


Future Developments 


There is no doubt that the methods and techniques 
outlined here will undergo great improvements as more 
workers enter the field and apply them to a wider variety 
of problems. In this country developments have so far 
taken place largely in research laboratories, and engineer- 
ing firms, on the whole, have been slow in taking up new 
methods, and not a little sceptical about such reports 
and resuits as have come to their notice. Thus, while 
it may fairly be said that, during the past few years. so 
far as fundamental research is concerned, the leeway lost 
to America in this subject during the war years has been 
made up, we still lag a long way behind in the matter 
of practical applications. 

There are signs, however, that the value of the 
photoclastic method is becoming more widely appre- 
ciated and the Stress-Analysis Group of the Institute of 
Physics, which has been largely responsible for the co- 
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Figure 14. Model of insulator cap cast in Marco Resin, show- 
ing part of slice cut in a plane of symmetry. 


ordination of work in this field, now has members 
representing a good proportion of the large engineering 
concerns in this country. 

As regards future development, there is one 
technique which is being investigated at the College of 
Aeronautics which offers a possible approach to the use 
of photoelastic methods in measuring the strains in 
actual components. In this method a layer of photo- 
elastic material is cast on to the surface of the component 
itself, and the optical effect in the layer when the 
component is loaded is observed by a reflection method. 
The photoelastic material thus provides a continuous 
strain-gauge over the surface. Satisfactory tests have 
been made on a number of test models, but so far as the 
author is aware, the accuracy and reliability of the 
method in general investigations have not yet been fully 
established. Difficulties of observations of the effect in 
some cases will probably limit the range of application 
of the method in actual investigations, but it offers the 
possibility of a means of measuring the amount and 
extent of plastic flow at the surface of a component, 
which might usefully supplement other methods of 
investigating the behaviour of the actual structural 
materials. 


Conclusion 

It is not suggested that photoelasticity will provide 
the solution to all stress distribution problems, or that 
a solution, when possible, will always be a quick and 
easy one. The aircraft designer especially will 
appreciate not only the need for using all available 
methods if he is to cover his wide variety of problems, 
but also the fact that the solution to many of these will 
be hard to achieve by any method. This paper has 
endeavoured to show, however, that there are certain 
vital problems in which photoelasticity will yield 
reliable information which no other method can give, 
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and other problems in which the photoelastic method 
will be quicker and easier than any other known at 
present. 
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Discussion 


R. L. LICKLEY (Fairey Aviation Co. Ltd., Fellow): 
Some of the photographs shown must have seemed like 
magic to people who had learned stressing before there 
was a photoelastic method. The Society was to be con- 
gratulated on having asked Colonel Jessop to lecture on 
that most important method, which was one which 
should be used to a greater degree in any work 
connected with stress and strain distribution. 

Some people doubted the effectiveness of the 
method, and he was very glad that a good portion of 
the paper had been devoted to answering the sort of 
questions which might be asked. The usefulness of the 
method depended on how the answers obtained from a 
piece of photelastic material could be related to the 
actual material which would be used in a structure: he 
felt that they could now say that those answers were 
fairly closely in agreement. In many cases to use the 
words “fairly closely” was to be conservative: it 
might almost be said the answers were exactly in 
agreement. 

Two matters in connection with the method which 
might be of great value to the aeronautical engineer 
were weight saving and fatigue loading. 

There were many occasions on which the engineer 
was unwilling to reduce the size or section of various 
parts of an aircraft because he was unable to calculate 
the stress to any degree of accuracy, and he could afford 
neither the time nor the expense of testing the parts: 
further, a part might not be availabie for testing for 
quite a long time. The use of one of the many types of 


photoelastic investigation which Colonel Jessop had 
demonstrated, however, would enable the engineer to 
find the places in a piece of material where the stresses 
were too low and also, where they were probably too 
high. The value of the method in helping to save 
weight was potentially great. 

There was no easy way to obtain a complete answer 
to fatigue loading and fatigue stressing, but they did 
know that in many cases fatigue failures arose from bad 
detail design, which permitted high stress concentration. 
Colonel Jessop’s illustrations, particularly those of fillets 


and holes, showed how the method indicated where the | 


high stresses occurred: and alternatives could be found 
quickly. On those two grounds alone the aeronautical 
engineer should use the method to a far greater degree 
than he did. 

The method had much value in teaching. When he 
was a Student he had found it difficult, and still did, to 
visualise strain and stress lines; it might not be easier 
to see them in colours, but it seemed easier. He thought 
that the universities might well devote more attention 
to the method and its use as a teaching medium. 

Could Colonel Jessop include, as an appendix to 
the paper, some information on materials ? One of the 
biggest problems in extending photoelastic work in 
the United Kingdom was probably the difficulty of the 
supply of useful materials which were reasonably 
inexpensive and at the same time easy to get. Two 
years ago, when he was last closely interested in that 
field, materials were neither easy to get nor were they 
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dge che. p. Therefore, if some data on that matter could Hickson’s photometric method of measuring light 
be i cluded in the paper it would be still more valuable. intensities. However for a straightforward application 
en Colonel Jessop had mentioned the work at the of results to design, the newer techniques had made it 
: College of Aeronautics with reflected light, which made possible to dispense with lenses by using a ground glass 
— possible the measurement of strains on the surfaces. screen for diffusing light. _There were several 
“a Some of the people connected with that work were advantages from the so-called diffuse light polariscope. 
- present at the meeting: he hoped that they would take such as no loss of time in setting up and aligning lenses, 
men- | part in the discussion, because that was another useful better rendering of fringes in direct vision, equal 
plied — approach to the problem. accuracy in counting the fringes, and compactness and 
cheapness. The principal objection to this, polariscope 
esses DR. R. B. HEYWOOD (Royal Aircraft Establishment): was caused by convergency of light, but this condition 
fae The potential advantage of photoelasticity had been did not arise when considering individual points as 
ie. demonstrated in that a wide range of stress problems distinct from areas of the model. Both the scope and 
: could be solved by employing modern techniques. In application of photoelasticity could be increased by this 
? he. effect the impressive claim had been made that all three simplification of technique and apparatus. 
June principal stresses and their three directions could be The photoelastic method, like all methods of stress 
found at each and every point in complicated three- analysis, could not be used to indicate the static 
craft dimensional models. That information was more than strength of structures. If structural loads were exces- 
1951, the designer was prepared to digest. so that in effect the sive, plastic flow would occur which would modify 
Cam- scope of photoelasticity had far outstripped the require- stress distributions. Thus the true elastic stress distri- 
ments of the designer. bution was not applicable. That was confirmed by 
hoto- | From the designer’s point of view, the impression Colonel Jessop’s implication in the results shown in 
netric | might be gained from a casual glance at the paper that Fig. 12 that maximum stresses of whatever value were 
+o } techniques were too complicated to be of practical use. reduced to a common yield stress. He thought the 
_ This was in contrast to the simplicity of the intrinsic author would agree that that was actually an over- 
wey problem which the designer wished to solve. Basically. simplification, since a part possessing a severe stress 
- he wished to know the stress at the point which was concentration was usually, although not always, weakest 
most likely to fail. Interest in stresses at other points in static strength. The photoelastic solution was at 
was only secondary. best of qualitative value only in its application to 
Practical experience showed that failure almost structures. . 
invariably started from the boundary of a part, rarely In certain other directions, such as in the determina- 
had from points within. The major problem was thus the tion of the first point and extent of yielding, and in the 
¥ determination of the stresses on external and internal confirmation of theories and mathematical solutions of 
sale: | boundaries. The solution was obtained directly from stress concentrations, photoelasticity could be of great 
wed fringe pattern, since the boundary stress was pro- value. 
Bis: ; portional to the fringe order. Thus no greater skill or The method could be applied with much greater pre- 
save | complexity was involved than merely counting fringe cision to fatigue problems in aircraft, for the fluctuating 
j orders. Even with three-dimensional frozen stress stresses due to gusts, and so on, were fully affected by 
nee models, the boundary stresses could still be found with stress concentrations when failure occurred in a large 
y did little further elaboration. Colonel Jessop’s method (1) number of cycles. That important application had been 
‘ bad _ in the Section * Three-Dimensional Techniques” was emphasised by Colonel Jessup. 
a : sufficient, there being no necessity to use scattered light One potentiality of the method which might have 
A ets Bas given in method (2). Admittedly there were a few received more attention in the paper, was that photo- 
e the ' cases in which stresses were required on a boundary elasticity should be used not only for the assessment of 
found that was in contact with an adjacent part, and then a stresses in a given design, but for the improvement of 
utical more advanced analysis became necessary. However, designs. A knowledge of stress distributions enabled 
egret Fit was significant of the type of problem usually certain new techniques and procedures to be employed. 
he f encountered that all except one of Colonel Jessop’s whereby perfection in shape could be achieved by 
as 2 examples had their strength determined uniquely by the ensuring that the stress distribution conformed to a 
id, to _ maximum stress at a free boundary and not at a loaded given idealised pattern. The example given in Fig. 9 
a 4 boundary, nor even at a point below the boundary. was not typical in that respect; the modest reduction in 
—_ ‘ | Thus the greatest contribution that photelasticity could stress from 6 to 5 units obtained by increasing the fillet 
canes ' give was by the use of simple techniques. radius could have been further reduced by other refine- 
— _ _ The simplicity of technique could be matched with ments in design. The photoelastic method could be 
Pe simple apparatus. In Coker and Filon’s classical used as a standard design procedure for evolving designs 
or ls | researches on photoelasticity it was necessary to use a of maximum strength, in addition to its more common 
rk in | complicated polariscope, employing several lenses in a application in the measurement of stresses. 
of the system lined up on an optical bench. Today the lens 
mnably type of polariscope was still the best type to use for CAPTAIN J. LAURENCE PRITCHARD (Honorary Fellow): 
Two F research investigations in which a complete and accurate This was an unsatisfying lecture to those to whom 
n that | stress analysis was required, particularly if the polari- photoelasticity was not well-known and to those to 
e they F scope were used in conjunction with Brown and whom it was. In the latter case it was preaching to the 
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converted, which was of little use to the former who 
would not know what the sermon was really about. 
The lecturer complained, “ In this country developments 
have so far taken place largely in research laboratories, 
and engineering firms, on the whole, have been slow 
in taking up new methods, and not a little sceptical 
about such reports and results as have come to their 
notice.” 

What could the lecturer himself expect when his own 
lecture was so definitely-indefinite ? 

“In comparison with other methods,” he said in his 
second paragraph, “ it is rapid and inexpensive, in many 
cases it gives much fuller and more exact information, 
and it is moreover the only method by which the mag- 
nitudes of the stresses in highly localised stress concen- 
trations can be determined with any degree of accuracy. 
This last point is of particular importance in modern 
aircraft design, for it is becoming increasingly urgent to 
have a more exact knowledge of the magnitude of these 
local stresses in order to avoid the possibility of their 
leading to fatigue failures.” 

That was definite enough, but when the accuracy of 
the photoelastic method was discussed later, the author 
hedged. “... accuracy... varies enormously with the 
nature of the problem . . . The main sources of error 
appear in the photoelastic material itself . . . No general 
rules for avoiding such error can be laid down . . .” 

The lecturer knew his job well, but sitting on the 
fence would not help those who read the lecture in 
Australia, Canada and the Outer Isles and would like 
to know more. It would be interesting to know. for 
example, what material or materials were used? Odd 
that a lecture dealing with the peculiarities of certain 
materials under stress did not give, specifically, even one 
material ! 

The main credit for the first practical developments 
for the determination of stress by optical methods must 
lie with E. G. Coker and L. N. C. Filon, following the 
earlier investigation of Clerk Maxwell. 

On 20th November 1918 Filon took the Chair at a 
lecture before the Society by Major A. R. Low on the 
work of the Chairman and Coker, who was in the 
audience and took part in the subsequent discussion. 
Coker himself, exactly a year previously, had lectured 
on Photoelasticity for Engineers at the Institution of 
the Automobile Engineers. 

At the end of his lecture. Major Low said his lecture 
should be regarded from two points of view. “In the 
first place as a method already, within limits, in full 
working order and immediately available as an instru- 
ment for obtaining numerical values in particular 
engineering applications. and as such of immediate 
commercial value. 

“In the second place, it is to be regarded as a 
method for whose further and fuller development 
research should be encouraged and subsidised.” 

Was it? One only knew that on the death of Filon 
and Coker the research work in the subject at University 
College appeared to die as well, and that between the 
wars the chief research papers were published in 
America, Germany and France. 


JOURNAL OF THE ROYAL AERONAUTICAL SOCIETY 


MARCH 153 


A Scotsman, Sir David Brewster, who did so much 
to found the British Association, was the first, in 1516, 
to draw attention that glass and similar isotropic trans- 
parent solids, exhibited the characteristic colours of 
crystalline plates under stress. Another Scotsman, 
Clerk Maxwell, made the first attempt to explore. by 
this method, the stress system in the plane of a strained 
plate: and Professors Carus Wilson, Coker and Filion 
carried out the further experimental investigation 
described in Major Low’s historic lecture before the 
Society just over thirty-four years ago. 

The liaison between research and practical achieve- 
ment in photoelastic methods was somewhat slower 
than it had been in atomic weapons. The gaps between 
research work at the Universities and the application of 
results in industry could be narrowed very appreciably 
by both sides doing a little research on each other's 
propaganda methods. The lecturer had the opportunity 
in the JOURNAL of the Society to put the subject across 
in a big way to aeronautical engineers. One hoped he 
would do what he undoubtedly could do in this direc- 
tion, with the full authority of his own work on the 
subject and the knowledge of being one of Filon’s 
pupils. 


V. M. HICKSON (Royal Aircraft Establishment): He 
did not agree that photoelasticity was the only method 
by which the magnitudes of the stress in highly localised 
stress concentrations could be determined with any 
degree of accuracy, especially in view of the subsequent 
statement that the peak stresses might be subject to an 
error of from 5 to 30 per cent. It was their experience 
that the claimed accuracy of brittle lacquer methods of 
15 per cent. held under most conditions, and they 


offered an alternative unless the surface radii were _ 


extremely small. Small radii were naturally avoided 
wherever possible in good design. There were many 
cases where the stress concentrations were inaccessible, 
or where contact stresses were involved. in which photo- 
elasticity was the only possible method. 

He also differed with the statement that the magni- 
tudes of the two stresses at the surface in the case of 
free boundaries in a plane of symmetry were best deter- 
mined by a combination of observations on a slice, and 
the scattered light method developed by Weller and 
considerably improved by Colonel Jessop. The oblique 
incidence method of Drucker. which he had described 
in two papers*, could be used with good accuracy in the 
case of unknown surface principal stress directions, and 
hence the method was even better adapted to the case 
quoted. It seemed to him that the combination of a 
slicing method with a scattering method was not so 
convenient as the slicing method alone, which offered a 
means of exploring the surface stress in the vicinity of a 
Stress concentration and thus covered cases (a), (b) and 


*Photoelastic Determination of Free Boundary Stresses on 
Frozen Stress Models by an Oblique Incidence Method. 
Brit, J. App. Phys., Vol. 2. pp. 261-269, September 1951, and 
Errors in Stress Determination at the Free Boundaries of 
Photoelastic Models, Brit. J. App. Phys., Vol. 3, pp. 176-181, 
June 1952. 
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(c). If the rate of change of stress in an inward direction 
froin the stress concentration were required, the method 
still applied. It did not apply to case (d), which was 
the determination of stresses of arbitrary direction in the 
interior, Since a knowledge of one principal stress direc- 
tion was necessary. He would suggest that in the 
majority of cases it was not necessary to know the 
interior stresses, since, in general, mean stress might be 
found from the applied loads and the conditions of 
equilibrium, 

li seemed to him that the proposal for using the 
scattered light method on cold models with a view to 
reducing the effect of ~ on the results had distinct pos- 
sibilities. He thought that the real virtue of scattered 
light was in this aspect. The non-correspondence of « in 
the model and prototype, coupled with large deforma- 
tions, Was a perennial problem to workers in the field. 
There were some theoretical solutions in which the 
effect. of ~ could be found—he was thinking in parti- 
cular of solutions in Timoshenko and of Neuber’s 
notches—and it appeared that the effect of « was 
greatest on the lesser of the two principal stresses at the 
surface. Those cases had been examined and he had 
three propositions for consideration : 


(i) The percentage effect of change of ~ was most 
significant in principal stresses on planes at 
right angles to the principal stress trajectories 
originating at the loaded boundary. 

(ii) The maximum effect at a point occurred in 
one or other of the lesser principal stresses. 

(ii) For practical purposes, the dependence of the 
stress on « at any point was linear. 


He emphasised that those propositions were based 
on known cases and could not have a theoretical proof. 

The conclusions were that the maximum stress was 
not affected to the greatest extent, and that extrapolation 
from results on two models of different ~ should be 
feasible. 

He agreed with Colonel Jessop’s remarks on the 
need for some research into photoelastic materials and 
would like to put forward a suggestion about frozen 
stress materials. The effective ~ for those materials was 
about 0:5, which might be attributed to the reduction of 
Young’s modulus as the temperature was raised, with- 
out a corresponding reduction in bulk modulus. 
Dr. Barwell* had given a theoretical expression for 
in a porous material, regarding it as a highly redundant 
Structure. In the limit the value was «/(2+©) i.e. 0-2 
if ~~ 0-5, which was verified approximately by tests on 
foamed materialst. If a porous plastic could be made 
with ratio of void-volume to total volume about a half, 
7 would be reduced to about 0:3. It would be necessary 
for the voids to be of sub-microscopic size to avoid 
scattering of light. Alternatively. if they were of micro- 
scopic size, the scattering would be useful in the 
Scattered light method. 

In the paper attention was drawn to the back 
reflection method applied to actual components. He 
*R. & M. 2467, 1951. 
+R. & M. 2686, 1952 (Pullen). 
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believed that the work of the College of Aeronautics 
was on plane specimens and the investigation was on 
yielding and plastic flow with consequent high values 
of strain. He was doubtful whether observations at 
boundaries on plane stress models were not subject to 
high error, since the necessarily thick photoelastic 
material was not in the same strain condition through 
the thickness. Moreover, if the method were applied 
to three-dimensional models, the main points of interest 
were stress concentrations with small concave radii, 
where conditions were extremely unfavourable for 
uniformity of stress through the thickness. In_ his 
opinion problems of that type were best covered by the 
brittle lacquer technique: its field of usefulness in the 
elastic range was exactly that covered by the surface 
layer technique, and its accuracy was probably higher. 

Colonel Jessop had done great service to the Aircraft 
Industry by drawing attention to the available photo- 
elastic methods which would be of use to the Industry, 
and also in endeavouring to assess the limitations. It 
was clear that the three methods of measuring strain and 
stress, strain gauges, photoelasticity and brittle lacquers, 
all had their applications to aircraft. and an organisation 
concerned with aircraft would not be said to be well 
equipped unless it could cope adequately with all three. 


F. A. KERRY (Saunders-Roe Ltd., Associate Fellow): 
In connection with the diagram showing the use of the 
photoelastic method for determining the distribution of 
stress around a series of holes in a cleat specimen, there 
would be a different distribution if the two parts of the 
joint were of different materials, say one part of light 
alloy and the other of steel. Was it possible to use 
photoelastic materials with different characteristics to 
represent this effect ? He supposed it might be done by 
making the thickness of the portion representing the 
steel about three times the correct thickness, but the 
effectiveness of this might be spoiled by the Poisson’s 
ratio effect. 


G. R. EDWARDS (the Chairman}: He was at a dis- 
advantage in relation to one or two comments which 
had been made about stressing aeroplanes in certain 
parts of the structures. All he could do was to support 
what Colonel Jessop had said, that what he did was 
quicker and cheaper and more accurate than the method 
used by Vickers-Armstrongs Ltd. 

On fatigue, the big bogey was the unloaded hole and 
other disturbances in the basic section of the spar. and 
he was quite confident that any device which could 
indicate in the very early design stages the sort of stress 
raisers which were being built into an aeroplane would 
make a big contribution, because they had not yet heard 
anything like the last about those problems, but only the 
beginnings of them. 


S. SCOTT HALL (Fellow): The pace of advance being 
what it was in aviation, there was a great risk of the 
Royal Aeronautical Society spending so much time con- 
sidering those advances and the applications of engin- 
eering to the art of flying that it did not give sufficient 
attention 


to fundamental advances in engineering 
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science itself and to the tools used to design and 
build aircraft; it was extremely important that now and 
again a lecturer such as Colonel Jessop should remind 
them of the advances in those fundamental techniques. 
The Society was indebted to Colonel Jessop, not only 
for his lecture, but also because he had instituted 
courses at University College on the subject. of which 
some present had taken advantage, and no doubt many 
more would do so in the future. 

Much had been said of Professor Filon and Professor 
Coker: in 1924 Professor Coker had lectured before 
the Institution of Aeronautical Engineers, one of the 
bodies now incorporated in the Royal Aeronautical 
Society. On that occasion he had said that the science 
and art of engineering design had been epitomised by 
Oliver Wendell Holmes, in “ The Deacon’s One Hoss 
Shay.” from which he had quoted, pointing out that the 
principles on which it was constructed were a counsel of 
perfection : — 

“*Fur, said the Deacon, *’tis mighty plain 

that the weakes’ place must stan’ the strain: 

Ah! the way t’ fix it uz I maintain is only jist 

To make that place uz strong uz the rest.’” 

That had been borne out by Colonel Jessop’s lecture 
very forcibly. 


D. A. THURGOOD (Vickers-Armstrongs Ltd.. Wey- 
bridge. Graduate), contributed: Colonel Jessop had 
rightly stressed the need for research into the develop- 
ment of new photoelastic materials, and it seemed that 
these developments must take place before the methods 
and techniques of testing and analysis underwent any 
great improvements. To the aircraft designer the three 
dimensional methods offered perhaps the brightest 
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prospects of enlarging his knowledge, or confirming his 
assumptions of the stress distributions within the 
structure. So far, however, nearly all three-dimensional 
work had been on what might roughly be termed 
* solid’ members, illustrated in the paper by the bracket 
in Fig. 9, and the rail member in Fig. 11. The three. 
dimensional model visualised by the aircraft designer 
was a reinforced cylindrical shell representative of wing 
or fuselage structure from which it was hoped to deter- 
mine spar boom and stringer end loads, spar web. rib 
and skin shear stresses in the case of the wing and frame 
bending stresses, stringer end loads and skin shear streses 
in the fuselage. In the construction of this type of model 
the casting technique could be employed, although not 
without difficulty, and there would appear to be con- 
siderable scope for the built-up type of model employing 
synthetic adhesives. Under test the loading could be 
arranged to give a reasonable number of fringes in the 
end load and bending members, but it was more than 
likely that less than one fringe would develop in the 
shear carrying members, which, for maintaining correct 
relative stiffnesses, were necessarily thin and the stress 
distribution could only be determined by laborious 
methods of compensation. It would appear that improve- 
ments in this direction could be affected by a range of 
photoelastic materials of similar mechanical properties, 
but with wide variation in the material fringe values. At 
present, two materials with different fringe values 
generally had different mechanical properties also. 

In his paper, Colonel Jessop mentioned that photo- 
elastic materials were not strictly “ plastic.” In view of 
the increasing application of plastics 
industry, might there not be a need for a “photo 
plastic material ? 


Colonel Jessop’s Replies 


In replying rather more fully to the speakers who 
contributed to the Discussion, he would like first to 
emphasise that the main purpose of the lecture was to 
bring to the notice of the Aircraft Industry the great 
advances which had been made in recent years in the 
photoelastic method of stress exploration, and to present 
some idea of the types of problem to which it could be 
applied and the sort of information a designer could 
expect from a photoelastic analysis. He had no inten- 
tion of “ preaching to the converted” (even though he 
might possibly feel that some of these were not quite 
confirmed in the true faith) and while he very much 
appreciated the presence of a number of practical photo- 
elasticians, his lecture was not addressed to them. 
Rather, he wished to convert such unbelievers as there 
might be among the senior men from the Aircraft 
Industry—the men in whose province it would be to 
authorise and support the establishment of a photo- 
elastic laboratory in their respective firms. (He would 
return to this subject in his reply to Captain Pritchard's 
comments.) Because of this he thought that no apology 
was needed for his omission to deal with most of the 
technical aspects of the subject. 


throughout | 


MR. LICKLEY: He was concentrating on and trying to ~ 


improve Marco resin. Of the other available materials — 


for three-dimensional work, the only satisfactory one in 
his experience so far was American Fosterite. It was 
very difficult to obtain that material in large pieces, how- 


ever, and it did not appear to be materially better except | 


that it had a slightly higher tensile limit than the Marco 
resin, which could be obtained in this country and which 
one could cast for oneself. He felt that progress in 
three-dimensional work, particularly in frozen stress 


work, was going to be rather on the lines of keeping to F 


low loads, small deformations and more _ accurate 


measurements of stresses, than on the lines which the | 


Americans tended to follow, i.e. of having bigger and 
bigger models, heavier and heavier loads and thick slices 
showing full fringe patterns. 
should stick to loads which conformed more nearly to 
the loads on the prototype. 
resin, but it still needed a lot of investigation. 


Written reply: The following was some general > 


information on materials. 


(a) For “ two-dimensional” investigations 
The three most useful materials they had found 


He thought that they 7 


He recommended Marco 
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general | 


found 


were Columbia Resin (C.R.39), Xylonite and Bakelite 
(B.i 61/893). For most general prrposes C.R.39 was 
satisiactory. It was obtainable in polished sheets up to 
2 ft I ft.. in thicknesses from in. to 4 in. 
Selected sheets were sensibly stress-free. The cost was 
roughly proportional to thickness. Sheets in. thick 
cost £3 per sq. ft. (Obtainable from Ashdowns Ltd.. 
St. Helens, Lancs.) 

Xyvlonite was much less optically sensitive than 
C.R.39, but was less brittle, was obtainable in very large 
sheets in a wide range of thicknesses, and was much 
cheaper. Sheets had to be carefully selected for 
freedom from stress. Built-up models were easily made 
by cementing with acetone. The cost varied with 
thickness. (Obtainable from B.X. Plastics Ltd., Larks 
Hall Road, London, E.4.) 

Bakelite had slightly greater optical sensitivity than 
C.R.39, and about three times its tensile strength. 
Obtainable in sheets I] in. x 7 in. x } in. The sheets 
needed machining and polishing, and also a rather long 
annealing process to remove initial stresses, and there 
was a considerable amount of waste in stressed edges. 
The cost was very high—about £10 10s. Od. per Ib., and 
an import licence was needed for purchase from 
America. (Now marketed by Catalin Ltd., Waltham 
Abbey. Essex.) 


(b) For “ three-dimensional” work (frozen-stress) 


Here the three most reliable materials they had tried 
were Bakelite, Fosterite and Marco Resin. The last 
was the one on which all their work was being done at 
present. 

The Bakelite was the same material as for two- 
dimensional work, but the sheets could be obtained in 
thicknesses up to 1} in. for use in frozen-stress work. 

Fosterite was another American material, obtainable 
from the Westinghouse Corporation, Regent Street, 
§.W.1. It could be supplied in cylinders up to 8 in. 
diameter and 2 or 3 ft. long. The rods contained initial 
casting stresses and required annealing. Annealing and 
“freezing” must be done in an inert medium (e.g. 
silicone oil). Cost high, approximately £7 10s. Od. per 
lb.. and an import licence was needed. 

Marco Resin (Sb.26C) was very similar in properties 
to Fosterite. Optical sensitivity and tensile strength 
both somewhat lower. The materials for casting the 
resin were supplied by Scott Bader & Co., Kingsway, 
W.C.1. Photoelastic models might be cast in their final 
shape in suitable moulds. Casting techniques. still 
needed improving, and up to the present difficulty was 
found in casting models whose smallest thickness was 
much above 2 in., but many models had been cast and 
investigated by “different workers, with completely 
satisfactory results. Cost worked out at about 
8s. per Ib. 

Fuller details about the properties of most of these 
materials could be found in any text book on photo- 
elasticity. Further information about Marco Resin 
could only be obtained at present from those photo- 
elasticians who had been experimenting with it. There 
Were many other materials which exhibited the photo- 
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elastic effect, and fresh ones were continually being 
tested. There was a reasonable hope that even better 
materials would be found than those at present in use. 


DR. HEYWOOD: It seemed to him that in cases where 
fillets occurred very close to load points in components 
which had irregular shape, they could not lay down any 
general law governing how the radius of the fillet would 
affect the stress; the kinds of stress and the applications 
of load might be so various that the change in fillet 
radius might have quite different results between one 
component and another. He advanced that view as 
purely his own; he was not an engineer, but a physicist, 
and consequently in the paper he had tried to keep clear 
of problems of fatigue, and so forth. He had dealt with 
the methods and had left it to the engineers to see how 
they could be applied to engineering problems: probably 
he would be wrong if he tried to indicate how they 
should be applied. 

Written reply: He was inclined to agree, on further 
consideration, that his paper might tend to give the 
impression that the photoelastic method was a more 
complicated business than it was in practice, in a very 
large number of cases. He regretted that he had not 
emphasised more how iarge was the proportion of very 
simple cases. The more complicated cases, however, 
did occur, and he thought that an aircraft designer 
would be pursuing a very short-sighted policy if he were 
to ignore the possibility of these and be content to cater 
only for the simpler type of investigation. The designer 
knew much better than he what his own problems were. 
and his object was to show that methods were now 
available for the solution of many of the more difficult 
problems which, a few years ago, were quite beyond the 
scope of any experimental method. 

He felt that Dr. Heywood was over-simplifying the 
problem, however, when he suggested that only the 
maximum stresses on a free boundary were of interest 
to the designer. In the case of fatigue stresses at a 
point of high stress-concentration on a boundary, for 
example, the extent of the plastic yield must depend 
upon the rate at which the stress decreased in all direc- 
tions around the point. '* was this factor which was 
largely responsible for the not uncommon failure of a 
component at some point other than the point of highest 
surface stress. 

Another example, which was illustrated in the 
lecture, showed the peak shear-stresses in a bar to occur 
on the boundary of a hole at a point in contact with a 
bolt. In that case the maximum tensile stress would 
probably also occur somewhere on the same loaded 
boundary, but its location and value were not given 
directly by the fringe-pattern, but needed the application 
of one of the integration methods. That type of prob- 
lem would occur in “three-dimensional” as well as 
“two-dimensional” cases, and surely it was of 
sufficiently frequent occurrence, and of sufficient 
importance to justify (if justification were needed) a 
mention of the methods by which such contact stresses 
could be evaluated. 

He thought that a discussion on the relative merits 
of various types of polariscope would be too technical 
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for this journal. He agreed, however, with Dr. Heywood 
on the advantage to be gained in many cases by direct 
observations in a diffused field of light: that was the 
original method of observation which Filon and he used 
some thirty years ago. It was still in use for suitable 
problems at University College Laboratory. The 
“projection” bench was converted to a “diffused 
screen” bench, simply by placing a sheet of tracing 
paper between condensing lens and polariser. Obser- 
vations were made with a low power travelling micro- 
scope. 


CAPTAIN PRITCHARD: He had worked with Filon in 
1919 and for several years afterwards on photoelasticity. 
and no one could have a greater admiration than he had 
for Filon’s work. He would have liked to have written 
more about Filon in the paper. As regards work in the 
U.S.A.. it was of interest to note that Professor Frocht. 
who had written the classic American work on photo- 
elasticity, had dedicated his book to Coker and Filon 
and had paid a very high tribute to their work. 


MR. HICKSON: He did not emphasise, when he had 
given the limits of 5-30 per cent. for a three-dimensional 
model, that that was the range from the simple three- 
dimensional problem of the surface stresses to the very 
complicated case of the internal stresses in a non- 
symmetrical model. The 30 per cent. applied to the 
latter, and not to a simple surface stress such as could 
be explored by brittle lacquer: for the simple surface 
stress the 5 per cent. was nearer the mark. He did not 
think the consistency of brittle lacquer results had been 
proved, and he was a little afraid of it. 

He was interested in the tentative propositions on the 
effect of Poisson’s ratio, and would like to look into 
that matter a little more. The biggest effect of Poisson’s 
ratio might be on the smaller stresses: but it would affect 
also the maximum shear stress, which might be in many 
cases the relevant factor. 

With regard to the surface layer method which was 
being developed at the College of Aeronautics, he quite 
agreed that on curved surfaces where stress concentra- 
tions occurred they could probably not use a surface 
layer of photoelastic material sufficiently thin to give 
anything like a reproduction of the surface strain. But 
he did think that on flat bars it might give useful 
information on the occurrence of plastic yield in the 
materials, and thus could be used for testing materials, 
though not for exploring stresses. 

Written reply: He agreed that he was in error in 
stating, without modification, that a combination of 
direct observation through a slice and observation by 
scattered light would give the required result “ most 
easily.” He should have said that that was one method, 
and the one which he had found most satisfactory in 
some stress concentration problems. He had no inten- 
tion of implying that any one method or combination 
of methods would be the best in all cases. In fact he 
was sure that different individual problems needed 
different methods, and that frequently one could not 
know which method of measurement would be most 
suitable until the frozen model had been examined in a 
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polariscope. He considered, therefore, that it would be 
unwise for an investigator to limit himself to the use of 
any one routine. 

With regard to Mr. Hickson’s tentative propositions 
on the effect of Poisson’s ratio, he felt that there was iiot 
yet sufficient evidence to justify the laying down of uny 
rules at all. Those propositions appeared to apply to 
the two cases referred to: Timoshenko’s solution for a 
spherical cavity in an infinite block under tension, ind 
Neuber’s solution for an infinite cylinder with hyper- 
bolic groove. 

Another exact solution had recently been produced 
by R. Capildeo and A. C. Stevenson, who had given him 
permission to quote results from it before its publica- 
tion. This was the solution for a sphere under 
diametral compression along a diameter, and the results 
for the stresses at the equator and at the centre (the only 
ones yet evaluated) did not confirm Mr. Hickson’s 
propositions. 

(a) In both cases the stress in the direction of the 
load suffered by far the greater change with a change 
in 

(b) In the case of the stresses at the centre it was 
the greater stress which suffered the greater change. 

(c) The departure from linearity in ™ was more 
marked than in either of the previously considered 
solutions. 

(d) The effect of a change in ~ was much greater 
than in the other solutions. The change in the stress- 
difference on the equator due to a change in « from 0-25 
to 0:5 was 30 per cent. of the greater stress, while that at 
the centre was 56 per cent. of the greater stress. 


It might be noted that in all three solutions «* appeared 
in each stress in the same form (a+ br)/(c+dr). The 
departure from linearity would not cause serious error 
in any of these solutions in the case of interpolation 
between two values of «, but errors of up to about 30 
per cent. might be introduced if one extrapolated, say to 
a value of ~=0-25 from results obtained with =0-35 
and «=0°5. 

While the sphere did not present the case of a stress- 
concentration at the surface, it was the only one of the 
three solutions which had completely finite boundaries, 
and it would. he felt, be rash to dismiss the evidence it 
provided as inapplicable to stress concentration prob- 
lems on a boundary. 

The suggestion of making a porous plastic which 
would have a Poisson’s ratio of about 0-3 was an 
interesting one, and he hoped Mr. Hickson would 
discuss it with the chemists and explore the technical 
possibilities of suitably “foaming” a resin during cast- 
ing. Like all the problems on producing and casting 
suitable resins it would probably be tackled much more 
efficiently by skilled industrial chemists than by 
engineers and physicists. 


MR. KERRY: Something of the sort suggested by Mr. 
Kerry could be done and it was the sort of thing that 
would have to be done by experimenting on models in 
two or three sections. In one case they could use a bar 
made of the photoelastic material and the cleat made of 
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a material which bore to the photoelastic material a 
suitanle strength ratio; and they could do another 
expe:iment with the photoelastic cleat and a bar of the 
other material, to study the stresses there. In his view, 
that method of approach would be preferable to using 
different thicknesses of materials. He did not know 
how the stress distribution might be affected by changes 
in the relative thicknesses. 

In dealing with that problem, no rivets or bolts were 
used, only pins. If rivets or bolts were put in, the prob- 
lem would immediately become three-dimensional. 
which was much more difficult to examine. He assumed 
that any pressure applied by bolts or rivets would first 
tend to induce frictional forces and so relieve the load 
on the rivets themselves; but it would also induce lateral 
compression, which would increase the total shear 
stresses around the hole, and that increase of pressure 
might be liable to cause more trouble than the absence 
of pressure on the bolts or rivets. 


MR. THURGOOD: He did not foresee the possibility 
of using photoelastic methods for the evaluation of 
stresses in models of complete wing or fuselage struc- 
tures. One of the main difficulties, he thought. would 


occur in the buckling of the thin skins under shear stress 
long before the applied loads were great enough to 
_ produce measurable optical effects in any of the mem- 
bers. In such complex structures maintenance of the 


exact scale ratio would probably be the only way of 
ensuring “similarity” between model and prototype. 
and that would mean that most of the members would 
be exceedingly thin. He thought that such cases were 
examples of the type of problem for which the electric 


_ strain-gauge method was most suitable. 


The “ plastics” to which he thought Mr. Thurgood 


_ referred were no more plastic in behaviour than were 
_ the photoelastic plastics. They consisted, in fact, of 


some kind of “filling” bound by synthetic resins, and 
most of the photoelastic resins were originally produced 
for just such industrial applications. He had not seen 
any data regarding the stress-strain relations of the 
various “ plastics,” but he thought they were likely to 
have a sensibly linear characteristic. It would be 


_ exceedingly useful if a range of photoelastic materials 


could be found which reproduced the stress-strain 


_ behaviour of the various structural materials, but he was 


afraid that was more than even the most optimistic 
photoelastician could hope for. 


CAPTAIN PRITCHARD, written reply: He had left his 


teply to Captain Pritchard until the last because 


Captain Pritchard had raised the question which he felt 
was of the greatest importance to the development of 


» photoelastic applitations in this country. He was not 
Sure whether he should take Captain Pritchard's 
| Criticism of the lecture as “ unsatisfying” as a compli- 


ment or the contrary. If he meant that he and his fellow 


_ aeronautical experts were left wanting to know more 


about the subject, then the lecture fulfilled one, at least. 
of its objects. He hoped, if that were the case, that the 
desire for further information would not evaporate 
without leading to definite action. He could assure any 
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would-be enquirers that the photoelastic laboratories at 
the Universities, as well as those at the R.A.E. and the 
N.P.L., were always ready to give any assistance they 
could. If. on the other hand, Captain Pritchard’s 
criticism was on the “ indefiniteness ” of the information 
given and the alleged “hedging” on the matter of 
accuracy, he could only say that he was endeavouring 
to present a true picture of the position, which could 
not be done by quoting a single figure for the probable 
error. An estimate of the order of accuracy to be 
expected in various types of problem was given later, 
and this had been further explained in his verbal reply 
to Mr. Hickson. He felt that no engineer could really 
expect more, especially in view of the fact that, problem 
for problem, the accuracy of the photoelastic method in 
the cases for which it was recommended was generally 
higher than that of calculations or of other experimental 
methods. 

Captain Pritchard deplored the slowness with which 
photoelastic methods had been taken up in industry, and 
suggested that University research departments and 
Industry were equally to blame. He disagreed with him. 
The responsibility for lack of active co-operation rested 
squarely on the shoulders of Industry. During the years 
1912 to 1930, Coker lectured on photoelasticity to all 
the Engineering Institutions and to many large 
industrial concerns. In 1918, as Captain Pritchard 
recalled, Major Low lectured to the Royal Aeronautical 
Society and called attention to the need to subsidise and 
encourage research in the subject. What happened ? 
The engineers, almost without exception, expressed their 
belief in the potentialities of the new method—and sent 
the few pioneers back to their laboratories (probably 
wishing them “ the best of luck ”) and saying they would 
be interested to hear of any new developments! If. 
during these vears, Industry had shown any initiative, if 
from their large resources, they had done what Industry 
did in America after Coker’s lecture tour there— 
established research units to hasten and further the 
development of the work—then in all probability those 
discoveries which eventually took place in America in 
the late 1930’s would have taken place in this country. 
and moreover research here would not have ceased tem- 
porarily on the deaths of Filon and Coker. 

It was a disaster that Filon’s death and Coker’s final 
illness occurred almost simultaneously just before the 
war. Many of them were away during the war, serving 
in some technical capacity or other in the Forces, and 
the staffs who remained with the Universities had no 
time for research. Industry, however, was more vitally 
interested than ever in stress problems, and was in a 
position to make use of and develop the new ideas from 
America. Some few large firms did, in fact, keep in 
touch with new developments and did photoelastic work 
on their own particular problems, but in general there 
was no sign that Industry was interested in the subject. 
Finally, it was left to the Universities after the war to 
effect co-operation between workers in this field by the 
formation of the Stress Analysis Group which was later 
affiliated to the Institute of Physics, and, be it noted, 
among all the Institutions and industrial concerns 
invited to co-operate, the Institute of Physics was the 
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only one which made any offer of financial or other 
assistance ! 

(It was not inappropriate to mention here that this 
Group, with a membership from academic, government 
and industrial research departments, was concerned 
with all forms of experimental stress-analysis, and its 
function was mainly the exchange of information and 
ideas on new methods and techniques. It was probably 
the best means by which stress-analysts could keep in 
touch with current work all over the country, and the 
Group would welcome still more members from 
Industry.) 

The position today was in many respects similar to 
that in the 1920's, for although there were now many 
more workers in this field there were still far too few 
and there was much more to be done. Two results 
followed from that—industry as a whole, and perhaps 
the Aircraft Industry especially, was not getting the 
advantages which photoelasticity offered; and photo- 
elastic methods were not receiving the impetus to their 
development which only work on a wide variety of 
problems could give. 

Captain Pritchard suggested that this had been an 
opportunity to “ put the subject across in a big way to 
aeronautical engineers.” He felt that the photoelastic 
method itself was already accepted in principle by 
practically every branch of the Aircraft Industry, and 
needed no “ putting across.” He would, however, act 
on the spirit of his suggestion by endeavouring to 
emphasise what were the relative roles of the academic 
and the industrial research establishments in the 
development of the subject. 


THE UNIVERSITIES AND TECHNICAL 
TWO MAIN FUNCTIONS : — 

(1) THE DEVELOPMENT OF THEORY AND NEW 
METHODS AND THE DEVISING OF SUITABLE TECHNIQUES. 

(2) THE TEACHING OF THE SUBJECT AT UNDER- 
GRADUATE AND POST-GRADUATE LEVELS. 

With those two functions the whole of the resources 
of those institutions in staff, technicians and workshop 
and laboratory facilities were fully occupied. During 
the course of their work it was possible, and even neces- 
sary. from time to time, for the departments to work 
on some particular engineering problem, but the investi- 
gation of such problems was not their job. On more 
than one occasion at the photoelastic laboratory at 
University College they had been asked if they could 
undertake comprehensive investigations to obtain data 
which could be made available to industry and which 
might, in some case, form a basis for the consideration 
of one of the British Standards Committees. Neither 
staff nor students would be justified in wasting their 
time on investigations of this kind which, of necessity, 
were largely routine work. Was not this a matter which 
should be dealt with by the establishment of a testing 
laboratory by some organisation such as the D.S./.R.? 
They already had the nucleus of this in the well 


COLLEGES HAD 
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equipped photoelastic research laboratory at the N-P.L, 
This might easily be expanded by the addition of tore 
apparatus, and there would be little difficulty in training 
a staff of laboratory assistants who could carry oui the 
work under direction. 

From the standpoint of teaching, also, the position 
was probably not generally understood by Industry. In 
the first place, a course of lectures on photoelasticity, 
or the reading of a book on the subject, was of little use 
unless accompanied by practical work by the student. 
The provision of enough photoelastic equipment to give 
students the bare minimum of practical work was 


probably the most serious bottle-neck in the training, | 
Unless a laboratory were to confine its post-graduate | 


teaching to one or two students a year, it must be 


content to give short courses, which, to be effective, must } 


be followed closely by much more extended practical 
work. That was where action was needed on the part 
of Industry. 


elastician, unemployed” until such time as it was 
decided to install a photoelastic laboratory. That, in 
fact, was what had happened to too many of the men 
who had attended the courses at University College. 
Their object in these short courses was to give sufficient 


training to enable a man to undertake the establishment | 
of a laboratory and the carrying out of stress- | 


explorations. In the early stages he might have to seek 
advice from time to time but he could still be doing 
useful work, and it was of primary importance that he 
should be given the facilities to acquire experience. 
THUS THE FUNCTION OF INDUSTRY WAS TO PROVIDE 
THE FACILITIES FOR THE CONTINUED TRAINING OF THEIR 


PHOTOELASTICIANS AND FOR THE SOLUTION OF THEIR OWN 


PROBLEMS. 

In the Aircraft Industry a few firms already had a 
photoelastic laboratory in use. Some were in process 
of installing one, others, he understood, were consider- 


It was not enough to send a man toa} 
vacation course, and then to retain him as a “ photo- | 


ing the matter. That was progress in the right direction, — 


but it was not enough. 


certain number of men could gain experience ? 


To sum up: The scientists had worked out the” 
theories, evolved the methods and devised techniques. — 
All those had been tested, and found capable of being 
of immediate value to industry. For the improvement ~ 


and development of methods and techniques, and for 
the extension of the field of application, what was 
needed was a large increase in the number of men 


working on industrial problems. It was up to Industry | 


to move in the matter and, he thought, the move must 


come from the heads of Industry; from the men who} 
were in a position not to be afraid of taking up some | 
thing new. Progress would not be made without theit © 


authority and support. 


Surely such a large and 
important industry could do something in the way of | 
co-ordinating effort. Could not the S.B.A.C. establish a | 
central laboratory where problems common to the| 
whole industry could be investigated and where a) 
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The Resistance to Air Flow of Wire Gauzes 


by 


W. J. D. ANNAND, B.Sc. A.F.R.Ae.S. 
(Rolls-Royce Ltd.) 


SUMMARY :—Wire gauzes are used in air duct systems and wind tunnels for many 


purposes. 


experimenters are correlated. 


Methods of estimating the pressure loss through a gauze currently available 
are often found in practice to be unsatisfactory. 
It is shown that the effect of Reynolds number on 
pressure drop is similar for a large range of gauze porosities. 


In this paper data obtained by several 


An empirical method 


of estimation of the pressure drop is given which can be used with confidence for gauzes 
of porosity (defined as the proportion of area not blocked by the wires) from 0-27 to 1-00, 
and over a range of Reynolds number, based on approach velocity 
and wire diameter, from 20 to 600. 


|. Introduction 

Wire gauzes are used in air duct systems and 
wind tunnels for a variety of important duties. Some 
examples are: — 

(i) to catch stones or other injurious material, as 
in engine air intakes. 
(ij) to produce uniformity of the velocity distribu- 
tion across a duct, as in many wind tunnels 
and other test apparatus. Non-uniformities 
are reduced because the gauze presents greater 
resistance to flow where the velocity is higher. 
Refs. 2 and 7 present some theoretical con- 
siderations and experiments on this application. 
The use of gauzes in diffusing ducts presents 
special features in that the effect of the gauze 
penetrates upstream, reducing the causes of 
non-uniformity as well as the non-uniformity 
itself. Extensive experiments are reported in 
Ref. 3. 


to reduce the turbulence of a stream. Much 
use is made of fine gauzes for this purpose in 
low-turbulence wind tunnels. The components 
of turbulence are modified in the contraction 
on entering the gauze and, provided that the 
gauze wires are fine and the porosity suitably 
chosen, the overall turbulence some distance 
downstream can be considerably reduced below 
that in the approaching stream. Experiments 
on this effect are reported (for example) in 
Ref. 2. and Ref. 7 gives some theoretical con- 
siderations. 


(iii) 


(i) oppositely, to produce artificially high turbu- 
lence. For this purpose, coarse grids of rela- 
tively large diameter rods are used; the large 
eddies in the rod wakes presumably produce 


by their breakdown the turbulence required. 


(:) to introduce known pressure drops into experi- 
mental systems. 
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For all these applications it is desirable to be able to 
estimate accurately the pressure drop which will be pro- 
duced by a given gauze—particularly so when the gauze 
is to be used to reduce the non-uniformity of a stream, 
since in this case’) a certain definite value of the ratio 
of pressure drop to approach dynamic head (4pV*) is 
necessary for the most complete removal of the irregu- 
larities. Nevertheless experience shows that the empiri- 
cal formulae usually employed are often considerably 
in error. Generally speaking, these have been based on 
very limited data, and it is the purpose of the present 
paper to review and correlate the available information 
and to derive an empirical relation which can be used 
with confidence throughout a wide range of conditions. 


Notation 


a coefficient in the formula for the estimation 
of k (see Section 5) 
mean wire diameter 


D,,D., “depths of weave” in the two directions (see 
Fig. 1) 
1,,1, pitches in the two directions of the weave 


(see Fig. 1) 
k_ pressure drop coefficient = pressure 
drop/4pV* 
Reynolds number=Vd/» 
porosity 
approach velocity 
kinematic viscosity 


< 


2. Definitions 


The main dimensions describing the geometry of a 
gauze are the mean wire diameter d, and the pitches /, 
and /, in the two directions of the weave (see Fig. 1). 
From these can be calculated the porosity 8, defined as 
the area of the holes in a unit area of gauze, as projected 
orthogonally onto a parallel plane. They are related by 
8=[1—(d/1,)) (1 —-(d/1,)). 
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Two other dimensions characterising the gauzes are 
the *“ depths of weave” D, and D, in the two directions. 
Fig. | makes clear the meaning of this term. Only one 
investigator has measured these dimensions for the 
gauzes tested. Although they might be expected to have 
some influence on the flow pattern through the gauze, 
the indications are that the effect on pressure loss is very 
small for gauzes of normal manufacture. 

Any one of the gauze dimensions can be used to 
form a Reynolds number characterising the flow. The 
diameter d has been chosen, giving a Reynolds number 
R=Vd/», where V is the approach velocity and v the 
kinematic viscosity. It will be shown later that this 
choice gives good correlation. 

The pressure drop is defined as the loss in stagnation 
pressure from a plane immediately upstream of the 
gauze to one sufficiently far downstream for the wakes 
of the wires to have become indistinguishable by mix- 
ing—that is, about 30 wire diameters downstream, to 
judge from certain measurements on a coarse grid re- 
ported in Ref. 4. In practice the pressure drop is usually 
measured by the difference in static pressure between 
planes some distance upstream and downstream of a 
gauze sample in a duct of constant section. 

The pressure drop coefficient A is defined as the ratio 
of this loss to the approach dynamic head $pV°. 


3. Sources of Data 


Test results have been drawn from four sources, in 
all of which each gauze was tested at a number 
of air speeds. 

(a} Eckert and Pfluger’’ tested 16 commercial 
gauzes, measuring the loss of static pressure between 
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dad i] 
FiGure |. Diagram showing characteristic dimensions 


of a gauze. 


tappings about 5 cm. upstream and downstream of the 
gauze in a tube of 45 cm. diameter. Their results are 
presented only in graphical form, and it is difficult to 
pick off the individual test measurements accurately, 
For this reason only the results for four gauzes have 
been made use of for analysis (unfortunately, not includ- 


ing the gauze of greatest blockage: the dimensions | 
quoted for this are inconsistent). 
(b) Taylor and Davies"? tested four gauzes, measur- } 


ing the loss of static pressure between tappings about 
2 in. upstream and downstream of the gauze in a tube 
of about 3:2 in. diameter. Their results are tabulated 
and the measuring accuracy is good. 


(c) Simmons and Cowdrey tested eight gauzes, 
measuring the loss of static pressure between tappings 
6 in. upstream and downstream of the gauze in a 6 in. 
square duct. Their results are tabulated and _ the 
measuring accuracy is very good. 


(d) Schubauer, Spangenberg and Klebanoff tested 
lJ gauzes, using a rather different arrangement from — 
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anv o: the foregoing. The gauzes were fastened across 
the e: | of a square duct and air was blown through 
them :2to the open atmosphere. The pressure drop was 
measured as the difference in pressure measured by 
static probes placed centrally upstream and downstream 
of the gauze. Here again the results are given only 
graphically on a rather small scale, and the accuracy 


with which points could be picked off is rather low. 


Results for six of the gauzes were used in the analysis. 

The measurements include the pressure drop due to 
wall friction over different lengths of ducting. It is 
impossible to estimate this loss accurately because the 
presence of the gauze completely alters the velocity 
profile near the wall, both upstream and downstream. 
However. it is improbable that the correction to the 
measured pressure drop coefficients would exceed 
0:03, so that the omission of the correction is not of 
importance. 

The dimensions of the gauzes included in the analysis 
are given in Table I. Most are of square mesh or nearly 
so; the greatest deviation from squareness is /,=0-80/,. 


4, Variation of Pressure Drop Coefficient 
with Reynolds Number 


The first step in the analysis was to plot all the 
pressure drop coefficients logarithmically against Rey- 
It was soon 
neticed that the forms of the curves obtained for the 
various gauzes were very similar, and on investigation 
it was found that identical curves, shifted only along the 
log k axis, could be fitted to all; or alternatively, by 
shifting each set of points parallel to the log k axis, all 


the points could be brought onto a single curve. This 


is illustrated in Fig. 2 where (log k— log kx-,,,,) is plotted 


TABLE I 
DETAILS OF GAUZES TESTED 


Weave depths Poro- 


| D,/d B 
| 0-277 


Wire Pitches 
Source diam. 
d:in, cane 


Ref. 5 00059 0-0125 


00055 0-0154 0-413 
00098 00547 0°672 
O-O118 00856 
Ref.6  0:-0087 | 0:0251 | 00251 | 1-00 | 1:25 | 2-61 | 0-429 
0-013 | 0:0417 | 0:0417 | 1:00 | 1-53 | 2:23 | 0:475 
0-024 0-126 | 2-05 0-649 
Ref. 4 00044 0:0099 0°96 0-318 
0-0053 | 0:0126 | 0:0125 | 0:99 0:333 
0-0068 0:0167  0:0167 | 1:00 0:351 
00067 0:0166 | 0°86 0:389 
0-0107  0:0333 | 0-95 0-451 
0:0250  0-0249 | 1-00 
0:0417 | 0:0333 | 0°80 0:649 
00147 0:0869 0:0832 | 0:96 0-682 
1 Ref.2 | 0:0135 | 9 0-025 0-212 
0-0170 0-050 
0-0055 0-0185 0494 
0-007 0-025 0-519 
0-0075 0-042 0:672 
0-025 


0:250 


0-810 
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DRAG COEFFICIENT 


§ 
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20 3:0 40 60 


LOG R 
FiGure 3. Drag coefficient of a circular cylinder. 


against log R for two gauzes from each source. In the 
complete plot, individual points would be indistinguish- 
able near the mean line; of 178 points representing tests 
on 22 gauzes, only 16 fall away from the mean line by 
more than 3 per cent. of k, and of these 16 all but five 
are points picked off from the graphs of Refs. 5 and 2. 

This degree of correlation is quite striking, and in 
some respects was unexpected. Some variation of the 
shape of the Reynolds number curve with ~ was 
expected, because the variation of k with R probably 
arises partly from variation of the discharge coefficient 
of the orifices formed by the meshes, and this variation 
would be expected to depend on the contraction ratio of 
these orifices—that is, on 1/8. Further, it would be 
expected that at very high ° the variation of the loss 
coefficient with R would approach that of the drag 
coefficient of an infinite circular cylinder; this curve is 
plotted in Fig. 2 for comparison and it will be seen to 
differ somewhat from the curve for k. Nevertheless, no 
significant alteration of the curve shape with 8 could 
be detected over the range covered (0:21 < 8<0-81). 

Again, in Ref. 2 observations are reported on the 
eddy formation behind several gauzes. It was found 
that at very low Reynolds number regular eddy forma- 
tion behind the wires of a gauze did not occur; at higher 
Reynolds number eddies were cast off with regular 
periodicity. These two states presumably correspond to 
the observed flow pattern changes round isolated circu- 
lar cylinders for which, at low Reynolds numbers, two 
vortices form just behind the cylinder and remain there. 
while at higher Reynolds numbers vortices are cast off 
alternately from opposite sides, forming the well-known 
“Karman trail.” It was found in the tests of Ref. 2, 
made on six gauzes, that the “eddy critical” Reynolds 
number, at which the change-over from one regime to 
the other took place, varied directly with 8. This corre- 
lation would be expected to apply also to the variation 
of k with R—that is, the curves of log k would be ex- 
pected to be similarly positioned when plotted against 
log (R/£) instead of log R. However, plotting shows 
that this is not the case (see Fig. 2). 

The overall variation of loss coefficient over the 
range of Reynolds number covered is 2:5: 1. By analogy 
with the drag coefficient curve for an infinite cylinder. 
the log k curve may be expected to continue to rise as 
log R is further reduced, and it may also be expected 
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04) TABLE III 
1 MEASURED VALUES OF k AT R= 100 
| Source k 0-71 2)/ R 
Q 02 Ref. 5 0°277 8:22 8:53 a 
0-413 4:06 3-46 
0°672 0°81 0:86 
0°743 0:57 0-58 R 
Ref. 6 0-429 3°31 3-15 
0°475 2:45 2-44 
will 0-585 1-47 1:36 
20 30 0-649 0:79 09 
Ref. 4 0318 6:10 | 
FiGuRE 4. Possible extension of curve of Fig. 2 to higher ; 0-333 5°66 5-68 If 
Reynolds numbers. 0-351 5-12 5-05 m 
0-389 4:12 399 th 
. 0-451 2°99 2°78 CC 
that large variations may occur at considerably higher 0-490 3.94 7-94 fi 
Reynolds numbers, despite the fact that the curve flat- 0-649 0-79 0-98 3 
tens out at the upper end of the range covered. In Fig. 0-682 0°81 0-81 ke 
3, the cylinder drag curve is plotted over a wider range Ref. 2 “0212, 13-21 *s 
of Reynolds number, and it will be seen that after 0436 7-74 
flattening out in the range of R between 100 and 250, 0-494 1:71 
the curve drops again at about 350, then rises, levels off, see es 
and finally drops sharply at R ~ 2 x 10°. Some idea of 


the actual trend of the gauze pressure drop curve at 
Reynolds number up to about 4.000 may be obtained 
from three measurements made by Simmons and 
Cowdrey®? on a grid of rods 0:25 in. diameter at 0-9 in. 
pitch. The observed values of k are given in Table II. 


TABLE II 
MEASUREMENTS ON A COARSE GRID OF 0°25 IN. 
DIAMETER RODS"), (3 = 0-522) 


Air speed V 


ft./sec. R k 
1310 1-34 
20 2620 1-45 

15 


30 3930 


In Fig. 4, log k — log kx~,,,, is plotted against log R:; the 
method given in the next Section was used to estimate 
k at R=100. A curve somewhat similar to the cylinder 
drag curve can be used to link these points to the 
mean curve established up to R=600, as shown by 
the dotted line. 


5. Variation of Pressure Drop Coefficient 
with Porosity 


Using the general curve of Fig. 2 to extrapolate 
where necessary, it is possible to pick off for each gauze 
the value of A corresponding to a given Reynolds num- 
ber. Table III gives the values so obtained at R= 100. 

Preliminary plotting showed that the values of k so 
obtained from the results of Ref. 2 were consistently 
lower than the others. It seems probable that this is 
a genuine difference arising from the different experi- 
mental arrangement used. Absence of constraint on the 
downstream side of the gauze may have affected the 
flow pattern, but it seems unlikely that this could have 
had much effect near the centre of the gauze; a more 


0-810 0°38 


likely explanation is that the method of measurement | 
of the pressure on the downstream side of the gauze was F 
at fault. The gauze under test was supported on flanges | 
attached to the duct outlet; it is known that a jet issuing 
from a duct with a flanged outlet into still air produces 
a small depression across the outlet, and it is probable 
that this depression is not uniform across the section. 
These results therefore were excluded from the investi- 
gation of the variation of k with /. 

Some assistance in deciding on the form of the rela- 
tion between k and / to be used as a basis for analysis | 
can be obtained from simple theoretical ideas. Consider} , 
a stream filament which just passes through a single — 
mesh. It suffers first a contraction down to the “ throat” F 
of the mesh, during which its mean velocity is increased | 
to V/; if it is assumed that losses during this process | 
are negligible and that the stream breaks away from the | 
wires at this section and then expands freely to the | 
original area, the stagnation pressure loss during this | 
process being that corresponding to a sudden exponents 
through the area ratio (1//), it is easily shown that the | 
loss coefficient k is equal to 8)°/ 5°. In practice, the 
flow pattern will not be as simple as this. There will be 
friction losses at the wire surface, with the formation of | 
boundary layers which will increase the effective con-~ 
traction, and the point of breakaway probably moves | 
round the wire as Reynolds number varies. However. 
it seems reasonable to investigate first the suitability of | 


TABLE IV 
ISOLATED MEASUREMENTS FROM REF. 1 


k corrected 


B Test R k toR=100 0-71 
0-372 102 4-67 4-68 441 
0-469 165 2-10 2:37 


0°554 252 1:20 1:47 1:60 
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the empty duct. 
contain a factor (1-8), and it may be written 
k=all 


The numerator therefore 


8)(1-cB)/8*. By plotting k8?/(1 — 8) against 
8. it was found that the best mean curve of this form is 


given by 


0-68 (1—8)(U+1-1 8) 


k 


- Calculation showed, however, that over the range of the 
~ tests (8 <{ 0-27) the values of k given by this formula are 
not significantly different from those given by the 
simpler equation 


0:71 (1 6°) 

In Tables III and IV and Fig. 5, this relation is com- 
pared with the measurements, including some isolated 
measurements from Ref. | not used in deriving the 
empirical formula. Agreement is quite good. The 


k= 


_ curve falls slightly high at £ less than 0°35 and slightly 
_ low for 8 between 0:35 and 0-45. 
could probably be obtained by fitting a more elaborate 


Better agreement 


equation containing higher powers of £, but it has not 


_ seemed to the author that the accuracy of the experi- 
mental data would justify this elaboration. 


| | 
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FiGuRE 5. k plotted against at R= 100. 
Curve drawn is k=0-71 (1—£?)/ 2. 


w. J. ANNAND THE RESISTANCE TO AIR FLOW OF WIRE GAUZES 
TABLE V 14 
VALUES OF a FOR VARIOUS VALUES OF R | | | 
2 30 40 60 7 80 90 
a 130 118 1:09 0:97 0°89 0:84 080 0°76 0°73 
R 100 200 300 400 500 600 . Laced 
the slightly more general formula k =a (1 + + cB*)/ “Tt 
If the formula is to be usable for gauzes of high it L 
must give the value k=0 for 6=1-0, corresponding to Los ® 


Ficure 6. Variation of multiplying factor a with 
Reynolds number. 


Some indication that the formula can be used down 
to 8=0-21 is given by the fact that a similar formula 
with a lower coefficient fits the results of Ref. 2 satis- 
factorily down to this value. 


6. Working Method of Estimation of k 
The working formula for the estimation of the 
pressure drop coefficient is then 
a(1— 


where a depends on the Reynolds number based on wire 
diameter, as shown in Fig. 6 and Table V. 

To assist in practical application, the nomogram 
given in Fig. 7 has been prepared. From this chart the 
multiplying factor a can be obtained directly from the 
values of the wire diameter, air speed, air temperature 
and air pressure. The method of use is indicated by the 
dotted lines; join diameter and air speed and find inter- 
section on reference line I; join air temperature and 
pressure and find intersection on reference line II; join 
the two intersection points and read off a. The point on 
reference line II corresponding to I.C.A.N. standard sea 
level conditions is ringed. 

It should, perhaps, be emphasised once more that 
the curve of variation of a with Reynolds number 
should not be extrapolated outside the range of the test 
data. The formula can be used with confidence for £ 
not less than 0:27 and is probably adequate down to 
8=0-20; below this value, there is no evidence as to the 
adequacy of the formula, and testing would be advisable. 

One other point concerning practical application 
should be remembered. That is, that accurate estima- 
tion of 8 and of k at low values of £ calls for very 
accurate measurement of the gauze dimensions. Table 
VI gives the percentage error in k produced by a one 


TABLE VI 


EFFECT ON ESTIMATE OF k OF ONE PER CENT. ERROR 
OF WIRE DIAMETER 


B 02 03 04 05 06 0-7 08 


Percentage error of estimated k 5:2 36 2°8 2:2 18 1:5 1:3 
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400 — 
4 2 
DOTTED LINE SHOWS METHOD’ OF USE 
4 a FOR TEMPERATURE = IS°C, PRESSURE = 6 LB./IN®, 
O56 + — SPEED: 40 FT./SEC., DIAMETER = 0-007 IN, 
0-09 THUS GIVING a =0-83 
0:08 
7 0-07 
O58 
loo 0:06 
= 90 0:05 
a 0-04 
w 30— \ et. w of the 
50804 a a = zt © 
= --0-008 find e 
110 0.007 20 {their 
4 =i... 

1-20 — 0006 the A 

1-30—4 0:005 40 travel 

100 stor 
0-003 
s— TEMPERATURE 
AND PRESSURE by the 
Figure 7. Nomogram for determination of multiplying factor a. interp! 
Eig! 
per cent. error in measurement of the wire diameter d 2. SCHUBAUER, G. B., SPANGENBERG, W. G. and KLEBANOFF, F pan. 4) 
at various values of 8. At 8=0-3, for example, one per P. S. 
cent. error in d produces 3-6 per cent. error in kK. Now 
one per cent. error in d may, for the wire sizes commonly 3. ScHuBAUER. G. B. and SPaANGENBERG, W. G. (1947). The # 

Ps . Effect of Screens in Wide-angle Diffusers. N.A.C.A. T.N. travell 
used in gauzes, be less than 0-000] in., so that the need 
f 1610, June 1947. OF 
or accuracy is evident. apt: 

y 4. Simmons, L. F. G. and Cowprey, C. F. (1945). Measure- able bi 
Provided that the limitations are observed and the ments of the Aerodynamic Forces Acting on Porous — the air 
necessary accuracy of gauze measurement obtained, the Screens. R. & M. 2276, August 1945. birds j 
method, according to a rough statistical estimate, should 5. Eckert, B. and PrLuGer, F. (1941). Resistance Coefficients F Wrote, 
give the pressure drop within 20 per cent. on 95 per cent. of Commercial Round Wire Grids (translated from Luft 5 day wt 
of occasions. fahrtforschung, Vol. 18, No. 4). N.A.C.A. T.M. 1003, of the 
January 1942. Britain 

6. Taytor, G. I. and Davies, R. M. (1944). The Aerodynamics the La 
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Experiments With Flat Planes 


(1866-1890) 


J. LAURENCE PRITCHARD 


(From time to time we hope to be able to publish extracts from the notes for the History 
of the Roval Aeronautical Society now being prepared by Captain Pritchard) 


"THE PIONEER rarely receives that meed of praise 

which is his due, until those days arrive when his 
pioneering has become a page or so of history. What 
meed of praise he may then receive is too often made 
savourless by critics unable to disassociate themselves 
from the knowledge of their own technical era. The 
halting speed of pioneering aerodynamic research can- 
not be judged by the supersonic speed of secret splittings 
of the atom. 

It was R. L. Stevenson who said that to travel hope- 
fully is a better thing than to arrive. In aviation we 
have arrived at that streamline flow of technical jargon 
which deludes so many they have at last reached the 
boundary layer of aeronautical knowledge, so that they 
find experiments with plates always a little flat. although 
their curiosity may still be aroused by flying saucers. 

Eighty-seven years ago there were many members of 
the Aeronautical Society of Great Britain who were 
travelling hopefully along roads which were badly lit by 
the first few glimmers of the dawning of science. And 
what fun they had in their Victorian travels on those 
dimly lighted roads! And what fun an amateur 
historian must have trying to travel for a little while 
with them before the roads are completely obliterated 
by the jargon of jets and the moonshine madness of 
interplanetary travelogues. 

Eighty-seven years ago? That was the year the 
Society was founded, 1866. Those who became mem- 
bers then had travelled alone, but from that year they 
were able to travel in company. 

It was in the 1850's that the great W. H. Wenham 
travelled up the Nile and there studied, from the point 
of view of an engineer, the slow flight of that unbeliev- 
able bird, the pelican; and that graceful ballet dancer in 
the air, the spoonbill; and the flight of a dozen other 
birds indigenous to the Egyptian air. On his return he 
Wrote, in 1859, that classic paper which saw the light of 
day when he read it at the first lecture meeting in 1866 
of the newly formed Aeronautical Society of Great 
Britain. In his paper, entitled “ Aerial Locomotion and 
the Laws by which Heavy Bodies impelled through the 
Air are Sustained,” Wenham remarked on_ the 
importance of a rigid leading edge; on the fact that the 
greatest part of the lift on a plane was in the first eighth 
to tenth part of the wing; and that man could never 
hope to fly under his own power. He advocated the 
long narrow plane, and to avoid great spans suggested 
the use of super-imposed planes. 
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Wenham, during the discussion following his lecture, 
drew attention to the work of Sir George Cayley, to 
whom as early as 1846, twenty years before the founda- 
tion of the Society, W. S. Henson had written on behalf 
of himself and John Stringfellow, “We therefore 
resolved to apply to you as the Father of Aerial 
Navigation . . .”, a description which has been more and 
more emphasised during the century which has 
followed, as the work of Cayley has become more 
widely appreciated. Wenham reminded his listeners 
that Cayley had in 1810 drawn attention to the move- 
ment of the centre of pressure at high speeds and small 
angles of attack. 

Indeed, Cayley, Wenham and others were well 
aware of the kind of data that was required, but 
experiment and research were painfully slow. F. W. 
Lanchester in his Aerodynamics, published in 1897, 
gives a brief account of some of the earliest attempts 
at aerodynamic research. To quote, “The earliest 
experimenters, Robins (1761), Hutton (1787) and 
Vince (1794-95, 1797-98) employed a primitive form of 
whirling table, the invention of which is attributed by 
Hutton to Robins.” 

These results and experiments of others were known 
to the early members of the Society but they led, in 
certain respects, to errors of belief about the pressures 
on flat plates at various angles of attack. Over a 
hundred years after Robins, the problems raised by flow 
of air over plates and bodies of different shapes 
remained unsolved and on 24th August 1866 the 
Council passed the following remarkable resolution : — 

“Tt was delegated to the Duke of Argyll, Lord 
Dufferin, Mr. James Nasmyth, Mr. J. McConnell, Mr. 
Glaisher and Mr. Wenham, to consider the law of 
resistance of inclined surfaces moving in elastic and 
non-elastic fluids, as air and water, and the resultant 
force obtained at right angles to the direction of motion. 
Experiment has now proved that as the velocity is 
increased, and the angle of inclination diminished, the 
centre of effort of the plane approaches forward, till 
nearly the whole effect of reaction or support is trans- 
ferred to the front edge. It therefore follows that in 
very acute angles with the line of motion, a given extent 
of reacting surface at high velocities will be most effec- 
tive if disposed in a narrow plane extended transversely 
to its direction. 

“ An investigation of this principle will determine the 
best ratio of surface to speed in the long narrow wings 
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L. PRITCHARD EXPERIMENTS WITH 


of birds of extended flights; in the blades of reciproca- 
ting and screw propellers; and in the sails of vessels for 
sailing near the wind.” 

The resolution was not acted on as a matter of 
urgency. After all, it was the Victorian age, and there 
were more important things to do than to follow the 
flying will-o’-wisp ! Some idea of the attitude of mind 
generally of the public at the time, an attitude which 
inevitably had some effect on how urgent the suggestion 
might be, may be obtained from the words used by 
Stewart Harrison on 15th July 1867. He was proposing 
a vote of thanks to the Duke of Argyll and Mr. Hurry 
for taking the Chair at the lecture meeting of the Society. 

“T have no doubt,” he said, “that the people out- 
side thought it very desirable that those present should 
be taken care of by their friends. I do not, but I do 
think that, as gentlemen who risked something of their 
reputation for soundness of judgment by taking the 
Chair when such flighty subjects were discussed, those 
present could not do less than acknowledge their kind- 
ness by a vote of thanks.” 

Meetings were not so frequent in those days as they 
are now, nor were the opportunities there for pressing 
forward with ideas about flying. However there were 
members of the Society always ready to keep the objects 
of the Society in the foreground on those annual 
occasions when a lecture was held. At the May meeting 
in 1869 a Captain Murray, R.N., said he thought the 
ciety should offer a prize of £100 “for a complete and 


_,eneral solution, on rigid mathematical principles” of 


| the following 
,; between the velocity of a current of air and its pressure 


problem:—‘“ Required relation 


on a plane surface of given size, shape and inclination.” 
Captain Murray offered to subscribe £10 towards the 


fund. 


At the same meeting a member, E. W. Young, said 
he had made some rough experiments upon the lifting 
effect of the wind upon an inclined surface. He had 
found that with a vane inclined at an angle of 12° to 
the horizontal the vertical lifting effect was about one 
half the pressure which it exerted against the same vane 
placed vertically. 

Almost exactly a year later the Council discussed 
again the absence of data and resolved that an Experi- 
mental Fund be set up and an Experimental Committee 
be appointed. In the Annual Report for that year, 
1870, it was stated, “ The first great aim of the Society 
is the connecting of the velocity of the air with its 
pressure on plane surfaces at various inclinations.” and 
continued, “The Aeronautical Society has _ been 
endeavouring for some time to inaugurate a systematic 
series of experiments on the connection between the 
pressure and velocity of air, and it is believed that these 
will afford the only data on which a true science of aero- 
nautics can be founded. For this purpose an instrument 
has been designed by Mr. F. H. Wenham, and approved 
by the Experimental Committee, which it is intended to 
subinit to the action of a fan not less than 30 in. in 
diameter, capable of delivering about 3,000 cubic feet 
of air per minute. A clear space of 15 feet or more in 


‘fron: of the fan will allow room for a square wooden 
‘tunnel to guide the blast, ascertain its velocity, and 
| 
| 
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insert the anemometer. Mr. John Browning, F.R.A:S.. 
will undertake the construction of the instrument, and 
Sir William Fairbairn, F.R.S., will afford all the aid in 
his power for its submission to the required test.” 

These were the first experiments to be carried out 
by the Society, and the first wind tunnel experiments for 
aeronautical purposes in the world. Sir William Fair- 
bairn was a famous engineer and, like Browning and 
Wenham, a member of the Council. The experiments 
were carried out at the engineering works of Messrs. 
Penn at Greenwich. 

The results of these experiments of over eighty years 
ago, were published in the Annual Report of the Society 
for 1871. They were not very satisfactory, for the 
sources of error were great and largely unknown to the 
experimenters. Angles of incidence below 15° were 
not tested and the conclusion was drawn that the lift 
varied with the sine of the angle of incidence. But 
from these experiments the Council believed that 
the power for flying was not so great as had been 
expected. The horse power of many of the larger birds 
had, indeed, been exaggerated. Wenham always 
said that the strength of a bird was less than many 
had asserted and said in his 1859 paper of the pelican, 
“it displays no symptom of being endowed with great 
strength, for, when only slightly wounded, it is easily 
captured, not having adequate power for effective 
resistance, but heavily flapping the huge wings, that 
should, as some imagine, give a stroke equal in vigour 
to the kick of a horse.” 

“Tt is desirable that these important experiments 
should be verified and continued,” added the Annual 
Report, and continued, “Up to the present time only 
flat surfaces have been experimented upon, and as all 
the sustaining surfaces in the wings of birds are curved 
or hollow, it remains to be proved what is the relative 
advantage of this form, for it is upon such data that all 
plans of construction must be based, and on which 
failure or success must depend, and to determine 
whether flight is practicable or not.” 

At the Annual Meeting held in 1872, James Glaisher, 
from the Chair, drew attention to the difficulties of 
measuring lift and drag and said. “ Anyone who had 
not considered with some care the nature of the pressure 
produced by the flow or rush of a fluid, elastic or incom- 
pressible, against a plane surface placed in its course, 
might imagine that the system of parallel forces was 
merely equivalent to a single resultant force acting at 
the centre of pressure, and capable of resolution accord- 
ing to the ordinary parallelogram law. But this of 
course is not the case, for the particles of the fluid which 
come in contact with the plane have somehow or other 
to get out of the way, by gliding along the surface of 
the plane (as they cannot get through it) and this pro- 
duces a complication in the neighbourhood of such a 
kind as cannot theoretically be predicted. One thing. 
however, is quite clear, and that is that the directions 
of all the small forces acting on the surface certainly 
are not parallel, and we must therefore have recourse to 
experiment.” 

The problem of the boundary layer and the difficul- 
ties of turbulent flow were then foreshadowed. 
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Research without the powerful incentive of the spirit 
or of necessity moves but slowly. The experiments of 
Wenham and Browning had attracted considerable 
attention, and Thomas Moy who, later, was himself to 
continue the experiments, said at the 1874 meeting that 
he hoped some clever geometrician or mathematician 
would construct a formula from Wenham’s results and 
demonstrate the true law of air resistance. At that date 
mathematicians had not been so well tamed as they are 
now. 

In May 1877 Thomas Moy and another heavier- 
than-air enthusiastic member of the Society, R. C. Jay. 
appealed to the Council for a grant from the Experi- 
mental Fund for apparatus for further experiment. It 
was not until a year later that their proposal was 
approved by Wenham on behalf of the Council and a 
grant of £20 made from the Fund. The proposal, in 
Moy’s own handwriting, is reproduced on page 148. It 
is one of the treasured possessions of the Society. It 
will be noticed that Moy estimated the cost at £15 and 
the Council granted £20. Financial control, inter- 
national. national or local had not then raised its 
octuple tentacles to curtail the enterprise of the 
perishing pioneer. 

Moy’s report was submitted to the Council on 4th 
April 1879. It read as follows :— 

“The apparatus answers its purpose very well up 
to a speed of 21 miles an hour. Beyond this speed it is 
very difficult to turn the handle fast enough; and it will 
be necessary to have a smaller grooved pulley on the 
upright shaft in order to increase the speed and obtain 
pressures above 20 miles an hour. 

“ The series of experiments tried during the autumn 
and winter of 1878 were for the purpose of finding the 
centre of pressure upon inclined rigid planes, as it had 
always been found. in practice, that the centre of 
pressure approached the front edge of such planes when 
acted upon by the air. 

“Numerous. planes were tried. various 
expedients adopted in order to work out this very 
delicate operation. At last, after many attempts and 
failures which I need not describe, a thin wooden plane 
4 inches wide by 12 inches long was mounted with an 
axle which could be made to compensate for any 
unequal weight existing between the fore and aft 
portion; a loose pin held it at the angle experimented 
upon, and upon the pressure on the after portion exceed- 
ing that of the front portion the pin was released and 
fell down. These experiments proved the very 
important fact that from the vertical to an angle of 45 
the centre of pressure was at the centre of the plane: 
but that from 45° to 10° the centre of pressure shifted 
forward to as much as 15 per cent. of the whole plane— 
at 30°. 20°. 15° and 10° the centre of pressure became 
uniformly 35 per cent. from the front edge, and 65 per 
cent. from the after edge. 

“T fully expected that varying speeds would give 
varying results as to this, and was suprised that I could 
not obtain them: I hope, however, soon to resume the 
experiments in order to obtain more light upon this 
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interesting and important point, and also to prove the 
results obtained by Mr. Browning’s anemometer and to 
extend the experiments to finer angles than were tried 
at Messrs. Penn’s establishment. 

“In conclusion I may add that flexible planes will 
form an important series in the experiments to be 
tried.” 

The actual report, in Moy’s handwriting, is filed in 
the Minute Book of the Council of that year. 

Of course these experiments were crude compared 
with present-day methods, and were not as accurate. 
But it took a major war to get better answers, and, full 


circle, war has made demands which have left members | 


of the Society urging more and more should be spent 


on aerodynamic research, for what Moy failed to find | 


with £15 with planes moving at low speeds, present wind 


tunnels have not succeeded in finding with £15.000,000 — 
Eighty years ago what members of the | 


at high speeds. 
Society were lacking in data and information on aero- 


dynamics they made up for in enthusiasm and ad hoc! 
attempts with flying models. Now too much power has | 
overshadowed the lack of information of aerodynamic | 
problems at high speeds, and supersonics have made us | 
deaf to the advantages of the slow but sure method of | 


progress. 

Let me finish with a quotation from that great 
American member of the Society, Octave Chanute, the 
friend of the Wright Brothers, writing in The Railroad 
and Engineering Journal of New York for October 
1891, 

“T had a chart plotted on which were delineated all 
the experiments on inclined surfaces which I could learn 
about—those of Hutton, Vince, Thibault, Duchemin. 
de Louvrie, Skye, the British Aeronautical Society and 
W. H. Dines: and on this chart I also had plotted the 
curves of the various formulas. 
great discrepancies, yet by patient analysis various 
probable sources of error were eliminated, and the con- 
clusion was reached that the formula last given, which 
I have seen variously attributed to Bossut or to 
Duchemin, was probably correct.” 


That formula, for those who do not know it, is 


1+ sin* x 
where P is the pressure on the inclined surface, P’ the 
pressure on the surface when at 90° to the air flow, and 
z is the angle the inclined plane makes with the current 
of air. But of all these investigations, up to the time 
Chanute wrote in 1891, W. H. Dines, in my opinion. 
showed the best experimental methods and _ results 
Past that date I do not feel inclined to go for the 
moment, for the very air above our heads has become 
secret and to talk about it may lead to a turbulence o/ 
the soul beyond our control. No longer can we say. 
with our national poet (habit almost made me sa\ 
nationalised poet): 


“The isle is full of noises, 
Sounds and sweet airs, that give delight and hurt not. 
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bers | 
pent | SUMMARY:—The problem of the stress analysis of circular fuselage frames has been 

i investigated by a number of authors: the analyses are, necessarily, rather mathematical 
find : in nature although in some cases the final results have been presented in extremely 
vind practical form. The present paper offers, in non-mathematical language, a short account 
000 of the fundamentals of the problem and a brief guide to the more important published 
the | literature. It also develops an approximate method for the analysis of frames using 
ero- Lagrange’s method for minimising a function of several inter-dependent variables. This 

hoc? method is shown to be particularly suitable for the analysis of frames with large cut- 

- has outs, about which little has previously been published. In addition, the underlying 
amic parameters which determine the stress distribution are deduced from the theory, and 
the deductions are compared with previously published work. In general, agreement 
le us is good. The paper is written throughout with a view to being of practical use in the 
dof actual stressing of frames. Numerical examples and explanations of mathematical 
methods are included in Appendices. 
oreat 

the 
road 1. Introduction g number defining section of frame 
paper falls naturally into three parts: — MOMents of 
dalll 1. A review of the whole field of currently accepted 
learng {tame theory, together with brief notes on some of the 

; more important of published investigations. : 
ss: ll. The development of an original method of —_ 
y and P nig K.K, constants 
d ‘the analysis. using Lagrange S undetermined multipliers for K — L1/(2UG) 
bited | Minimising strain energy, the method being particularly K.—1/QE) 

adapted to the solution of frames with a large cut-out. 
— Ill. A discussion of the underlying parameters L assumed length of fuselage in 
con- i which distorted shear reverts to 
which? hich determine the stress distribution, using the results Cane 14400) 
<a of previously published work and predictions from the L (Section 19 only) average frame 

theory developed in Part II. spacing 

is length of one section, measured in 

Notation skin line 

1) (or 6 1) for Ref. M bending moment 
4. See Section 19) M, applied couple 
P’ the *A_ (Ref. 5 only) area of cross section M, free bending moment in floor beam 
of frame. See Section 20 (Appendix IT) 
urrent “BO Err rr M, redundant moment (Appendix 1) 
age for Ref. 4. See Section 19) M,, bending moment at centre of m 
iy: *A/B=Gitr'/(EIL) section of frame due to shear in 
er C=1/Gr’) skin in sections | to m 
esults ES/(GS’) M,,’ bending moment at centre of 
or the a,,, product of /and “ offset from centre section due to external applied 
ecome of m'" section to tangent at skin forces P. T or M, 
NCE of line of n'" section.” See Section 14. m number defining section of frame 
ve Say: C,.d,, Offsets defined and used in Appen- M1,.Ms....M, coefficients used in Appendix II 
ne dix V N internal redundancy (axial load) in 
c.d.¢ parameters used in Refs. 1, 8 and Appendix I 
referred to in Section 19 n number defining section of frame 

t not. E Young's modulus P applied radial load 
shear modulus 

Receis cd September 1952 *For A. A B and B, L is to be taken as mean frame spacing. 
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p number defining section of frame 
Q, defined and used in Appendix III 
q_ shear per unit length, or “shear 
flow” in skin 
Rs radius of skin centre line in Ref. 1. 
5. 8 (see Section 19) 
r radius of skin centre line 
r constant used in Ref. 8 (see Section 
19) 
r’ radius of centre line of frame, if 
different from r 
S cross-sectional area of frame (i.e. 
end load carrying area) 
S’ cross-sectional area of frame web 
(i.e. shear carrying area) 
t skin thickness 
t’ thickness of equivalent end load 
carrying sheet replacing skin and 
stringers 
U strain energy 
V internal redundancy used in Ap- 
pendix I 
2W distributed floor load in Appen- 
dix V 
internal redundancies used in Ap- 
pendices II, V 


PART I 
A REVIEW OF CURRENTLY ACCEPTED FRAME THEORY 


2. Statement of the Problem 


When one of the circular frames of a cantilevered 
monocoque fuselage of constant diameter is loaded (see 
Fig. 1), the load is reacted by shear in the skin, together 
with associated end loads. The reacting shear gives 
rise to stresses in the frame, i.c. bending moment, axial 
load and web shear. For a complete stress analysis all 
of these must be determined. 

The solution falls naturally into two parts: (a) the 
determination of the distribution of the reacting shear 
in the skin; (b) the determination of the bending 
moment, et cetera, in the frame. If the skin shear dis- 
tribution is known, (b) consists fundamentally of the 
evaluation of two redundancies (see Section 9). If the 
skin shear distribution is not known, the values of these 
redundancies are generally obtained as a by-product 
from the solution of (a). (>) is thus a relatively simple 
problem, but (a) is in general highly complex. The 
determination of the skin shear distribution constitutes 
the crux of the problem. 


3. Types of Circular Frames 
Circular fuselage frames can be (arbitrarily) classi- 
fied into five types. all of which are of common 
occurrence in practical aeronautics :— 
(i) Complete frames of constant cross section, 
with no internal reinforcing. 
(ii) As (i) but with internal reinforcing in the form 
of floor beams and supporting struts. 
(iii) Frames with a cut-out. 
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X,,X, (Section 17 only) restraining forces 
from wing 


flow’ 
...n'" sections of frame 


x distance along floor beam (Appen- © 


dices II, V) 
via sft angles defining points of frame 
A, undetermined multipliers for use 

with Lagrange’s method of mini- 


shear per unit length or “shear © 
in skin at first, second,...p", 


mising a function of several inter- 


dependent variables (see Appendix 
VI) 
» eccentricity =(r—r’)/r 
Poisson’s ratio 
A, = (maximum bending moment, shear 
flow, etc.. for a given value of A)/ 
(maximum bending moment, shear 
flow, etc., for A > 10°) 
X, = (maximum bending moment, shear 
flow, etc., for a given value of 1))/ 


(maximum bending moment, shear | 


flow, etc., for =0) 


€ used and defined in Appendix VII.) 


(iv) Frames with a variable cross section. 

(v) Complete or partial bulkheads. 

Unless otherwise stated. Sections 5 to 10 refer to 
type (i). Types (ii) to (v) are dealt with in Sections I1 
onwards. 


4. Types of Frame Loading 
The following types of loading are considered : — 
(i) Radial loads applied direct to the frame. 
(ii) Tangential loads applied direct to the frame. 
(iii) Moment loads applied direct to the frame. 
(iv) Loads applied to the frame through a_ floor 
beam. 
Other systems of loading can be treated as a com- 
bination of these four. 
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Ficure 1. Cantilevered monocoque fuselage loaded 


at frame A. 
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Fict ke 2. Radially loaded 
frame showing skin shear 


reactions. frame). 


5. Distribution of the Reacting Shear 

When the applied load is radial, simple engineering 
bending theory requires that the reacting shear should 
be distributed sinusoidally; i.e. at a point A in Fig. 2 
defined by an angle 6 measured from the diameter 
through the point of application of the load, the shear 
per unit length or “ shear flow” q is given by g=K sin @. 

The value of K is readily obtained by considering 


.) overall balance in a direction parallel to the applied 


load. 
Hence 2 | K sin? 6rdé=—P 
vielding K- = and q= sin (1) 


If a tangential load T or a moment load M is applied 
at B (Fig. 2), simple engineering theory similarly gives 
the shear flow at A as 


g= ra +cos 4) : (2) 


respectively. These shear distributions are illustrated 
diagrammatically in Figs. 3, 4 and 5. 

Engineering experience suggests that in the neigh- 
bourhood of the loaded frame there may be some 
distortion of the pattern of the skin shear distribution. 
Both experiment and calculation confirm this, if the 
loaded frame is flexible. Moreover the distortion of the 
shear pattern gives relief to the bending stresses in the 
frame. If. however, the frame is rigid there is no 
distortion of the shear pattern. 

The terms “rigid” and “ flexible * are comparative. 
not absolute. and require definition. A suggested 


definition will be found in Section 14(j). 


When there is variation of skin thickness round the 
periphery of the loaded rigid frame, the skin shear dis- 
tribution given by equations (1)-(3) will be modified. 
This modified shear can still be calculated by engineer- 
ing theory (the shear stress being given by some 
expression such as FSy/(7B) in the notation of Case, 
Strenth of Materials, p. 202). But it is still the “basic” 
Shear. Distortion of the skin shear pattern due to this 


OF CIRCUL 


Figure 4, Skin shear distribu- 
tion due to tangential load 


Ficure 8. Skin shear distribu- 
tion due to radial load (rigid 


Ficure 5. Skin shear distri- 
bution due to moment load 
(rigid frame). 


(rigid frame). 
cause must not be confused with distortion due to 
frame flexibility. 


6. Nature and Magnitude of the Deforma- 
tion of the Skin Shear Distribution 
Pattern (and Bending Moment Pat- 
tern) at a Flexible Loaded Frame. 


Figures 6 and 7 show typical skin shear and frame 
bending moment distributions for a flexible frame. 
loaded radially, compared with similar distributions 
for a rigid frame. The figures are diagrammatic only 
and represent the type of result obtained by various 
investigators. Actual curves can be seen in the quoted 
references. 

The ratio q,/q. in Fig. 6 gives a measure of the dis- 
tortion of the shear pattern. Hoff ” calls it the “ shear 
concentration factor” and examines its variation with 
parameters representative of frame and skin stiffnesses. 
Hoff also plots the “moment reduction factor ” (repre- 
sented by m,/m, in Fig. 7) against the same parameters. 
Similar results are obtained with tangential and moment 
loadings. 

For the most flexible frame examined Hoff found 
shear concentrations as high as 35, and moment reduc- 
tion factors as low as 0:33 (see Ref. 4, Figs. 5 and 6). 
Values outside these limits could occur with still more 
flexible frames, but practical considerations will put a 
limit to the flexibility of frames actually used in 
aircraft structures. 


7. Reversion of the Distorted Shear to 
Basic Form 


If a flexible frame A well away from YX, the 
supported end of the cantilever (see Fig. 1). is loaded 
radially, the skin shear distribution in its immediate 
neighbourhood will be distorted in the manner of 
Section 6. But at X the shear will have taken up its 
basic form again. The change occurs in steps at the 
intervening frames B,C, D, E,...... The number 
of frames required to “revert” the distorted shear 
pattern to basic depends on their flexibility. One 
fairly rigid frame might suffice. With very flexible 
frames, as many as six might be required. (See Ref. 8. 
Fig. 7; Ref. 11. Fig. 4; Ref. 2, Fig. 4.) 

Stresses will be set up in the succeeding frames 
B,C, D, and so on, depending on the difference in the 
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FIGURE 6. Typical skin shear patterns for flexible and 
rigid frame. 


Shear distribution on the two sides of such a frame. 
These stresses are plotted for various cases in Refs. 2 
and 7. In general they are small, but not necessarily 
negligible. 

If no intermediate frames occurred, there would still 
be a tendency for the shear to revert to basic before X, 
each ring of skin acting as a very flexible frame. There 
can, therefore, be a change of shear pattern without 
intermediate frames. But when frames are present, and 
not too widely spaced. strain energy considerations 
demand that this change shall occur chiefly at the 
frames. Hence in current frame theory. the shear dis- 
tribution is treated as constant between two successive 
frames. 


If a frame is loaded with tangential or moment 
loading, similar changes in the shear pattern occur. 
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4 Figure 7. Typical bending moment distributions for flexible 


and rigid frames. 
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8. Factors that Determine the Shear 


Distribution 

These are: — 

(i) The second moment of area of the frame cross 
section (/), 

Fuselage skin thickness (f), 

(iii) Radius of fuselage (r), 

(iv) Frame spacing, 

(v) Eccentricity », i.e. offset of centre line of skin © 
from centre line of frame. expressed as a frac- 
tion of r, the skin radius, 

(\7) Stringer area and spacing and, to a less extent, 
bending and torsional stiffness of stringers, 

(vii) The number of “fields” between the loaded 
frame and the point of support, 

(viii) The number of “fields” between the loaded © 
frame and the free end of the fuselage. 

(ix) The extensional and shearing rigidities of the | 


frame. 


The term “field” is to be taken as indicating the 


portion of fuselage between two successive frames. 


The relative importance of these factors is discussed’ | 


in Part III. 


9. Determination of the Bending Moment) 


in a Complete Frame of Constant 
Cross Section with a Given Distri- 
bution of the Reacting Shear 
Figure 8 shows a frame loaded radially at a point A. 
An imaginary cut is made at a diametrically opposite 
point B. and internal forces N (axial load), V (shear) 
and-M., (moment) are applied to the cut edges. 


If the frame is symmetrical about AB, the shear dis- | 


tribution will also be symmetrical. Then, from sym- 
metry. V is zero, and the problem resolves itself into 


the determination of two internal forces only, N and 


M,. These constitute the redundancies of the system. 
If a tangential load T or a moment load M,, is 
applied to the frame A, the redundancies again reduce 
to two in symmetric frames. (For 7, make the cut at 
0°. Nis P/2 and —P/2. Redundancies are M, and V. 
For M, make the cut at 0°. M, is M,/2 and - M,,/2. 
The redundancies are N and V.) 


One of the first published solutions to this problem 
was due to Wise“. Assuming the distribution of skin, 
shear given by simple engineering theory (see equations 
(1) to (3) of the present paper), he evaluated the 
redundancies for a radial, tangential or moment load, 
and hence calculated the bending moment, axial load> 
and web shear at any point in the frame. His results 
were published in the Royal Aeronautical Society's! 
Stressed Skin Data Sheets (03.06 series), but since he 
used a shear distribution based on engineering theory. 
his results apply to rigid frames only. 

Wise obtained his solution by using a “method of 
consistent deflections.” The problem is equally amen- 
able to routine strain energy analysis, and a solution 
by this method is outlined in Appendix I. 
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If the original shear distribution is given in 
griphical, not analytic, form, the strain energy method 
of Appendix I still applies, but the integrations must 
be performed graphically. This case could arise when 
the skin in the neighbourhood of a rigid frame is vary- 
ing in thickness round the frame, the shear being deter- 
mined by engineering theory but not being capable of 
analytic expression (see Section 5). In other cases. as 
already stated, determination of the shear pattern 
generally yields the values of the redundancies and 
further strain energy analysis is unnecessary. 

Morse“ suggests that the calculation of the bending 
moment from a known shear distribution is facilitated 
by applying the redundancies at the “elastic centre,” 
which he defines in terms of the elastic properties of the 
frame. This seems to simplify the calculation in certain 
cases. 


Figure 8. Radially 
loaded frame, cut to 
show redundancies. 


10. Determination of the Skin Shear Dis- 
tribution at a Loaded Flexible 
Frame 


The determination of the skin shear distribution at 
a loaded flexible frame is a matter of considerable diffi- 
culty. It has been attacked by several investigators 
from varying points of view, but the method most 
favoured has been to represent the shear distribution in 
the skin by a trigonometric series (e.g. the shear flow 
at a point defined by @ is taken as a,+a, sin@+ 
a, sin 20+ a, sin 30+ .). 
The strain energy of the loaded frame and adjacent skin 
is then written down in terms of the assumed shear flow, 
and minimised with respect to each of the coefficients 
and so on, in turn. Alterna- 
tively differential equations are set up which ensure 
compatibility between the skin and the frame, and these 
are solved by any convenient artifice. 

Analyses of this type are long and cumbersome but 
some of the investigators have presented their results 
in an eminently practical form. Others have not done 
so and it is difficult to make practical use of their work. 

Brief comments on some of the more important of 
the published analyses are appended. 
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REFERENCE | 

This was the first published report showing the effect 
of fuselage frame flexibility. By ensuring compatibility 
between skin and frame deflections the authors estab- 
lished a differential equation whose solution gave them 
the shear flow. No allowance was made for the effect 
of eccentricity or stringer flexibility under end load (see 
Sections 19 and 20 of the present paper). The report 
contains a discussion of underlying physical concepts, 
and nests of curves similar to those discribed in Ref. 7 
are presented. Unfortunately, these will be substan- 
tially in error in many practical cases, i.e. for low values 
of A (see Section 19 and Fig. 29 of the present paper). 


REFERENCES 2, 3 AND 4 

Hoff and his associates analysed frames using strain 
energy, the shear flow being represented by a Fourier 
series and the strain energy minimised with respect to 
the successive coefficients. Ref. 3 is an advance on 
Ref. 4 in that more of the factors affecting shear dis- 
tribution are included in the analysis, and allowance is 
made for radial reactions between the skin and the 
frame. Ref. 4 contains curves showing the effect of 
relevant factors on skin shear and frame bending 
moment. These are referred to again in Sections 19-21. 

Hoff deals only with radial loading and gives only 
frame bending moment and skin shear. His curves 
therefore do not have the same completeness as those in 
Ref. 7. Skin shear and frame bending moment can be 
calculated for tangential and moment loadings from 
the curves for the radial loading by using Goodey’s 
method (see Section 16 of the present paper), but axial 
load and web shear in the frame can only be deduced 
after the evaluation of two redundancies, as explained in 
Section 9. As the values of these redundancies must 
have been obtained at some time during the calculation. 
the practical utility of these curves would have been 
enhanced if their values had been plotted or tabulated. 
REFERENCE 5 

Beskin used a strain energy method. representing 
the skin shear flow by a Fourier series. He included 
the effects of frame eccentricity and also of shearing and 
extensional rigidities of the frame. He produced curves 
for frame bending moment. frame axial load and frame 
web shear for radial, tangential and moment loads, with 
parametric variation with skin eccentricity. These 
curves give the best available picture of the effect of 
eccentricity on frame stress distribution. They are dis- 
cussed more fully in Section 20 and a cross plot is 
presented in Fig. 30. 

REFERENCES 6 AND 7 

Duberg and Kempner? self-equilibrating 
systems of stress which they added to the basic stress 
system to give compatible displacements, relating these 
self-equilibrating systems successive bays by 
recurrence formulae. They did not allow for eccen- 
tricity. Hoff. in Ref. 3, claims that their solution agrees 
with his earlier solution. 

Reference 7 contains families of curves based on 
Ref. 6. These form the readiest available source of 
information on stresses in flexible frames. The skin 
shear and end load, the frame bending moment. axial 
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load and web shear are plotted round the frame for 
radial. tangential and moment loads with parametric 
variation against two factors representative of frame 
and stringer stiffnesses. There are some 45 pages of 
these nests of curves. The loading of the adjacent 
frame is also shown. The curves require correction for 
the effect of skin eccentricity (see Fig. 30). 

REFERENCE 8 

Goodey’s analysis is perhaps the most fundamental 
of the published analyses. He starts by considering the 
effects of all the factors that control the skin shear dis- 
tribution. including all those listed in Section 8 of this 
paper and also including the Poisson effect and the 
normal reaction between skin and frame. Later, in the 
interests of simplicity, he drops a number of these: 
notably he ignores stringer bending and torsion and 
uses zero eccentricity. 

Despite the fundamental nature of his analysis, his 
results have not been of great practical use owing to the 
complex mathematics and the paucity of published 
curves. Ref. 11 is a welcome step towards making his 
results more readily applicable. 


FiGuRE 9. Typical cir- 
cular fuselage frame 
with propped floor 
beam carrying distri- 
buted load 2W. 


REFERENCE 1] 

This report sets out to give curves based on Ref. 8 
which can be used for the practical stress analysis of 
frames. The curves apply only to radial loading. but 
Goodey’s method for deducing the corresponding values 
for tangential and moment loading (see Section 15) is 
described in detail. 

The curves give skin shear only. To determine the 
bending moment, et cetera, in the frame, two redun- 
dancies must be evaluated, and this necessitates strain 
energy analysis using the given shear distribution. As 
explained under Ref. 4, the value of the report would 
have been enhanced if these redundancies had been 
tabulated. The curves also require correction for 
eccentricity. 

Curves showing the variation of the skin shear with 
I and ¢ are also presented. These are discussed in 
Section 19. 


Frames of Constant Section with 
Internal Reinforcing 
Figure 9 shows a typical frame of constant section 
with a propped floor beam, loaded with a distributed 
floor load 2W. There will be inter-reactions between 
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the frame and the beam, and these are shown in Figs. 10 
and 11. 

Using a set of curves such as the Royal Aeronautical 
Society Stressed Skin Data Sheets, 03.06 series, for a 
rigid frame or Ref. 7 for a flexible frame, the stress 
analysis is comparatively simple, being essentially the 
determination of three redundancies, M,, P and X. The 


Figure 10. Frame of 
Fig. 9 showing forces 
applied to frame by 
floor beam and props. 


problem can be solved by conventional strain energy 
methods, and one such solution is given in Appendix II. 
It contains nothing new, but is includeed for complete- 
ness and easy reference. A solution of a similar prob- 
lem for a rigid frame (i.e. using skin shear from basic 
engineering theory) was published in Ref. 10. 

12. Cut Frames and Frames with Varying 

Section. General Remarks 


Fuselage frames with large cut-outs are of common 
occurrence in aircraft, e.g. at hatches, in doorways and 
in bomb bays, but hardly any published work touches 
on the problem of their stress analysis. Ref. 12, for 
instance, gives the effect on the shear in the fuselage of 
a torsional loading remote from a long cut-out. But it 
does not examine the part played by the frames them- 
selves in the redistribution of this shear, nor does it 
deal with a load applied at a frame with a cut-out, which 
is the problem under discussion here. 

Similarly, little work has been published on frames 
of widely varying section, although Lewis") has made 
some attempt to assess the significance of changes in I. 

As a first approach, some attempt might be made to 
modify published solutions for complete frames and to 
apply them to frames with cut-outs. Most of these 
solutions, however, use trigonometric series, which are 
relatively easy to handle between limits such as 0° and 


W-Pcos W-P 
OW FiGuRE 11. Forces 
} applied to floor beam 
by frame and _ props. 
(Horizontal compo- 
0 nents omitted.) 
Pcos < 


‘> or 0 and =, but are much more intractable over 
intermediate ranges such as 0° and 130°. 

Some other artifice must therefore be adopted. 
The author has developed a simple ab initio method of 
frame analysis which is particularly suitable for frames 
with large cut-outs. It can be applied with sweeping 
idealisations to give a quick solution for design. It can 
be extended to include greater detail, although it soon 
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STRESSING OF 


becomes unwieldy. 
cai results. 

The method is described in Part II. No attempt is 
m.de to present the results in the form of comprehen- 


It yields some interesting theoreti- 


CIRCULAR FRAMES 


sive nests of curves, as the number of curves required 
would be prohibitive. A rapid method of calculation 
for a given case seems the best that can be hoped for at 
present. 


PART II 
ANALYSIS OF FUSELAGE FRAMES BY THE STRIP METHOD 


13. Introduction 


The method consists essentially of dividing the 
periphery of the skin into a number of equal strips, 
denoting the average shear flow on each of these strips 
by ¥, ... and minimising strain energy by 
Lagrange’s method. 

As the shear in an actual fuselage is usually applied 
to the frame in discrete steps, at first sight the strip 
method would appear to give a closer approximation 
to reality than analyses which represent the shear flow 
by continuous functions. To keep the number of 
variables down to reasonable limits, however. it is 
generally necessary to make the number of steps in the 
calculation less than the actual number of steps that 
occur in practice, so this advantage is unfortunately 
lost. 


A detailed solution for a frame with a cut-out. 
loaded radially, is set out in Section 14, and the method 
of solution for tangential and moment loadings is 
indicated. Numerical examples are worked in Appen- 
dices HI and IV. 


14. Analysis of a Frame with a Cut-out 
Loaded Radially at a Point Dia- 
metrically Opposite to the Centre 
of the Cut-out 


Figure 13 shows a frame with a cut-out. loaded with 
a radial load P at a point diametrically opposite to the 
ceiitre of the cut-out. Let the half frame AB be divided 
into p sections of equal length /. Let the average values 
of the shear flow in the skin over the n'" strip be x,,. and 
the angle to the centre of this strip be z,: x, and z, refer 
to the strip nearest to the cut-out (see Fig. 13). 


Peos oc Pecos 
b 
Figure 12. Effect of 
P on bending moment 
floor beam (See 


Appendix II). 


of 
ra b-2a a | 
Pcose< 


Prose 


(he bending moment at the centre of the m'" strip 
due to shear on the n' strip (n } m) will be a,x, where 
a, is the product of the length of the skin in one section 
and the geometrical offset from the centre of the m' 
Section of the frame to the line of action of the shear 
force on the n'" section. 

let the total bending moment at the centre of the 


m'" section be M,,. Then 


n=m 


M,, ay, nXn e e (4) 
n=1 
eM, _ _ 
also =An»(n =0 (n>™m) (5) 
The strain energy of the half frame in bending will be 


where /,, is the average value of J for the m' section. 
Assume now that the skin shear distribution per- 
sists unchanged for a length L along the fuselage (see 


Ficure 13. Cut frame 
(varying section) divi- 
ded into p segments 
each loaded with a 
discrete increment of 
shear (Section 14). 


(/) of this section), and then reverts to basic sinusoidal 
form. The strain energy in shear of the skin in half 
the fuselage over the length L will be 
21G 
The total strain energy of one half of the frame and 
fuselage skin is thus given by 


m =p M,? 
U=K, 3 x,7+K. —~— +k, (6) 
n=1 m=1 
where K,, K,, K, are terms independent of x,, ... 


K, being LI/(21G), K., being //(2E) and K, representing 
the strain energy in shear of the remainder of the half 
fuselage beyond the length L. (This ignores the energy 
of the stringers due to end load. This is considered in 
(g) of this section.) 

Another relation follows from consideration of the 
static equilibrium in a direction parallel to the line of 
application of P, 


i.e. 
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From (6) and (7), minimising by Lagrange’s Method of 
Undetermined Multipliers (see Appendix 6):— 


M,, 
mn 


m=p 
je. Asin 2,=2K,x,+2K, “Tan, 
From consideration of this result the following points 
arise : — 

(a) If J,, is varying round the frame, it remains 
inside the summation sign, but if it is constant it can be 
taken outside. Similarly if ¢ is varying round the peri- 
meter, it must be removed from K, and placed inside 
the summation sign for the shear term. 

(b) When /,,—> 0% (rigid frame), from (9), 


Asin @,=2K,x,+2K, (mtn) (8) 


(9) 


i.e. the shear distribution is sinusoidal. Thus for a 
rigid cut frame loaded as in Fig. 13, the shear at any 


1000 lb 


pl 
FiGuRE 14. Typical 
partial bulkhead 
loaded radially at top 
centre line. 


point z may be written as K sin z where K is given by 


2r | K sin ad2=P, i.e K r(B—4sin 26)” 


(c) For a complete rigid frame, the shear distribu- 
tion will also be sinusoidal. since the redundancies M,, 
and N of Section 9 give terms in the bending part of 
equation (8), and these vanish when /,,—> ©. More- 
over, putting == in the expression for K, the shear 
flow at point z becomes [P/(=r)]sinz which is in 
complete agreement with simple engineering theory 
(Section 5). 

Values of the shear flow for complete or partial 
rigid frames loaded with tangential or moment loadings 
can similarly be obtained, and are also found to agree 
with simple engineering theory (see Section 15(e)). 

(d) When the frame is rigid over a considerable 
portion of its perimeter, this portion being symmetri- 
cally disposed about A (Fig. 13), and extending from 
n=g~1 to n=p (say), then from (9), 


Asin z,=2K,x,+2K, —— 


npg, 
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i.e. over the rigid portion of the frame. shear flow is 
sinusoidal. This result is of great practical value, as 
the shear flow from n= g+1 to n=p can now be repre- 
sented by K sin z where K is the only unknown. This 
simplifies the solution (but see Appendix VII). 

This method is thus particularly suitable for the 
solution of partial bulkheads of the type shown in 
Fig. 14. A numerical example of such a solution is 
worked in Appendix III. In the example given there 
are seven equations to solve, five of which can be solved 
separately from the other two. Thus the solution is 
obtained without excessive work. 

(ec) When /,, is finite. p equations of the form (9) 
are obtained, giving with equation (7) a total of p~1 
equations. There are p+1 unknowns, namely x, to x, 
and Lagrange’s multiplier 4. The solution is thus 
determinate. 

(f) For a more accurate solution the procedure is as 
follows. The frame is solved with a small number of 
strips as in Appendix III. Then a new block of 
equations is written down, increasing the number of 
strips as required. From a “faired-in”™ graph of the 
first solution, approximate solutions of these equations 
may be read off and more accurate values can be 
obtained quite speedily by relaxation or by applying 
Von Seidel’s method (see Appendix VI). 

(g) For very flexible frames, the shear flow obtained 
by this method shows a tendency to oscillate between 
very large positive and negative values at successive 
stations, giving a quite unrealistic solution. The 
explanation is that rapidly fluctuating shear of this kind 
would be associated with high end load in the skin and 
stringers. The strain energy due to this end load has 
been ignored. If it is included, the peaks of the shear 
curve will be planed off. Computational details are 
given in Appendix V, and an indication as to when this 
energy of skin and stringers in end load should be 
included in the calculation is given in Section 19. (The 
validity of ignoring stringer energy at a rigid frame is 
also discussed). 


oO 
° 


Typical “ oscillating “ shear flow for very flexible 
frames (especially cut frames). See Section 14(g). 
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The tendency for the shear to oscillate still persists 
fo’ a very flexible frame, and this is borne out by 
the curves in Ref. 7 for high values of the parameters A 
and A/B (see also Section 19). 

(h) The effect of eccentricity has been allowed for 
auiomatically in this method of analysis by the use of 
a. (see Sections 8 and 20). 

(i) The choice of a length LZ in which the shear 
reverts to basic presents some difficulty. If there is a 
fairly massive frame only one or two fields away from 
the loaded frame, L should be taken as the distance to 
this frame. If succeeding frames are flexible, L should 
not exceed three times the average frame spacing (i.e. 
effectively allowing the reversion to take place in steps 
over six fields). Unless the bending and shear terms in 
equation (9) happen to come fairly equal in magnitude. 
however, the final result will not be very sensitive to the 
assumed value of L (see Appendix III] for numerical 
examples). 

Alternatively, instead of choosing an arbitrary value 
for L. shear flows in succeeding fields could be denoted 
by Wy Z,. Z,, and the stram 
energy of the skin in shear and the frames in bending 
in these fields could be written down and added to the 
original energy. The work will be greatly increased, as 
the original equations will now contain terms in x and y, 
and there will be an additional set of equations with 
terms in vy and z. The final solution will be more 
accurate. but the method has lost its chief virtue. 
namely its simplicity. 

(j) It is now possible to define a “rigid” frame and 
a “flexible” frame. A rigid frame is a frame whose 


FiGuRE 16. Symmetri- 
cally placed radial. 
tangential and 
moment loadings. 


Struin energy in bending is small compared with the 
Strain energy in shear of the skin in adjacent bays. 
Conversely, a flexible frame is one in which this is not 
the case. This definition agrees with Ref. 1, page 15. 
It is given a quantitative character in Section 19 of the 
present paper. 

(A) It must be pointed out that when partial or 
complete bulkheads are examined by the methods of 
this paper, no information is given on the distribution 
of the shear stresses in the bulkhead web. This could 
be determined by ordinary engineering methods, or by 
more refined analysis, but in either case the skin shear 
distribution given by this paper could be used as a 
Starting point. since it is unlikely that the web energy 
would be high enough materially to affect the skin shear 
distribution, 

(/) As in all strain energy calculations, the main 
block of equations will form an array symmetrical 
about the leading diagonal. Special methods exist for 
the solution of such equations. For the convenience of 
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readers, some of these methods are outlined in 
Appendix VI. 

(m) It might be argued that a substitution method 
could be used in preference to Lagrange’s undetermined 
multipliers. This is true, but the elegance and sym- 
metry are lost, and the theoretical results emerge less 
clearly. 


15. Cut Frames with Other Forms of 
Loading and with Internal 
Restraints 


(a) Two symmetrically placed radial loads 

If the radial loads P are applied at A, A’ (Fig. 16). 
the bending moment M,, at all sections where z,, < 
will contain a term in P. Equation (4) will be modified 
to:— 


n=m 


M,.= Qnix,—Prsin (¢— 2,,) 


n=1 


(10) 


n=m 
In Lo, & 
where r’ is radius of centre line of frame at section /. 
The equilibrium equation now becomes 


n=p 

Equation (9) is unchanged. 

Reference to Appendix III shows that the bulk of 
the computational: work from the previous (radial) load- 
ing stands unchanged. The constants on the right hand 
side of the block of equations must be computed afresh 
and equation (11), giving the value of A, is different. 
These are the only alterations. The solution to this 
loading can thus be obtained expeditiously from the 
solution to the preceding loading. 


(b) Two symmetrically placed tangential loads (see 
Fig. 16) 
At sections where z,, > @ the expression for M,, will 
contain a term involving T. Equation (4) will be 
modified to 


an >>, Mn= cos (¢ — z,,)] 
n=m (12) 
Gam 
n=) 
The equilibrium equation becomes 
—sing= x, sin z, 


n=1 


Equation (9) is unchanged. Also, as in (a), the solu- 
tion follows readily from the numerical solution of the 
case dealt with in Section 14. 


(c) Two symmetrically placed moment loads (see 
Fig. 16) 
Equation (4) becomes 
n=m 
9, M,, AnnXn- M, 
n=1 
(14) 
n m 
Zm { M,, An nXn 
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The equilibrium equation (7) becomes :— 


> Kise 


n=1 


(15) 


Equation (9) is unchanged, and the remarks in (a) and 
(b) still apply. 
(d) Asymmetric Loadings 

With an asymmetric loading of any type. the method 
of procedure will be as in (a), (b) and (c), but twice the 
number of unknowns will be required, and integration 
must be carried out over the whole perimeter (instead 
of half as previously). For the loading shown in 
Fig. 17 the equilibrium equations will be :— 


n= 2p P 
X, Sin 2, 
n= l 
n= 25 
XCos' a, =0 (16) 
n=1 
Pb 
and 
n=1 lr 
where b, is offset of P from centre of frame. Three 


Lagrange multipliers will be required and equation (9) 
will become : — 


2K = 


man Ey, 


A, sin 2, +A, Cos z, +A, = (17) 
There will be 2p +2 equations to solve, and the method 
is beginning to lose its simplicity. It will be better 
therefore to break the loading down into a combination 
of a symmetrical load and a moment load of the type 
given in (e) if this is possible without too great a loss 
of accuracy. 
(e) A moment load or a tangential load applied at the 
point diametrically opposite the centre of the 
cut-out (Fig. 18) 

The shear on the two sides of the frame is now anti- 
symmetric, so only p variables are required. ™,, is 
unchanged, but there will be two equilibrium equations 
for each loading. 


Moment Loading : — 


n=p M 
0 
2rl 
(18) 
Xcos'a,,=0 
Tangential Loading : — 
n=p 
| 
(18a) 
n=p 
= Xn COS 2, = | 
n=1 2! 


Equation (9) will thus contain two Langrange operators, 
and for both moment and tangential loadings it will 
become : — 


m=p M 
m 


A+ PCOS @,=2K,X,+2K, — ann (19) 


m=n 


The case when the frame is completely filled with a 
shear web (see Fig. 27), is of considerable practical 
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value. The bending terms disappear from equation 
(19), which could then be rewritten 


X,=a+bcos 2, 


where a and + are constants. These are easily deter- 
mined from considerations of overall equilibrium as 
follows. 


Pp 


xp! pti 


FiGURE 17. Asymmet- 


\ ric loading showing 
1 reacting shear (2p un- 
knowns. See Section 
B | 15(d)). 
x, \ 


Moment Loading : — 


(a+bcos z)rdz 


B 

and | (acos z+ b cos*?z)rdz=0, 
0 

giving af+bsin B= 


F [- 
asin 8 + —0, 


which yields 


M, 8+4 sin 26 
M, sin l 
b= 
2r° 
sin’ — 


Tangential Loading : — 
B 
| (a+bcos 2)r°?dz=Tr 
v 
B 


T 
| (acos a+bcos? z)dz= — , 
2r 


0 


which similarly yields 


[ 
T (sins 28) 
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When B=7, i.e. a complete bulkhead, these relations 
yield: — 


M ment Loading. Shear flow is 


+ Cos 2). 


Tengential Loading. Shear flow is 


These results agree with equation (3). 


16. An Alternative Method of Solution for 
Moment and Tangential Loading 


if the distribution of the reacting shear is known 
for a radial load P applied at a point, the distribution 
for a tangential load T or a moment load M, applied 
at the same point may be deduced directly. The method 
is due to Goodey (Ref. 8, Fig. 5), and is also referred 
to by Lewis (Ref. 11, pp. 6-8). The shear flow for an 
applied radial load is shown to be the differential co- 
efficient of the shear flow for an applied tangential load. 
Conversely, the shear flow for a tangential load may 
be deduced from that for a radial load by a single 
integration with adjustment for the arbitrary constant. 

In practice, the shear flow for the radial load is 
plotted against 4, integrated, multiplied by the ratio T/P 
and the axis adjusted for the constant of integration. 
This constant can be estimated by considerations of 
overall equilibrium. More simply, it can be determined 
by inspection if the value of the shear flow at any one 
point on the perimeter can be deduced from general 
considerations. 

An exactly similar procedure can be applied to 
bending moments, axial loads and web shears. (This 
was also pointed out by Goodey™’.) The exercise of 


Figure 18. Moment 

or tangential loading 

\ at B showing anti- 

symmetric shear (p 

| unknowns only. See 
Section 15(e)). 


a little ingenuity generally suffices to determine the 
arbitrary constants without much calculation. 

For example, quoting from Wise (Ref. 9, equation 
(16) ), the axial load in the frame due to a radial load P 
is given by 

N= 7 008 4+ (= — 4) sin |. 
Multiplying by T/P and integrating with respect to 6, 
gives for a tangential load T, 
| 


T 


where A is an arbitrary constant. But at 6=0°, the 
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_ 
point of application of the load, it is obvious from Fig. 6 


of Ref. 9 that by symmetry N=—4T (or at 6=7, 
-+3T). This gives A=0, and hence 


N= E sin (= |. 


agreeing with equation (18) of the same reference. 

For a moment load M, the shear flow for a tan- 
gential load is differentiated twice with respect to 6, 
and the result is added to the original shear flow for the 
tangential load. On multiplying by the ratio M,/(rT), 
the shear flow for a moment load is obtained. A similar 
procedure can be applied to bending moments, axial 
loads, and web shears in the frame. 

Once again, quoting from Wise (Ref. 9, equation 
(18)) the bending moment M due to a tangential load T 
is given by 


M= 4)(1 — cos @)- | 


Tr 


hence 


[ — 4) cos 6— sin 
and for the moment load M,, 


M= [ 2sin# |, 
agreeing with equation (13) of the same reference. 

The method appears unsuitable for use with shears 
obiained graphically by consideration of a few points 
only, for differential coefficients calculated as slopes of 
such curves are notoriously unreliable. (If however 
the second differential coefficient is small compared 
with the original term, this objection disappears.) 

In practical calculations the first differential co- 
efficient of the shear from the radial load would be used 
in place of the second differential coefficient of the shear 
from the tangential load. 


17. Restraint from the Wing 

The effect of restraint from the wing can readily be 
estimated by standard engineering methods, applied to 
the results of an analysis based on any of the methods 
described in this paper. One simple procedure is as 
follows. 

If the restraint occurs at points defined by 6,, 6,, use 
the selected method to determine the bending moment 
in the frame for unit load applied at 6, in the direction 
of the restraint, and similarly for unit load applied at 
6,. Determine also the bending moment due to the 
given system of external loading. 

Double integration then yields the deflection at 6, 
and also at 4, for. (a) an assumed loading of X, at 4,. 
(b) an assumed loading of X, at 6., (c) the given external 
loading. Equating the total deflection at 6, and 64, 


separately to zero gives two equations for X, and X,. 

The effect of a restraint on the “leg” of a partial 
bulkhead is to relieve the bending stresses at the top of 
the leg (P and P’ in Fig. 14). Normally only a portion 
of this calculated relief could be safely applied as the 
restraint is unlikely to be absolute. 
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PART III i 
A DISCUSSION OF THE UNDERLYING PARAMETERS WHICH G 


DETERMINE THE 


18. Introduction 

Analysis by the “strip method” of Part II throws 
considerable light on the nature of the underlying para- 
meters which determine the stress distribution in the 
skin and the frame. In the succeeding sections, these 
parameters are discussed and the predictions of the 
“strip” theory are compared with the findings of other 
investigators. 

In general. agreement is good, but some previously 
published results are shown to be in error, as important 
factors have been ignored or results have been plotted 
against incorrect parameters. 


19. Variation with J, r, and 

It was shown in Section 14 that the degree of dis- 
tortion of the skin shear distribution pattern depended 
on the ratio of the strain energy of the frame in bending 
to the strain energy of the skin in shear. From Section 
14 these are respectively 


Dimensionally 
M,, ocl x geometrical offset x 


Hence. dimensionally, strain energy of frame in bending 


and strain energy of skin in shear 

° 

The ratio of these yields a non-dimensional parameter 
Gtr'/(EIL). For small values of this parameter, the 
skin shear will be “basic.” For high values the skin 
shear will be distorted. 

This parameter is identical with the parameter A/B 
of Ref. 7. Moreover the curves in Ref. 7 show that 
for A/B=0 the skin shear is basic, and for increasing 
values the skin shear pattern becomes increasingly 
distorted. There is thus good agreement between Ref. 
7 and the theory of Part IT on this point. 

There is another such parameter common to Ref. 7 
and the present paper; this is representative of the 
importance of the skin and stringers in carrying end 
load. From Appendix IV, the strain energy of the 
skin and stringers over a length L of the fuselage is 

Elt’ Ert’ 

Likewise the strain energy of the frame in bending 
is proportional to x,*r°/(E/). The ratio of these gives 
the non-dimensional parameter L*//(r't’), which is the 
reciprocal of the parameter A in Ref. 7. Small values 
of A thus indicate high stringer energy, with consequent 
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relief to shear distortion (see Section 14(g) and Appendix 
IV). Large values of A represent low stringer energy 
with little relief to shear distortion. 

So far it has been implicitly assumed that when the 
frame bending energy is negligible (rigid frame) the T 
stress distribution is controlled by the skin shear energy 
and that provided that the stringers are distributed 
evenly round the frame, their energy has no effect on 
the shear distribution. This is in line with engineering © 
experience, and with the results of Ref. 7, but it is more | 
satisfactory to demonstrate mathematically that this is 
indeed the case. Such a demonstration is given in 
Appendix VII. (It is shown to be only approximately 
true for an incomplete frame but strictly true for a 
complete frame.) 

The following is thus obtained. At small values of | y} 
A/B, there is little shear distortion, and the shear | jj, 
pattern is relatively insensitive to variations in A (pro- | 
vided that the end load carrying material is distributed | 
evenly round the frame). At high values of A/B, shear | jp; 
distortion is considerable, and increases with increasing |) ., 
values of A. This is fully in agreement with the curves 
plotted in Ref. 7. 

With the aid of Ref. 7. the following, therefore, can 7 gy, 
be deduced : — Hi 

(i) When A/B is small (less than 1-0, say), the FX ; 
frame can be treated as rigid. 

(ii) When A4/B is larger (not greater than 10, say), 

the frame is slightly flexible. Stringer energy 
(or more strictly, energy of skin and stringers | 
in end load) can be neglected in strain energy | 
calculations. Analysis based on engineering | 
theory may be used if desired, the maximum | ; 
skin shear being then underestimated (by 
about 10 per cent. when A/B is 10) and the 
maximum bending moment being over-" 
estimated (by about 15 per cent. when A/By ;, 

is 10). ent 

(iii) When A/B is large (of the order of 200 up- § of | 

wards, say), the frame is truly flexible. | of . 
Stringer energy is important (unless A is large, | agr 
of the order of 10°, say). Analysis based on 
engineering theory will lead to substantial | ; 
errors. 

A quantitative character has thus been given to the | 
definition of Section i4(j). It should be noted that these > 
criteria apply primarily to complete frames and_ that 
frames with large cut-outs will tend to develop 
“flexibility at smaller values of 4/B. 

The non-dimensional parameters A and A/B appear | 
to be of so fundamental a nature that they should figure | 
(directly or indirectly) in all sound frame analyses. It} 
is interesting to look for them in other published work. | 

Starting with Ref. 8, p. 848; in the derivation of 


strir 


equation (D7) Goodey speaks of two non-dimensional F ied 
parameters “c” and “e” which represent 
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“fundamental characteristics of the structure.” Using 
Goodey’s own notation, these are respectively 
R°K.t 


3 l R = 
c= 3m (1 + where m= 1, 
K = rt being in the notation of the present paper) 


e= 3m? (2K,(1+%)—»). 
These yield 


SRY 6R°t’ 
ES = A, 
. 
and e = since = 21 +»), 
3R? VE 


where » is Poisson’s ratio and r a constant representing 
a comparison between frame stiffnesses, r being unity 
Examina- 
tion of equation (D12) in Ref. 8 shows that the shear 
flow is dependent on e/c and also on c, but as there are 


no curves showing the variation of the shear flow with 
these constants, it is not possible to deduce their effect 
so readily as in Ref. 7. 


Equations (D7) and (D12) of Ref. 8 refer to a short 
length of fuselage. For a long fuselage, Goodey intro- 
duces a parameter A which is dependent on his c and e. 
His final expression for the shear flow is dependent on 
\ and hence also on c and e, but once again the exact 
method of dependence is not easy to determine without 
complicated algebraic manipulation or parametric 
plotting. 

The parameters A and B also appear in Ref. 4. Hoff 
points out their importance, but with rather different 
Also, Hoff’s A 


the A and B of Ref. 7 and the present paper. 
Reference | also uses a non-dimensional parameter 
d which is defined as kKR*/(E/) (notation of Ref. 1). k 


is later defined as RrG/L, giving d=GtR'/(EIL), 


entirely equivalent to the 4/B of the present paper or 
of Ref. 7. Fig. | of Ref. 1 also gives numerical values 
of d, above which a frame is “ flexible.” These values 
agree well with those given in this paper. 

The parameter A does not, however, appear. This 


indicates that the results of Ref. 1 apply only to 
fuselages where the energy of the skin and stringers in 


end load is negligible. The effect of neglecting A is 


shown in Fig. 29, which is a cross plot from the curves 


of Ref. 7. The quantity A, is defined as 


Maximum bending moment (or skin shear, etc.) 
for a given value of A 


ippear A _ = 

figure | Maximum bending moment (or skin shear, etc.) 
s. Ith 

10° may be taken as “strain energy of skin and 
ion 


sional 


Stringers negligible,” i.e. skin and stringers stiff in end 


‘load. The curves show clearly, that unless this energy 


is negligible, large errors can occur by neglecting it. 
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Beskin® used a parameter AL/[(1 9) (in his 
notation A is the total cross-sectional area of the frame). 
Since I=An’?R*, the parameter may be written 
IL/{R*ty*(1—»)] or approximately 
Ignoring the eccentricity, this gives a numerical multiple 
of the reciprocal of the parameter A /B. 

Beskin uses another parameter Rh /(L*t) (again in 
his own notation, h being equivalent to ¢’ in this paper). 
This is a numerical multiple of the B of Ref. 7. 
Although Beskin states that this second parameter has 
little effect on the stress distribution, since he is 
primarily investigating the effects of eccentricity, his 
statement must be construed as meaning that, in the 
range he investigated, the effects of eccentricity on the 
stress distribution were not sensitive to variations in B. 
(Further discussion of the effect of eccentricity will be 
found in Section 20.) 


It is gratifying, therefore, that so simple an analysis 
as the present paper offers should agree with other more 
elaborate investigations, including one so basic as 
Goodey’s, as to the nature of the fundamental constants 
that determine the skin shear. 

Lewis‘'') produces curves for skin shear based on 
Goodey’s analysis in Ref. 8. These include some 
interesting curves showing the variation of skin shear 
with 7, L, t and r (in the notation of the present paper) 
and also showing the effect of keeping the J of the 
loaded frame constant but varying the J of adjacent 
frames. He concludes that within reasonable limits the 
skin shear is not sensitive to quite large variations in 
these quantities. This is in general agreement with the 
present paper but seems to require some qualification. 

Thus, consider variation with J. The sensitivity is 
small if 4/B is less than 1-0 (say), as in this range the 
bending energy of the frame is trivial. It can also be 
small when A /B is large (greater than 1,000, say) as in 
this range stringer energy can become important com- 
pared to bending energy. Maximum sensitivity would 
therefore be expected in intermediate ranges. For 
example, in the neighbourhood of A/B=200, local 
variations in J of the type used in Appendix III could 
conceivably be highly significant. Also, when 4/B is 
large (greater than 1,000, say) they could again be sig- 
nificant if stringer energy and bending energy happened 
to be of about the same magnitude. But the point 
requires examination using actual numerical examples. 

Similar remarks apply to variations in L, ¢ and f. 
The range of greatest sensitivity to L (average frame 
spacing) would probably be A/B greater than 1,000 
with A not greater than 10', i.e. flexible frames with high 
stringer energy. 

Section 6 of the “ Conclusions ” and Fig. 21 (of Ref. 
11) also seem to require modification. Fig. 21 presents 
a plot of / against r (in the notation of the present paper) 
and states that engineering beam theory may be used 
if J falls above the line. It has already been shown, 


however, that the fundamental parameter that deter- 
mines this is Gtr'/(EIL). To plot / against r is there- 
for misleading. J should be plotted against Grr'/(EIL) 
or against c/e (Goodey’s constants). Selecting numeri- 
cal examples at random gives the values in Table I. 


j 
| 
ix: 
| 
i 
x 
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From Fig. 21 of Ref. 11, all the frames would have 
been treated as rigid, but Ref. 7 and Table I indicate 
that this is incorrect except for the first. 

In the discussion of the present section, L has been 
used in two senses—to represent frame spacing, when 
referring to the work of other investigators, and to 
represent “the length along the fuselage in which the 
shear pattern reverts to basic” when referring to the 
results previously deduced in the present paper. There 
is an apparent discrepancy, but this is not serious as 
the argument of the present section is based on dimen- 
sions, and since the distorted shear reverts to basic in 
a small integral number of bays, the two senses are 
dimensionally similar. For calculating numerical values 
of A and 4/B, however, L should be taken as meaning 
“frame spacing.” 

One other point arises in connection with the effect 
of stringers on shear distribution. Hoff has examined 
the effect of the applied load being reacted partly by 


TABLE I 

Bending 

moment 

Shear cone reduction 
A_Gtr factor from factor from 

- BIL 1L3 Ref. 7 Ref. 7 
20 20 0:036 0:07 8 14 70 (Fig. 6) 0°80 (Fig. 3) 
50 50 0:048 0:10 10 1.250 2-0 (Fig. 21) 0°55 (Fig. 18) 
70 70 0:064 0:12 12 720 3,400 2-8 (Fig. 21) 0°45 (Fig. 18) 


(Values of A and A/B are approximate only. Figure references are to Ref. 7.) 


radial and tangential bending of the stringers. Figs. 13 
to 15 of Ref. 4 show the magnitude of this effect. It is 
small if A is less than 10° . (Note that A has values 
six times greater numerically than in the present paper.) 

Goodey? also considered the effect of bending and 
torsion of the stringers, although he ultimately ignored 
it in the interests of simplicity. 

The present analysis could easily be extended to 
include this effect. It is simpler, however, to ignore it 
in actual calculation and to correct the peak values 
afterwards, using Hoff’s curves. It will rarely be 
significant. 


20. Variation with Eccentricity (’) 

Most of the published analyses on fuselage frames 
treat the skin line as coincident with the centre line of 
the frame. Beskin™, however, has investigated the 
effect of eccentricity, i.e. the offset of the skin line from 
the frame centre line. 

Beskin calls the skin radius R and the radius of the 
frame centre line (1 —) R, where » is the eccentricity. 
He takes A as the area of cross section of the frame. 
and An?R* as the moment of inertia of the cross section 
of the frame. As already explained in Section 19, he 
finds that the effect of eccentricity depends on two para- 
meters, sensitivity to one of these being small. 

A cross plot* of Beskin’s curves is given in Fig. 30. 
This shows the effect of eccentricity on maximum frame 
stresses, plotted against the parameter A/B (notation 


*Minor irregularities have not been “ faired in.” They would 
probably disappear if higher harmonics were considered in the 
original solution. 
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of the present paper), values of A/B having been 
obtained from Beskin’s parameter AL/[(1 1) tR*] 
(Beskin’s notation) by using the conversion : — 


A/B=* Reciprocal of Beskin’s parameter ” 
x 


X,, is defined as 


Maximum bending moment (or web shear, etc.) 


‘ for a given value of 7 
' Maximum bending moment (or web shear, etc.) 
for 


Surprisingly, the cross plot shows that for a rigid 
frame loaded radially the maximum bending moment is 
independent of the eccentricity. This case can readily 
be analysed by routine strain energy methods, and such 
an analysis is given (for convenience) in the second part 
of Appendix I of the present paper. It will be seen 
that the maximum bending moment comes out as 
3Pr (1 — )/(4=) or approximately 0-24 (1 — 1) Pr, whereas 
Beskin’s curves give 0:24 PR, both r and R being skin 
radius. The author is not able, at present, to offer an 
explanation of the discrepancy. 

Beskin does not show the effect of eccentricity on 
skin shear distribution. Hoff’, however, presents 
curves showing the effect of eccentricity on moment 
reduction factor and skin shear concentration factors, 
plotted against his parameter A. Thus, from Figs. 5 
and 6 of Ref. 4, when A=10', varying » from zero to 
0-075 makes little difference to the shear concentration 
factor, but increases the moment reduction factor by 
about 10 per cent. For A=10*, varying 4 from zero to 
0-075 decreases the shear concentration factor to about 
half, and increases the moment reduction factor by 
about 25 per cent. 


To interpret these results, it is necessary to decide | 


whether the skin radius is kept constant and the frame 
radius diminished, or vice versa. 
radius constant, the results mean that increasing 
eccentricity decreases the skin shear distortion and 
increases the frame bending moment. 


Predictions of the effect of eccentricity from the 


theory developed in Part II are rather unsatisfactory. | 


Thus, keeping the skin radius constant and allowing the 
frame radius to decrease, but also allowing the frame 
moment of inertia to increase as in Beskin’s article, has 
the following effects: (a), terms in equation (9) from 
skin shear will remain unaltered, (4), terms in equation 
(9) from bending energy will be decreased on two 
counts—the bending moment is reduced as already 
explained and the moment of inertia is increased. The 
shear terms therefore become relatively more important 
with increasing eccentricity and hence there is less shear 
distortion, with a consequent increase in the bending 
moment. Whether the final effect is a net increase or a 
net decrease seems difficult to predict. But the decrease 
in shear distortion at any rate is in agreement with 
Hoff’s curves. Moreover, the effect should be more 
marked at high values of A than at low values, for low 
values of A indicate high stringer energy, in which case 
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small changes in frame bending energy are less signifi- 
cant. This, too, is in keeping with Hoff’s curves. 


21. Variation with Number of Fields 


Hoff has investigated the effect of varying the 
number of fields between the loaded frame and the 
support. He finds maximum distortion of the skin 
shear pattern with one field. and decreasing distortion 
with increasing number of fields. The maximum bend- 
ing Moment increases with the number of fields. The 
curves showing these effects tend to flatten off after 
three or four fields. 

Other investigators have pointed out that increasing 
the number of fields above three has little effect on the 
maximum shear and bending moment (see Ref. 7. Fig. 6 
and Ref. I1, p. 5). 

Goodey’s analysis“ includes the effect of fields on 
the side of the loaded frame remote from the support 
(i.e. to the right of P in Fig. | of the present paper). 
Compatibility demands that there shall be self- 
equilibrating stress systems in this field, but their mag- 
nitude appears to be small in general. 

The present theory explains the significance of the 
number of fields as follows. Variation of the number 
of fields gives variation in L, the length in which the 
loaded shear reverts to basic. When 4/B is small, skin 
shear energy is paramount, and the distribution is not 
sensitive to L (see, however, Appendix VII). When 
A /B is in an intermediate range, where shear and bend- 
ing energies are comparable, the shear distribution will 
be sensitive to L and therefore to the number of fields. 
When A/B is large, stringer energy is becoming impor- 
tant compared with bending energy and once again the 
final shear distribution will depend on L and therefore 
on the number of fields. 

This does not explain the generally observed 
phenomenon that sensitivity is small once the number 
of fields exceeds three or four. It is suspected that this 
is more of a coincidence, based on the range inside 
which practical aircraft designers work, than a funda- 
mental theoretical necessity. It is brought about by 
the fact that the bending energy produced in succeeding 
frames by change of the shear flow at these frames is less 
than the additional shear energy in the skin if the dis- 
torted shear had continued unchanged. 

The preceding considerations suggest one point of 
practical importance. If a short length of fuselage is 
used on test to represent a longer fuselage, provided 
that there are at least six fields and the loaded frame is 
two fields from the open end, the maximum measured 
shear flows and bending stresses should adequately 
represent the maximum values actually occurring on 
the longer fuselage. 


22. Extensional and Shearing Rigidities of 
the Frames 


The present analysis has so far taken no account of 
the extensional and shearing rigidities of the frames. 
These will be important when the strain energy of the 
frame in axial load or of the web in shear is large 
compared with the bending energy of the frame. 
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Dimensionally, strain energy from axial load 
depends on (/x,)*[//(ES)], i.e. r°x,?/(ES), and due to 
web shear (/x,,)°[//(GS’)], i.e. r°x,.7/(GS’). From Section 
19 strain energy in bending depends on r°x,”/(EI). The 
relevant ratios are thus //(Sr*) for axial load and 
1/{((G/E)S’r?] for web shear. Hoff? calls the first ratio 
C, and the second is equivalent to his CD, where he 
defines D as S/(GS’/E). 

The present theory thus predicts that the axial and 
shearing rigidities will be important as C and CD 
increase. In Ref. 4, Figs. 9 and 10, Hoff shows the 
effect of variation of C (keeping D constant) and D 
(keeping C constant). The shear concentration surpris- 
ingly goes down with increasing C, which appears to be 
contrary to the present theory. The apparent dis- 
crepancy arises perhaps from the manner of variation 
of the parameters, which Hoff does not specify, except 
that he keeps A and B constant. If / and r are kept 
constant while § varies, increase of C increases the total 
strain energy of the frame and hence the shear concen- 
tration factor. If, however, J, r and S all vary it might 
be possible for an increase of C to be accompanied by 
a decrease in the total strain energy of the frame com- 
pared with the strain energy of the skin in shear, and 
this would give a decrease in the shear concentration. 

In general, the effect of including shearing and 
extensional rigidities of the frame is small, but where 
there is reason to believe these effects are significant. 
terms giving the strain energy of the frame in axial load 
and the web in shear in terms of x, to x, can easily be 
included in equations (6) to (9) of Section 14. 
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APPENDIX I 


STRAIN ENERGY ANALYSIS OF A COMPLETE CIRCULAR 

FRAME OF CONSTANT CROSS SECTION, LOADED RADIALLY 

AT ONE POINT, THE REACTING SHEAR BEING DISTRIBUTED 
SINUSOIDALLY 


The radial load acts at A (Fig. 8), the frame is 
assumed cut at B and redundancies N, V, M, are intro- 
duced. By symmetry V=0. At the point defined by 6, 
the shear flow is [P/(=r)]sin@. Hence at the point de- 
fined by z, the bending moment in the frame is 


M=M,-Nr(1+cos z)+ | (2 sin [1 —cos (4 — z)] dé 


=M,-—WNr(1+cos 2)+ Pr +.cos 2—4 (x — x) sin 2} 


where positive M gives tension in the inner flange. 
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cM giving 
Then and aN r(1+cos 2). 
From accepted strain energy theory, | Mdz=0 and | M cos zx dz=0, 
0 0 
| d2=0 and | which yields 
0 


Combining these four conditions gives 


) Maz=0 and M cos xdz=0 
0 
and hence 
M, —Nrz id 
Nrz  3Pr 
. _ 
1.€. = iz M,= an 
(N_ positive, tension. M, positive, tension on inner 
flange) 


Finally M= =4cos 2—(=—2)sin a} 


which agrees with equation (16) of Ref. 9. 


NOTE (i) Where the shear flow is not capable of 
analytic expression, numerical values of M and M cos z 
must be tabulated round the frame and the integrals 


| Mdz and M cos zdz 


must be evaluated graphically. The work is consider- 
able. as to get numerical values of M and M cos z at any 
one station requires graphical integration also. Terms 
involving M, and N can, however, be evaluated analyti- 
cally throughout the calculation. 

(ii) In the preceding, the skin line has implicitly been 
assumed coincident with the centre line of the frame. 
If this is not the case (i.e. if there is “ eccentricity ” (see 
Section 20), taking the skin radius as r and the frame 
centre line radius as r’, and applying the redundancies 
N, M at the frame centre line. gives 


+r’ cos z)+ 


| (sin r’ cos (4 — rd4 


= M,, — Nr’ (1 +cos z)- 


Proceeding as before : — 


6M _ ] oM r (1 +cos z) 
OM, aN SZ), 
[ ..aM 
| Moy 0 and 0 


+P(r- =6 


=0 


Then, putting 


and m= (+4 cos z+ (= z)sin 2). 


The maximum value of M (at z=0) is 


3Pr 
M= — »). 


This result is discussed in Section 20. 


APPENDIX II 


SOLUTION OF A COMPLETE CIRCULAR FRAME OF CONSTANT 
CROSS SECTION WITH A PROPPED FLOOR BEAM, LOADS 
BEING APPLIED THROUGH THE FLOOR BEAM 


(See Figs. 9, 10 and 11 and Section 11) 


The solution consists essentially of the determina- 
tion of three redundancies, M,, P and X by strain 
energy analysis, using curves such as those in the Royal 
Aeronautical Society Stressed Skin Data Sheets 06.03 
series for a rigid frame, or N.A.C.A. T.N.1310 (see Ref. 
7) for a flexible frame. 

The bending moment at any point in the frame is 
given by 


(W —Pcos z)+m.X +m,M, +m,P., 
where #,, m,, m,, m, are obtained from the curves. 
Then 


(m,W +m.,X +m,M,,+m.PY 

0 
b 


2) 


1s 


0 


Where Usrane ANd denote respectively the strain 
energy of the frame and the beam, s is the distance 
measured round the frame from some fixed point, M, 
is the free bending moment in the floor beam from the 
distributed load (i.e. M, is (W/b)(bx-.x*)), m, is a 
constant evaluated in Note (a) of this Appendix, and 
mM, COS 2. 

The foregoing ignores the effect of the horizontal 
components Pcos z and X on the bending moment. In 
general this is small; it can be included if desired. 
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forced with a propped floor beam (See Fig. 9 and Appendix II). 


Then if U is total strain energy, U = Usrame + Utwams 
and 


Qtr 27r 


0 
0 0 
+M, | m,m,ds +P - 
0 0 
+M, + P{m. — 
0 0 
b b 
| M,dx-P | m,dx - | —0, 
beam 
oar 
| mym.ds+ X [ man.ds + 
0 0 
Mf) m,m,ds 4 m.7ds 
0 0 


b b b 


[ [mcd M, | mats | 0. 


0 0 0 


I, am 


These integrals are evaluated graphically and the 
equations solved for X, M, and P. By a further 
application of the curves already referred to, the frame 
can now be completely solved. Typical skin shear and 
frame bending moment distributions are shown in Figs. 
19 and 20. 


reinforced with a propped floor beam (See Fig. 9 and 
Appendix II). 


NOTES. (a) Integrals containing m 
From Fig. 12, 
xcosa 
M1, = A COS 


from x=0 
from x=a 


to x=a, 
to x=6/2. 
whence 


a 
| x cos zdx +2 | acos zdx—a(b— a)cos z, 


0 


mdx =2 cos? zdx—2 | cos? zdx 
0 0 
= (3b - 4a} cos? z. 
| Mudx (ox x*)dx= 
b a 
Mm, dx | (bx? — dx + 


2 

cos : {a (bx dx 
a 

_ Wacos z 


6b 


(b) m. will consist of four parts arising respectively 
from the radial and tangential components of XY on the 
left side and right side of the frame. In general two of 
these will be negligible. Similarly there will be four 
components in each of and 

(c) Values of m,, m,, m, and m, should be tabu- 
lated from 0° to 180° and also values of their squares 


b+ b*). 


1953 
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and cross products. Integration can then be performed 
graphically. It is advisable to tabulate a considerable 
number of points, as there are numerous peaks on the 
cross product curves, and if some of these are over- 
looked considerable errors can result. 

(d) A solution for a concentrated load at the centre 
of the floor beam is readily obtained from the preceding, 
the only change being in the value of M,. Integrals 
containing M, must be recalculated, and the equations 
solved again, but the great bulk of the work stands 
unchanged. 

(e) A direct solution for an asymmetric floor load 
would be more laborious, as there would be seven 
redundancies. It can, however, be broken up into a 
symmetric plus an anti-symmetric load, each of which 
involves only three redundancies. This reduces the 
computation work. 


APPENDIX III 


NUMERICAL SOLUTION FOR A PARTIAL BULKHEAD 
(See Fig. 14 and Section 14(d)) 


r=70 in., 2=55°, B=135°, L=55 in., K,=LI1/(2tG)= 
0-279 x 10-7, x10~-*. Shear over 
PP’ =x, sn 

Length of strips=16°, i.e. /=19-5 in., t=0-048 in.. 


a. 1=42-2 a, /1=27-0 1=64 1=26 


 ay/L=1S°7  ay/ l= 77 | ay/1l=4-0 
a.,1=161  a/l= ay/l= 5:2 

90 a,/l= 58 

a.. 6°5 


i= 25, i=80,.. 1,=27,. 1,=46 
M,=19°5 (2-6x,) 

M,=19°'5 (6-4x,+ 40x.) 

M,=19-5(15-1x,+ 7:7x.+ 5:2 x;) 
8:4x,+5-8x,) 
M.=19-5 (42:2 x, +. 27-4.x. + 16-1 x, x, +6°5 x;). 


The strain energy of the skin over PP’ is 


55° 


2x 0-048 x 4x 108 J 27a 


= 101 x10-*x,? =2-47 x 10-*x,?. 


The equilibrium equation is 


55° 


70 x, sin? zdz+ 19-5 sin z, = 500, 
a=1 
ie. 0-88x,+ x, sin z,=25-7. 


Minimising by Lagrange’s method (see Appendix VI). 


For x, 0-88’ =4-94 x 10-*x,, 
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FOr x; tO 


Asin 2,=5-58 x 10-*x, x 10-7 3 dun» (a) 
Let Qn = m 19-5 
Dividing equation (a) by 19°5° x 10 
Q,+7:°53x,=1350A sinz, . (b) 


Values of Q, are now calculated in tabular form 
(Table II). 


This gives the following equations, by substitution into 
equation (b):— 
96-2x, +28-4x, +141x, +5-95x,= 1080) 


S1-8x, -39-2x, +17-1x, + 8-74x, +3-86x,= 1260 
28-4x, +17°53x, + 4:95x, +2:27x,= 13450 
I4dx, ~ 8-74x, + 495x, +10°S4x, +1:27x, = 1326 
5-95x, + 386x, + 2-27x, + 1-27x, +8-45x,=1205A 

165K 


0:799x, 0°934x,+ 0:996x,+ 0°982x, +0°891x, + 0°88x, = 25-7. 
Solution of the first five equations yields 


X,= 54-7 x,=20"1 A, x, =89°6 A, 120-0 A, 
tA. 


Substitution of these values in the last equation yields 
A=0°0727, and hence x,= —3-07,. x.=1-46, x,=—5-02, 
x,=6-°73, x,;=7:15, x,=10-0. (All these have been cal- 
culated to slide rule accuracy only.) 

From these values the bending moment is obtained 
by substituting in the expressions for P, to P;. Axial 
load and web shear at z,, are given respectively by 


n=m 


n=m 
19°5 >. x, cos (4, 2,,) and x, sin —2@,,); 


n=1 =1 


with correspondingly adjusted formulae over the 
region PP’. 


LOADING AS FIG. 


SKIN SHEAR FROM 
APPENDIX IIT 
S (L=55") 
/ SKIN SHEAR 
=| FROM RIGID FRAME 
0° 0 


Figure 21. Skin shear flow as given by Appendix III for 


partial bulkhead, compared with skin shear flow from 
engineering theory. 
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A small refinement must be used in the calculation 
of the axial load. In the summation, the term giving 
the increment of axial load at the centre of the m' strip 
due to shear on the m'" strip must clearly be halved. 
(This is also true of bending moment and web shear, 
but the increments are respectively trivial and zero.) 

It will be noticed that the shear at the cut end of the 
frame is not zero. This presupposes the existence of an 
edge member capable of carrying end load. In practice 
such a member will almost certainly be present, but if 
its cross-sectional area a, is small, its strain energy 
(given by L*x,*/(6Ea,)) should be included in the cal- 
culation. This will have the effect of reducing x,. 

Figure 21 shows the skin shear plotted against 4. 
Skin shear from engineering theory is included for com- 
parison. The effect of varying L is also shown. When 
L is halved, the maximum shear is increased by about 


10 per cent. When L is doubled or quadrupled, the 
maximum shear is decreased by about 9 per cent. and 
17 per cent., respectively. (For clarity, these latter 
values were not plotted.) The considerations of Section 
14(i), however, enable quite narrow limits to be put 
on the allowable variation of L, and the numerical 
results given thus support the statement in Section 14(i) 
that the maximum shear will not be very sensitive, in 
general, to variations of L within these limits. 


APPENDIX IV 


EFFECT OF THE INCLUSION OF THE STRAIN ENERGY OF 
THE SKIN AND STRINGERS IN END LOAD 

The calculation of Appendix III is modified as 

follows. Let t/=0-06 in. Assume that the total end 

load carrying area per strip is concentrated at the 


TABLE Il 
n=! 
26x $8 =2-7 
64 8 40x 8 
15:1 
Q = 88°7x, +51-8x 
64x 8 =3-2 40x 8 =2°0 
27°0x = 15-7 1S-7x =91 8:4 x 
27:4 
42-2 x 46 27:4x 46 16-2 16:1 
=51-8x, +31-7x 
Is? 15 =2-°7 5 
2700x 1S-7x 84x = 
42:2x 46 = 14:8 27 4x 46 9°5 16°1 x 
O =28-4x, Ix, 
5:8 5:8 
9-0 9-0 ; 
8:3 27°4x 46 5:36 x 
QO, =14-Lx, +8°74x 
6°5 6°5 
42-2 x 46 5-95 27:4 x 46 86 16:1 x 


) 
n 
7. — 
27: 
ds 5-8 x =5-8 
al- 9-0 x =§-25 6°5 x 5 
+14-Lx, +5-95x, 
ial 
= 
17 
the 58x =3-38 — 
27:4 27:4 
+8°74x, +3-86x, 
8:4 
Q 5-8 x 57 = 1°80 — 
90x =3-15 65x =2:27 
+4-95x, +2:27x, 
~ 
5-8 x B25 
9-0 9-0 
B +3-01x, +1:27x, 
: 
ae 
90x =1-27 65x =0 
0. =5-95x, +3-86x, +2-27x, +1-27x, +0-92x. 


JOURNAL OF THE ROYA 
junctions of the strips. Assume that the edge member at 
the cut-out is large enough for its strain energy to be 
negligible. 

The end load will build up at one “strip junction ” 
from zero to L(x,, — x,,-,) over a length L. 

The strain energy for one strip junction is thus 
L 


== 2°38 


In the x, region, the procedure is rather different. 
End load in a strip of width réz varies from zero at one 
end to 


x, (L sin (2+462)— Lsin z]=.x,L cos 6z 
(see Fig. 23) at the other end, giving the strain energy 
per elementary strip 


L 
(x,Lcosz6zdL  L*x,? cos? z6z 
2Er5zt’ 6Ert’ 


The total strain energy from end load in the x, 
region is thus 


47 


L*x," | 2 2 E sin 
6Ert’ cos* zdz=6°6 x 
= 4-34 x 10-*x,”. 
(The x, region has been taken as ending at 47° and the 
last 8° has been included with the x, - x, junction.) 
The total additional strain energy due to skin and 
Stringers is thus 


10 ‘x,?+2-38 x 10>* (0-777 x, — + 
+(x; +(%; — — 


0 


From this the increments of ¢CU/éx,, CU/ex,, ... . 
can now be calculated, but examination of Appendix III 
shows that what are really needed are the increments of 
(19-S Ox, (19-5) Ox, 


LOADING AS FIG. 14 


, WITH ALLOWANCE FOR 
z | (t= 55) STRINGER ENERGY 
FOR STRINGER ENERGY 

Zz 

0° B 
Ficure 22. The effect of including stringer energy in the 


calculation of the shear flow in Fig. 21 (See Appendix 1V). 


MARCH 1953 


L AERONAUTICAL SOCIETY 


Lx. sin of 


Fol 


tL 


sin & +5) 


Xe sine Hol 


Shear loading of a strip of skin in “x,~ 
(See Appendix IV). 


FiGuRE 23. region 


et cetera, and accordingly these are calculated. 


Incree CU _ 


= 135 x 10° x 4-76 x (x, x3) 
mentof (19-5)* ex, l x 10% x (x 


=6°4(x, 
= 12-8 (x,) 6-4(x, +,) 
10° = eu — = similarly 
(19°5)* dx, (19-5) ex, 
= 12-8 x, --6-4x,-4-98 x, 
(117+ 3-86) x, 


— 


The equations from Appendix III now become, with the 
direct additions of the terms just calculated, 


102-6x, +45°4x, +28°4x, +14:1x, + 5-95x; 1080, 
45-4x, +52:0x, +10°7x, + 8°74x,+ 3-86x, = 1260\ 
28°4x, + 2°77x, = 1345\ 
14-1x, + 8°74x,— 1:45x,+23°34x,— 5-13x, = 1326) 


5°95x,+ 3°86x,+ 2°27x,— 5°13x,+21-25x,— 4°98x,= 1205A 
— 4:98x, + 1180) 
0°799x, + 0:934x, + 0°996x, + 0°982x, + 0°891x, + 0-88x, = 25-7. 


The solution of these gives x, = — 2°37, x, =1-0, x,=4°5, 
x,—6°7, x,=7:65, x,=9-9, \=0-0652. (All the numeri- 
cal work is to slide rule accuracy only.) 


The shears are plotted graphically in Fig. 22. The 
effect of including stringer energy has been a slight re- 
duction in the maximum shear (i.e. from x,=10-0 to 
9-9). For more flexible frames, a much greater effect 
would have been obtained. See Sections 14(g) and 19. 

In the preceding, it has been tacitly assumed that 
after a length L the end load in the skin and stringers 
is redistributed in such a way that it becomes indepen- 
dent of x,, x,, and so on, and hence does not figure in 
the strain energy calculation. This redistribution could 
take place by bending the next frame out of its plane 
or by additional shear in the skins. The energy of these 
items has not been considered. It would constitute an 
additional complication, but its inclusion would be 
straightforward. It would increase the relief to the 
maximum shear. 
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DALLISON "STRESSING OF 
APPENDIX V 

SOLUTION OF A FRAME WITH A CUT-OUT AND INTERNAL 
REINFORCING 


(See Figs. 24 to 26) 


The centre line of the cross beam is defined by angle 
The edge of the cut-out is defined by angle 8. Other 
notation is as previously. 


Bending moment in the beam M, =500x—M,,. 
Strain energy of half beam = Upean, 
b 
(500 M,):d 


Bending moment at z,, due to X,, X, and M, is M’,, 


where M’ n— =X - X dm M.,. 

and Cn=rsin (y— 2), —r-cos (y — em). 

ax. =~ = | (provided that z,, > y). 


Ficure 24. Cut frame 

reinforced with floor 

beam (See Appendix 
V). 


The frame is divided into nine strips and / is taken as 
infinite for m=1. Total strain energy of the half frame 
in bending is 


n= 

Strain energy of the skin in shear is WG ot" 
Total strain neil U is given by 

3 
The equilibrium equations are 
n=1 
—X,cosy+X, sin y= 500 (b) 
Then, minimising by Lagrange’s method, 
b 
500 x— M,)dx "S (My +M’,,)=0 
0 
(c) 


“CIRCULAR FRAMES 


+M",) Cn = — COS 

(Mn+ MW) gin 


Equations (c) to (f) yield 12 equations in x,, x,...x,, X,. 
X,, and M, with terms in A and » and a constant on the 
right hand side. Solving in terms of and » and the 
constant, and substituting these values in (a) and (b) 
gives numerical values for A and » and also for x,. 


Figure 25. Frame of 

Fig. 24 showing forces 

applied to frame by 
floor beam (See 
Appendix V). 


Xo Xr 


APPENDIX VI 


(i) MINIMISING A FUNCTION OF SEVERAL VARIABLES 
USING LAGRANGE’S UNDETERMINED MULTIPLIERS 


and let x,. x,, X,,...%, be related by the expression 


When U has a stationary value (i.e. is a maximum or a 
minimum) it may be shown that:— 


oU 
ax, 
aU 
Ox, OX. 
aU 
Ox, OX, 


where \ is an “unknown multiplier.” 


50 1000 500 
FiGURE 26. Forces # b 
applied to floor beam 
by frame (See Appen- 
dix V). Mo Mo 


These equations are sufficient to determine x,. x,.... 
in terms of the unknown multiplier A. 
these values in X..... 


Xn 
Substitution of 
x,) yields and hence x,. 


If are related by two relations, e.g. 
and 


= (| 
Mo Mo 
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ex, Cx, OX, 
Ox, Ox, 


These equations yield x,, x,,...x, in terms of A and p. 
Substitution in @(x,, x.....x%,)=0 and u(x, x.....x,)=0 
gives A and and hence x,, x,.,...Xn- 


Example 
Suppose U=K,x,°+ (1) 


At a stationary value. 


cu cu 

2K ix. —a,’\=2K.x 
oU 

Bs. =@,’\=2K,Xx, ax, =a,\=2K,x,. 


Substitution in (2) yields 


Afa,? a," 


giving finally 


x,=ae/[K, K | 


K 


An example using both A and «u will be found in 
Appendix V. 


x, x 

1 Coefficients of (1) 3-16 4-28 

a 3:16 4:28 

2 Coefficients of (2) 4:28 3:46 
(1,2) 

b ) — 4-28 - 5-797 

c 2+5 0 -2°337 

3 Coefficients of (3) 5:37 1-07 
(3) .. 

d *™ 5237 =7 273 
(2-3) 
(2.2) xc) 0 6-203 

f 3+d+e 0 0 

4 Coefficients of (4) 2°10 —2:92 
(1,4) 

g “ait x (a) —2°10 —2:844 
(2, 4) 

h (2.2) x (c) 0 5-764 

(3. 4) 

i = (3.3) x (f) 0 0 


4+eth+i 0 0 
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(ji) CROUT’S METHOD FOR THE SOLUTION OF EQUATIONS 

This method can be applied to any equations but is 
specially suitable for arrays which are symmetrical 
about the leading diagonal. Details are too intricate to 
describe here, but can be found in Ref. 14. The method 
is best adapted for use with a computing machine. 


(iii) GAUSS’S METHOD FOR THE SOLUTION OF EQUATIONS 
(See Michielson, Journal of the Aeronautical Sciences, 
November 1949, p. 693) 

This method applies to arrays of equations which 
are symmetrical about the leading diagonal. It repre- 
sents a great saving of time and labour over the 
traditional methods. 

The table below is self-explanatory except for the 
following: —(2,3) means the coefficient of the third 
term in the second row. (This is normal matrix notation, 
the first suffix applying to rows and the second to 
columns). (3) means the complete third row of coeffi- 
cients, and so on. 


Equations for solution 
(1) 3-16x,+4-28 x, +5-37x,+2°10x,= 
(2) 2°18 
(3) 5°37x,+1:07 x,—3-05x,+1-62x,= 3-41 
(4) 2-10.x, x, + 1-62 x, — 4-85 x, = — 19-63. 


13 z 0 
F _ 134-770 _ 
rom (j) 33.593 =4-012 


x,=3-012 (from (f)) 
xX, =2:227 (from (c)) 
x,=0-692 (from (a)). 
Row (a) is apparently superfluous, being identical with 


Right hand 


Check 
210 51-23, 
5-37 2:10 36-32 51-23 
2:18 8-07. 
7.273 2844 49-193 69-387 
6-203 5-764  —47-013 61-317 
— 3-05 1-62 3-41 8-42 
9-126 3-569 61-721 87-059 
16-464 15-299 124-785 162-751 
4-288 13-350 66-474 84-112 
HBS 19-63 — 23-68 
15-299 14-216 115-953 151-233 
13-350 —41-563 ~206:956  —261-869 
0 —33-593  -—134-770  —168-363 
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FIGURE 27. Incom- 
plete bulkhead loaded 
with tangential or 
moment load on axis 
of symmetry (See 
Section 15(e)). 


Row I. It is introduced, however, to preserve the sym- 
metry of the notation. The check column is optional. For 
Rows 1, 2. 3. 4 it is obtained by adding the coefficients 
(including the constant term) in these rows. (If neces- 
sary the coefficients in any one column may be divided 
by 10 (or 100), to bring them to the same order as the 
other columns. This must then be done in every row). 
Operations carried out on the numbers in the rows are 
also carried out on the check column. From time to 
time a check addition is done to see that the sum of the 


coefficients in the row is still equal to the figure in the ~ 


check column. 
(iv) VON SEIDEL’S METHOD FOR THE SOLUTION OF 
EQUATIONS 

This method is applicable only when the coefficients 
of the leading diagonal are considerably larger than the 
other coefficients. The leading diagonal terms are 
retained on the left hand side of the equation and the 
other terms carried across to the right hand side. Taking 
X., X3, ¥, aS zero, a value of x, is calculated from the 
first modified equation. Using this value, and taking 
x, and x, as zero, a value of x, is next calculated from 
the second modified equation. Similarly with x, and 
X,, always using the latest available values of each un- 
known. Convergence is rapid if the given condition is 
satisfied. 


3:17 x, + 2°14 x, + 1:85 x, + 1-62 x, = 4:18 
2°14 x, + 4:23 x, + 1°64 x, + 1:21 x, = 3-17 
1-85 x, + 1-64 x, + 3:16 x, + 0-93 x, = 2-43 
1-62 x, + 1:21 x, + 0-93 x, + 2:94 x, = 1-16. 


X,=0°8085 —9-4139 x, —0-3578 x,—0-3133 x, 
X,=0°7494 —0°5059 x,—0-3977 x, —0-2861 x, 
X, =0°7690 — 0-5854 x, —0°5190 x, —0:2943 x, 
X,=0°3945 —0-5510 x,—0-4116 x,—0-3163 x,. 


x, 0°808511 0-696693 0-631610 0-611111 
x, 0°340375 0:416221 0:412692 0:400594 
x, 0-119000 0:212400 0:252083 0-263075 
x, —0°228677 -—0-227824 -0-203062  -—0-190265 


STRESSING OF CIRCULAR FRAMES 


Figure 28. Radially 
loaded frame cut to 
show redundancies— 
skin line offset from 
median line of frame 
(See Appendix 


A solution correct to one place of decimals was thus 
obtained in the third column, to two places in the fifth 
column and to three places in the seventh column. Had 
the approximate solution correct to one place been 
available, a solution to two places could thus have been 
obtained in about three columns, and to three places in 
about five columns. 


(v) SOLUTION OF EQUATIONS BY RELAXATION METHODS 

The method (due to Southwell) can be applied to 
equations of any type, not necessarily symmetrical about 
the leading diagonal. It is particularly adapted to 
obtaining a quick solution with the use of the slide rule, 
when an approximate solution is available. It can, how- 
ever, be used ab initio if no approximate solution is 
available. 

Using the same equations as in (iii) the solution ab 
initio is as follows. Choosing the variable with the 
largest coefficient in equation (1), and assuming the 
others are zero, write down its value. Applying this 
value to the other equations, write down the values of 
the left hand sides, and the residues obtained by sub- 
tracting these values from the right hand sides. 

The largest residue is then selected and the corres- 
ponding equation examined to find the largest coeffi- 
cient. A value is then assigned to the variable with this 
coefficient to make this residue smaller, and the other 
residues are correspondingly reduced (or increased). 

It can be shown that it is advisable to over-compen- 
sate, i.e. to select a value for the variable greater than 
that required to liquidate the residue. It therefore pays 
throughout to work to the next highest simple number 
(i.e. the next highest integer, tenth, hundredth, and so 
on). This makes computation extremely easy. 

The final values of the variables are obtained by 
adding together all the various values used. The 
example which follows will make this clear. 

Where an approximate solution is available, the 
calculation starts with this value and the work is greatly 
shortened. 


0-608176 0-610102 0:610561 0:610607 
0-394156 0-391671 0-391721 0-391795 
0:264369 0:263243 0:263043 0262994 
—0-186407 -—0-186090 -0-186267 —0-186340. 


P 
iS 
al | 
Vv 
BM, 


x,=0°9 
4 4-653 —0-473 0-642 —1-115 
1:926 +1:244 1-269 — 0-025 
Ny | 1:665 +0-765 0-492 +0:273 
| 1-458 -—0-298 0-363 —0-661 
x,=—073 x, =0:07 
x, | —0-486 —0-150 0-150 —0-300 
Xx, | —0°363 +0-274 0:296 —0-022 
—0:279 -—0-026 0-115 -—0-141 
x, | -—0°882 0-266 0-085 +0-:181 
x,= —0-03 x,= —0-02 
—0-155 —0-003 —0:037  +0-034 
Xs —0:064 +0:049 —0:033 +0-082 
—0-056 —0-067 —0:063 —0-004 
x; —0:049 +0-033 —0:019 +0-052 
x,=0-01 x,= —0-005 
¥; 0-016 —0-006 —0:009 0-003 
Xs 0-012 +0-001 —0:008 0-009 
Xs } 0-009 0-014 —0:016 0-002 
0-029 —0-008 —0:005 —0-003 


Final solution 
x,= 09-0-2 = 0-6l, 
x,= 0-34+0:07+0-02 = 
x,= 0-265, 
x,= —0°3+0:10+0-01 = -0:19, 


agreeing with the preceding solution. 


APPENDIX VII 


EFFECT OF STRINGER ENERGY ON THE SHEAR 
DISTRIBUTION AT A RIGID FRAME 


(See Section 19) 


From Appendix III, the strain energy in end load 
of skin and stringers (uniformly distributed round the 
frame) is 

L* "NS " 2 
When this energy is included, equation (6) of Section 14 
becomes, for a rigid frame 


n =P m = 
21G n=1 6Elt m=1 
or in calculus notation 


B B 


0 0 


with the condition from overall equilibrium (for a radial 
load P) 


B 
inzdz 
xX S81 Q 


0 


(2 being as in Fig. 27). 
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| 
—1:034 -0:081 | 0-555 —0-636 
—0:428 +0-403 0-492 —0-089 
-0:370 +0-643 0-948 0-305 
—0-324 -0-337 | 0:279 —0-616 
x,= —0:06 | x, 
—0:310 +0010 | 0-162 —0-152 
-0:128 +0-106 0-121 0-015 
—0-111 —0-030 —0-123 
—0:097 +0:278 | 0:294 —0-016 
x. =0-020 | x,= —0-010 
0:043 ~0:009  -—0-019 0-010 
0-085 —0:003 -—0016 0-013 
0-033 —0:037 ~-0:032 —0-005 
0:024 +0:028 | 0:009 0-037 


Moreover, if there is an end load carrying boom at £ 
in addition to the uniformly distributed end load carry- 
ing material, a term L*x,°/(6£A,) must be added to the 
strain energy, where Ax, is the end load carrying area 
of the boom and x, is the value of x at 2. 

Now x is a function of z, say F(z). The problem is 
to find F when U has a stationary value. Using the 
methods of the calculus of variations (see, for instance. 
Jeffreys and Jeffreys. Methods of Mathematical Physics, 
p. 315 et seq.):— 


8 2 
sy — | ] of ] 
E + dx dz\ép 
where p= sin z, 
= i L3 


and \ is undetermined. (Note: 6x, 6U represent small 
changes in U and x fora fixed z). Hence 
B B 
2K. pax | sin z 2K, + 
BEA, 
the last term being dropped when there is no boom. 
Now for a complete frame, ==. Also x has stationary 
values at both ends (z=0 and z=7=) and hence 6x=0 at 
z—0 and z=7. (The fixed values are both zero by 
symmetry). For a frame with a cut-out, 8 = = and two 
cases arise : — 
(a) There is no end load carrying boom at /. 
(b) There is an end load carrying boom at / in 
addition to the uniformly distributed end load 
carrying material, as already described. 
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In the first case x has fixed values (i.e. zero) at both 
ends (z=0 and z=). Hence 6x=0 at both ends. In 
the second case x has a fixed value (zero) at z=0, but 
at z=, x is indeterminate. In this case, it readily 
follows that at 


2K.p+ 3EA, 
i.e. that at 2=8, 
dx rx, 


dz  3EA,x2K, 


In all three cases it also follows that 


2K ,x—Asin2z—2K, =(), 
The general solution is, 
x— A, sinh z/(3B)+ A. cosh (3B)+ 


where B is K,/K, and is identical with the parameter B 
of Section 19. Applying the boundary conditions 
already discussed for the three types of frame con- 
sidered : — 


| | AL 
+ 
x | 
| 
A 
| | FOR ALL VALUES OF 
i 70 10 
KEY - A» 2x10? 
2x404 
BM=MAXIMUN FRAME BENDING MOMENT] | 
AL MAXIMUM FRAME AXIAL LOAD | 
WS = MAXIMUM FRAME WEB SHEAR a 
20 SS *MAXIMUM FUSELAGE SKIN SHEAR | — 
| | 
| 
| | (AXIAL LOAD+ 1-0 FOR ALL VALUES A) | 
Gers | | 
EIL | 
7 To 108 104 105 


Moment 


(BENDING MOMENT © 0 FOR ALL VALUES OF A) 


FIGURE 29. 


STRESSING OF “CIRCULAR ‘FRAMES: 


Complete frame 
A,=A,=0 and x varies as sin z, 


i.e. shear is sinusoidal. 


Frame with cut-out, no boom at free edge 


A,=0, A,= 


i.e. 
sinh z/(3B) 
sinh 8 J (3B) 


Shear is clearly far from sinusoidal. This case, however, 
is unlikely to occur in practice. 


x=C (sin 2~sin 


Frame with cut-out, free edge reinforced with boom 


—C(A,cos 8+?’r sin 8) 


A,=0, A,= (3B)cosh 8 (3B)+?rr sinh (3B) 
i.e. 
x=C( sin 


B+ rr sin B 


~ Any (3B) cosh BY (3B)+ t sinh BY (3B). ) 


EFFEGT [OF ECCENTRICITY ON MAXIMUM FRAME STRESSES 
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Gtr? | | 
19 102 103 104 
MomeNT. 
= 
= 
(BENOING MOMENT= 1-0 FOR ALL VALUES’ OF 7) 
Gtr } 
0 102 05 10+ 


FiGureE 30. 


q 
4 Res 
ee 
& 
EFFECT OF STRINGER FLEXIBILITY ON MAXIMUM FRAME (STRESS 4 | 
AND SKIN SHEAR | | 
BM (WS ETC) FOR A >2 x 106 | pak 
: | 
| | 
ll 
q 
i 
~ 
|__| | | | | 
| | | | | At 


MARCH 1953 


1% =VOL. 57 JOURNAL OF THE ROYAL AERONAUTICAL SOCIETY 


Since B is large (of the order of 100 upwards for prac- 4. Horr, N. J., Botey, B. A., and SALERNO, V. L. (1948). 
tical frames). for all but small values of 2. Shear Stress Concentration and Moment Reduction Fac- 
tors for Reinforced Monocoque Cylinders subjected to 


x=C (sin 2— ¥ Concentrated Radial Loads. Institute of Aeronautical 
; : Sciences, Preprint No. 151, July 1948. 
where is? Sin Ancos 8 5. BESKIN, L. (1946). Local Stress Distribution in Cylindrical 
VQB)Agttr Shells. Journal of Applied Mechanics, June 1946. 
/( +p?) 6. DuserG J. E.. and Kempner, J. (1947). Stress Analysis by 
é R e Formula of Reinforced Circular Cylinders 
Putting Ay=pt’r, gives €< a y 
GB) under Lateral Loads. N.A.C.A. T.N. 1219, 1947. 
: 5 7. Kempner, J. and DusBerG, J. E. (1947). Charts for Stress é( 
Taking B= 100, p=2, €<017 Analysis of Reinforced Circular Cylinders under Lateral 
Taking B— 100. p=3, &<0-065 Loads. N.A.C.A. T.N. 1310, 1947. 
8. Goopey, W. J. (1946). The Stresses in a Circular Fuselage. 
This gives for p=2, x=C(sinz—0-17) as z—>8 Journal of the Royal Aeronautical Society, November 
1946 
and fo =3, x=C (sin z—0-065) as 2— 8. 
‘ P C( si 9, Wise, J. A. (1939). Analysis of Circular Rings for Mono- 
' coque Fuselages. Journal of the Aeronautical Sciences, 
The second term in the bracket becomes smaller as z le 1939. } 
becomes smaller. Hence, in general, provided that the 10. Taytor, J. S. and Git, G. S. (1948). Analysis of a Circu- t 
boom area is reasonably large, the shear distribution lar Ring with Propped Floor Beam. Journal of the | on 
will be approximately sinusoidal, and the deviation from Aeronautical Sciences. April 1948. zx 
sinusoidal can be neglected. 11. Lewis, S. R. (1950). Diffusion of Loads in Non-rigid Cir- 0 
cular Frames. College of Aeronautics Report No. 33. t 
February 1950. te 
12. LANGHAAR, H. L. and Situ, C. R. (1947). Stresses in a 
REFERENCES Cylindrical Semi-monocoque Open Beams. Journal of the 
1. WiGnNot, J. E.. Comps, H. S., and ENsrupb, A. F. (1944). Aeronautical Sciences, April 1947. 
Analysis of Circular Shell-supported Frames. N.A.C.A. 13. Morse, W. L. (1952). Analysis of Fuselage Frames. Air- 
T.N. 929, 1944. craft Engineering, February 1952. 
2. Horr, N. J. (1944). Stresses in a Reinforced Monocoque 14. Crout, P. D. (1941). A Short Method of Evaluating 
Cylinder under Concentrated Symmetric Transverse Loads. Determinants and Solving Systems of Linear Equations 
Journal of Applied Mechanics, Vol. 11, No. 4, Dec. 1944. with Real or Complex Coefficients. Supplement to Elec- 
3. Horr, N. J. (1947). Thin-walled Monocoques. Anglo- trical Engineering Transactions, A.1.E.E., Vol. 60, Dec. 
American Aeronautical Conference, Sept. 1947, pp. 313- 1941. (Abridged as Marchant Method MM182, Sept. 1941. 
362. The Royal Aeronautical Society. Marchant Calculating Co.. Oakland, California.) 
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TECHNICAL NOTES 


The increasing flood of technical information makes it 
difficult to publish even abstracts of such information and 
it has become still more difficult to publish technical items 
of information which do not reach the status of full- 
length papers. 


Under this general heading of TECHNICAL NOTES it 
is hoped that members of the Society, and indeed all those 
who read THE JouRNAL, will help in contributing just those 
very items of useful information which normally might 
not see the light of the printed page. 


There will be no set form for these Notes. They may 
be illustrated or not. They can be short papers of a 
thousand or two thousand words; or abstracts and com- 
ments upon unpublished papers or upon published papers 
with a limited circulation; or notes of some interim results 
of research in hand and notes for further research; or in 
the form of a letter raising technical points or asking 
technical questions; or letters commenting on Notes 
already published. 


These pages are intended to be a technical forum where 
members can argue and discuss; where they can inform 
and be informed; and pass on their own experience, or air 
their technical difficulties. Contributions will be published 
as soon as possible after they are received and will be 
eligible for JOURNAL Premium Awards. 


No member need ever say that his particular subject 
does not get discussed in the pages of THE JouRNAL. He 
can always raise it himself, for these pages are his pages. 


Three Notes are published this month. A treatment of 
vortices by potential theory, a consideration of the velocity 
components for swept-back wings at incidence, which gives 
an example of the errors introduced by the usual practice 
as compared with the exact method of resolution; and a 
theoretical discussion of minimum comfortable cruising 
speed. It is felt in some quarters that a minimum com- 
fortable cruising speed for turbo-jet aircraft does not exist, 
and it will be interesting to have the comments and experi- 
ences of operators on Mr. Toms’ Note.—THE Epttor. 


Resolution of Velocity for the Swept-Back Wing at Incidence 


by 


C. H. McKOEN, Ph.D., D.I.C. 
(Sir W. G. Armstrong Whitworth Aircraft Co. Ltd.) 


HE COMMONLY ACCEPTED practice of resolv- 

ing a velocity V at incidence z into components 
V sin \ parallel to the leading edge and V cos A normal 
to the leading edge at incidence z sec .\ is only approxi- 
mately true and the exact method of resolution is 
developed here. 

In dealing with the aerodynamics of swept-back 
wings the free stream velocity may be resolved into a 
constant component parallel to the leading edge and a 
variable component normal to the leading edge. This 
is only possible if the constant percentage chord lines 
are all parallel to the leading edge, that is there is no 
taper, thickness taper or twist on the wing. This note 
gives the exact form and the small angle approximations 
of the velocity components and of the incidence of the 
component normal to the leading edge. 

There are two cases to be considered, the swept-back 
wing with and without dihedral, at incidence. There is 
no effect with dihedral alone as the incident velocity 
remains in the plane of the wing and may be resolved 
parallel to and normal to the leading edge in the plane 
of the wing. 


1. DEFINITION OF AIRCRAFT INCIDENCE 

Aircraft incidence is applied by rotation about an 
axis perpendicular to the free stream direction, that is, 
parallel to AB (see Fig. 1). Incidence cannot be applied 
to a swept-back wing by rotation about the trailing edge, 
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as this would imply breaking the wing along the plane 
of symmetry. 


2. RESOLUTION OF VELOCITY FOR THE SWEPT-BACK WING 
WITHOUT DIHEDRAL 


If z is the aircraft incidence and A is the angle of 
sweepback, then, using axes as in Fig. 1:— 

Velocity V resolves into 

(i) Vcosz parallel to the x-axis 

(ii) V sin parallel to the z-axis. 

V cos z component resolves into 

(i) VcoszsinA parallel to the leading edge 

(ii) VcoszcosA normal to the leading edge. 

Velocities V cos zcos A normal to the leading edge 
and V sin z parallel to the z-axis compound into 

V,=V ¥ (cos? zcos? A +sin? z) normal to the leading 
edge at an angle of incidence 8, where tan = 
tan zsec A (see Fig. 1). 

Thus velocity V resolves into 

(i) VcoszsinA parallel to the leading edge 

(ii) VJ (cos? zcos? A+sin* z) normal to the lead- 

ing edge at incidence 8=tan~' (tan asec A). 
For small angles of incidence this becomes 
(i) VsinA parallel to the leading edge 


(ii) VcosA normal to the leading edge at an 
incidence A. 
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IETY 


FiGure 2. Swept wing with dihedral at zero incidence. 


Velocity V at angle of incidence z resolves into 


(i) Vcosz parallel to the x-axis 
(ii) V sin parallel to the z-axis. 


V cos z parallel to the x-axis resolves into 


(i) Vcoszsin A parallel to the 
leading edge in the plane of 

(ii) Vcoszcos.\ normal to the { the wing. 
leading edge 


Ficure 1. Resolution of velocity for swept wing without 
dihedral at incidence. 
V sin z parallel to the z-axis resolves into 


3. THE EFFECT OF DIHEDRAL ON THE PLANE OF THE WING (1) V sin zsiny parallel to the 

Dihedral is put on a swept-back wing by rotation P leading edge in the plane 
about the line of symmetry. The chord lines in the (ii) Vsin zcosy normal to the { DEF of Fig. 2. 
downwind direction remain horizontal but any other leading edge 


line in the plane of the wing is not parallel to its original 
direction after dihedral has been put on the wing. We 
require to know the angle between the plane of the wing 
and the vertical plane through the leading edge for a 
swept-back wing with dihedral, at zero incidence. If 
this angle is 4, then in Fig. 2, angle DEG=8. 
DE is the normal to the leading edge in the vertical 
plane. 
GE is the normal to the leading edge in the plane of 
the wing. 
y is the angle of dihedral of the leading edge. 


From the triangle GED, 
GD? = GE? + 2GE ED cos 

or FG? - (FGsin A sec y)? =(FG cos A)? + 
(FG sin \ tan y)? — 2FG? sin A cos .\ tan y cos 4 


whence cos #=tan \tany. (1) 


4. RESOLUTION OF VELOCITY FOR SWEPT-BACK WING 


Ficure 3. Resolution of velocity for swept wing with 
Using axes as in Fig. 3: — dihedral at incidence. 


Recei 
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V cos zcos A normal to the leading edge in the plane 
of the wing and V sin zcos y normal to the leading edge 
in the plane DEF, compound into 

V,=V¥ {cos? zcos? .\ + sin? z cos” y 
—2 sin zcos zcos A cos y cos 6} 

=V {cos? zcos? A +sin* z cos? y - 
—2sin zcos zsin A sin y} 
normal to the leading edge at an angle of incidence 9, 


sin cos y sin 4 


where tang= — 

cos zcos .\ — sin z cos 7 cos 4 

Thus velocity V at incidence z resolves into 

V (cos zsin.\ +sin zsiny) parallel to the lead- 
ing edge 

(ii) V, normal to the leading edge at incidence 9. 


For small angles of dihedral » this reduces to 


(i) V (cos zsin.\+ysinz) parallel to the leading 
edge 
(ii) V/(cos* zcos* \ +sin® z- 2y sin «cos zsin A) 
at incidence 4, 
where 


tan o=tan zsec \ +y tan \ [tan z sec A]. 


TAYLOR 


5. NUMERICAL EXAMPLE 

To illustrate the application of these formulae, the 
magnitude and incidence of the velocity component 
normal to the leading edge have been calculated for a 
wing with 45° sweepback with no dihedral and with 5° 
dihedral at various aircraft incidences. These values 
are given in Tables I and II compared to the values 
obtained from the approximate expressions cos A and 
sec A, 


TABLE I 
MAGNITUDE OF VELOCITY COMPONENT 
cos .\ V,/V fory=0 for y=5° 
0 0-7071 0:7071 0:7071 
5 0:7071 0-7098 0-7022 
10 0:7071 0°7177 0-7027 
15 0°7071 0°7304 0-7086 
20 0-7071 0°7473 0°7197 
TABLE Il 
INCIDENCE OF VELOCITY COMPONENT 
B° for y=5° 
0 “0 ~ 
5 7-071 7-053 7-058 
10 14-142 14-002 14-152 
15 20°754 21-191 


20 28:284 27°236 28-055 


The Potential and Field of Flow (Induced Velocity) of a Plane Vortex 


J. LOCKWOOD TAYLOR, D.Sc. 
(Norges Tekniske Hogskole, Trondheim) 


ORTICES are not usually treated by potential 

theory, but the latter is so well known and of such 
general application that there seem to be advantages in 
doing so. 


Consider first the potential function 
(=) where r= (x7 4+ y? + 27) 


ie. distance from O. 


(-) 


FiGuRE 1. 


It can be shown by differentiation that this satisfies 


as required for irrotational (perfect fluid) motion. 


Also, along the x axis, +x=r, 


() . X positive 


. X negative. 


These are recognised as the potential of a straight vortex 
line, with the sign of the vorticity changing at O. 


Similarly along the z axis, +z=r 


x 
y 


(*) . negative 


i.e. there is also a vortex line along Oz. with change o 
sign at O. : 


| 
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The cruciform vortex represented by 


yr 


can be converted into an H-vortex by taking 


dx 
o=tan-?( ) -tan '(° 
vr yr 


( 2X - 
or — }tan a 6x = 5x 
Cx yr r(x?+y’) 
(Actually the two vertical lines of the H have coalesced 
into a line of doublets). 


| 
| 


> 


(-) (+) 


FIGURE 2. 


When a series of such short vortices, length 6x or 6s, 
are combined by placing them end to end in a plane 
curve (or straight line), a vortex line of any form is 
obtained. The vertical line vortices now cancel out 
(provided the curve closes) and we have finally 


x $s 
r(x? + 


o= 
replacing 6x by 4s to show it is measured along the curve 
and remembering that x and y are now tangential and 
normal to the curve at each point. and z is perpendicular 
to the plane of the curve. 
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The induced velocity at a point in the plane of the 
curve 


= now=x*+y° 
O27 


which is tne well-known Biot-Savart law. 
[Similarly for any other point if y is measured normal to 
és and r is the distance from 4s]. 

An alternative treatment, using the expression for 


$s 


FIGURE 3. 


the potential of a rectangular vortex derived in The 
Aeronautical Quarterly (November 1951), namely : 


o=tan tan 
Xa); 
+tan- ~tan ') 
Zr, 


where the suffixes refer to co-ordinates measured from 
the corners of the rectangle in the plane of the vortex, 
i.e. now using cruciform vortices in this plane. 


When the rectangle is small 


( xy 
Oxey ar 


where the double brackets indicate insertion of limits 
corresponding with the boundary of the area, which is of 
course the vortex line. 

This can be converted into a line integral round the 
vortex, but the expression previously given is probably 
more convenient. 


and in all 


HE MINIMUM comfortable cruising speed of an 

aeroplane depends on the complexity of charac- 
teristics for the design in question, which may be 
broadly grouped under the headings of stability and 
control, drag and propulsion. 

The store of knowledge relating to such speed is 
composed, in the main, from experience with aeroplanes 
powered by supercharged reciprocating engines, for 
which, at the low powers appropriate to the conditions 
considered, the brake horsepower and drag horsepower 
both increase with height at roughly the same rate 
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Minimum Comfortable Cruising Speed 


C. F. TOMS, B.Sc.(Eng.), A.F.R.Ae.S. 
(Percival Aircraft Ltd.) 


(constant aeroplane weight and equivalent air speed 
being assumed). For a turbine-engined aeroplane the 
brake horsepower (or thrust) decreases with altitude and 
it might be supposed that this difference would result in 
a difference between the minimum comfortable speeds 
of similar aeroplanes driven by turbine engines (either 
jet or propeller) and by reciprocating engines. It 
appears, however, that a much more important differ- 
ence arises from the change of thrust-speed charac- 
teristics as between jet and propeller-driven (whether 
turbine or reciprocating) aeroplanes. 

The minimum comfortable cruising speed may be 
broadly defined as a mean speed which, together with 
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TECHNICAL NOTES: Cr TOMS” 


the associated mean height, can be maintained under 
operational conditions, having regard to atmospheric 
and internally induced disturbances, with only moderate 
departures from the mean values and moderate times 
for recovery. 

The analysis which follows is applicable to flight 
under normal or automatic control since it is concerned 
with the recovery characteristics of an aeroplane on the 
assumption that the governing factor is the manner of 
the dependence of thrust and drag on speed and height, 
stability considerations in the classical sense being set 
aside. 


NOTATION 
v, rate of climb (ft./sec.) 
v forward speed (true) (ft./sec.) 
thrust (Ib.) 

D total drag (Ib.)=T 

W aeroplane weight = lift, L 
pair density (slugs/ft.*) 

p, Standard sea level air density 
o relative air density 

D, induced drag (Ib.) 

D, non-induced drag (Ib.) 


f 

b wing span (ft.) 
S wing area (ft.*) 
B 


=De: 


P., engine power at sea level (h.p.) 
T,, thrust at sea level (1b.) 
n exponent of o for power- or thrust-altitude 
relationship 
propeller propulsive efficiency 
h height (ft.) 
v dy 


ni- 


“md 


true speed for minimum drag (ft./sec.) 
K_ induced drag factor 

aspect ratio 

g acceleration due to gravity (ft./sec.*) 
vs true stalling speed (ft. /sec.) 

time (sec.) 


(7)... 


|. ANALYSIS 
The rate of climb of an aeroplane is given by 
v 
ve d(T -D) 


Ve T-D 
di 
or dv.= (dT - dD) +2, 


Hence 


~W 


Considering initially level flight («.—0), 
dv. dD 
dv dz 
Now 
D=D,+D,=fep,v? + - 
TPyv 
dD _ dD, aD, 
dv dv dv 
aD, do dv 
and since 
o v 
dD,__ do _ dv 
dD 2 Ido 
2 1 dedh 


Considering the idealised cases of propeller-turbine 
engines (with negligible jet thrust) and pure jet engines 
(with thrust independent of speed) separately : — 
Propeller-turbine 


de dn du 
om Vv 


aT do dh. (3) 
dv. v 
Pure jet 
T= 
dT nDdeo dh 
dv dhdv (4) 


Now in the standard atmosphere, over a wide height 
range, 


de 
34x 10-°. 
6) 
Also D, — 
leading to D, - Di= (= (6) 
Again (7) 


A combination a slight rearrangement of equa- 
tions (1) to (7) gives the following expression for the 
derivative of 7. with respect to v:— 


Propeller-turbine 
dv, 


dy 3(7) 


+34~x 10-' 


Pure jet 


As equation (8) except for the absence of the term 
It is seen that a variable dh/dv remains and in con- 
nection therewith it is necessary to consider, the kinds 
of disturbance which an aeroplane may suffer; here a 


(8) 
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major difficulty arises. For the ideal conditions of still 
air the values of dh and dv in a small disturbance will 
be related in a simple manner. The possibility of small 
perturbations in height at constant air speed, and vice 
versa, in a disturbed atmosphere is obvious. A more 
remote, but not negligible, possibility is that of speed 
and height changes of the same sign occuring together. 
This latter possibility does not lend itself to any speci- 
fication other than an arbitrary one, and will not be con- 
sidered further here. 
For still air conditions and small disturbances it is 

readily shown that 

dh v 

and substitution in equation (8) gives, very approxi- 
mately 


min 


+. m'— | (2- | \ 


for propeller-turbine engines. 

The value of dz../dz for given values of m, 7, z and n 
is directly proportional to (D/L),,in. Thus it is implied 
that. of two aeroplanes operating at a given speed the 
less aerodynamically efficient will be the more “ speed- 
and height-stable.” Since propeller-turbine and jet 
aeroplanes are respectively somewhat, and considerably, 
more efficient than those powered by reciprocating 
engines (aerodynamically) it is implied that they are 
on this account the less stable in the present sense. other 
conditions being the same. 

Since for a given aeroplane 7 varies with m, when 
evaluating equation (11) (and its counterpart for jet 
engines) #z,,, is submitted for 7, and various values of 
7,7 Can then be considered. 

Figure | shows the variation of dv../dv with m from 
equation (11) for two values of (D/L),,in, two values of 
7, and of three values of n (- 0-5, 0:5 and 1-0; the first 
of these represents approximately the conditions for a 
supercharged reciprocating engine at constant r.p.m. 
and manifold pressure and the latter two the probable 
range of values for turbine engines). Also shown is 
dv../dv from equation (11) after putting «*=0. It is 
seen that these terms are relatively unimportant and 
produce only small variations in the value of dz./dv. 

To complete the presentation, equation (11) is 
plotted for two values of (1 - z) and the corresponding 
equation for jet engines for 7?=0 only. 

Since putting 7«°—0 introduces little error, and this 
amounts to putting dh/dv=0, it follows that the 
behaviour at low speeds is not significantly affected by 
the relationship between dh and dz in a disturbance, at 
least between the limits dh/dv =O and -7/g, and hence 
equation (11) covers the case of independent distur- 
bances to / and 7. 


2. CONCLUSIONS 

On the assumption that the motion of an aircraft 
after disturbance at low speed conditions is governed 
only by the variations of thrust and drag with speed and 
altitude (i.e. that inertia and longitudinal stability in the 
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FIGURE 1. 


classical sense are unimportant factors) it is shown that 
the minimum comfortable speeds for propeller-driven 
aeroplanes are substantially independent of the type of 
engines installed, but that, owing to the changed thrust- 
speed relationship, jet aeroplanes have slightly higher 
minimum comfortable speeds. It is also shown that the 
“aerodynamic efficiency” of an aeroplane, as defined 
by (D/L)nin has a noticeable effect on the speed in 
question: an increase in the “aerodynamic efficiency ~ 
(reduction in (D/L),...) increases the minimum speed. 

If the negative value of dv./dv which represents 
acceptable low-speed characteristics for reciprocating- 
engined aeroplanes be taken as that corresponding 
to m= 1-05, then the corresponding values for a turbo- 
prop aeroplane are about 1:10 (allowing for a lower 
value of (D/L),i,) and for a jet aeroplane about 0-2 
greater. 

It is advisable to check that the selected speed does 
not too closely approach the stalling speed; it can be 


shown that 
vy 1:33 ) 
\K 


It is appreciated that other considerations may con- 
tribute to a selection of a suitable value of m for flights 
of long duration; comments and information from 
manufacturers and operators of jet aeroplanes in par- 
ticular would be welcomed. 
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STRUCTURAL PRINCIPLES AND DATA. Handbook of 
deronautics, No. 1. Pitman. 1953. Published under the 
authority of the Royal Aeronautical Society. Series Editors: 
C. G. Burge and J. L. Pritchard. Co-Editors, Structures: 
D. M. A. Leggett and M. Langley. Fourth Edition. 322 pp. 
/Ilustrated. 45s. net. 

The writing of a book is an affair of personal aim and 
personal effort: team work may be good for the schoolboy 
and sometimes, perhaps, for the scientist, but when it 
comes to writing a book, a team of authors is often in 
difficulties from the start. A realisation of this may be 
common ground in older fields of authorship, but it is 
currently viewed in aeronautical circles as slightly heretical 
or at least a little unsporting. Not least among aeronautical 
efforts at team writing was the first edition, in 1931, of 
the Handbook under review, and it was indeed the study 
of this at the time that first revealed the difficulties of the 
process to the present reviewer. 

Then, as no doubt now, the editors did their best to 
enrol eminent and able authors and to instruct them in 
the aim of the work and their part in its achievement. 
But, nevertheless, the authors produced contributions that 
were an expression of themselves and their aeronautical 
interests without much reference to the rest, or else they 
never really put their heart in the matter and contented 
themselves with piling table upon table of rapidly dated 
data. The best contributions—those that have lasted in 
value—were those that took no notice of the others, a 
commonplace to anyone who cares to study scientific 
articles in early editions of the Encyclopedia Britannica. 

This new volume, the first of several that are to 
constitute a 4th edition of the original work, opens with a 
list of four very well known and able editors. One 
wonders, is there to be a new and stronger effort to co- 
ordinate the authors? But no, the two parts of this book 
have as their authors men who cannot be expected to 
submit to a surfeit of co-ordination: let us defer judgment 
till we have studied each part. 

Part I, with the title * Structural Airworthiness,” is by 
W. Tye, the Chief Technical Officer of the Air Registration 
Board. As one would expect from its title, this part deals 
primarily with the loads that arise on aeroplane structures 
and the specification of values for these to ensure safe, 
yet efficient, design. The chapters deal in succession with 
manoeuvring loads, gust loads, tail loads, aileron loads, 
undercarriage and engine mounting loads, crash loads and 
Structural safety questions in general. The whole form 
an admirable, though condensed, statement of current 
airworthiness methods in this field, the best chapters, 
naturally encugh, being those dealing with recent develop- 
ments and the least satisfactory. perhaps. that on engine 
mountings, an almost static art. 

Part II, with the title “ Structural Analysis,” is by J. H. 
Argyris, now of Imperial College. and P. C. Dunne, both 
of whom worked together until recent years on the 
Stressed-Skin Data Sheets of the Royal Aeronautical 
Society. This part deals mainly with »svects of structural 


analysis arising in modern aircraft structures and naturally 
assumes a background of engineering knowledge of roughly 
degree standard. Most of the chapters are concerned with 
basic information, with titles referring to stress-strain 
relations, beams and struts, buckling of plates. tension field 
shear panels, stress diffusion, and curved beams and rings, 
and only two deal specifically with larger aircraft com- 
ponents, entitled Wing Stressing and Fuselages with 
Flexible Rings.” The whole is followed by an extensive 
bibliography. The approach throughout is the critical 
scientific one typical of the authors, and no condensation 
of the matter is allowed to interfere with this. 

Looking now at the two parts of this book in the light 
of our opening remarks, it is clear that, whatever the 
intention, the authors have nor been unduly co-ordinated, 
and this the reviewer regards as, in the long run, a good 
thing. But as a result, and because of the different natures 
of the subjects and authors, the two parts will probably 
have different real uses. By its nature, no one can write 
a reference handbook on airworthiness that will long be 
used as such by designers; inevitably, their only safe hand- 
book is that giving the official up-to-date requirements 
themselves. Part I must in the long run have educational 
rather than design uses, and for this it is admirable. 
Part II is also not really a section of a design handbook; 
it is a most interesting, critical text book in a specialist 
field and brings together for the first time much recent 
work in the Argyris-Dunne tradition, most of it, naturally 
enough, relating to fundamental rather than empirical 
knowledge. 

The whole is well and accurately printed, and though 
the price is rather high, the publishers are to be congratu- 
lated on avoiding any competition with current American 
prices.—A.G.P. 


STRUCTURE OF METALS. C. S. Barrett. McGraw-Hill. 
1952. 2nd edition. 661 pp. Diagrams. 72s. 6d. net. 


The appearance of a second edition of this book, first 
published over nine years ago, is evidence that although 
its scope is limited, it satisfies a definite need. The author, 
a physicist of the highest erudition, described his original 
work as being intended to serve both as a text and as a 
reference book. This description applies equaily to the 
present edition, which, although nearly one hundred pages 
longer, is largely identical with its progenitor. 

To the engineering metallurgist, structure has a wide 
significance. Such a man thinks of it as relating to the 
build-up of a casting, a forging or of an extrusion—that 
is to say, to a complete metallurgical entity. To Professor 
Barrett, however, the structure of metals is largely a 
question of the arrangement of a relatively small number 
of atoms in a unit often not so large as a single crystal. 
Since, however, in one of iron, having a volume of one 
cubic inch, there would be approximately 1-4 x 10** atoms, 
one sees that the author’s field is often a restricted one. 

To a reader whose primary interest is limited to such 
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a microcosm, this carefully written and well reproduced 
work with its excellent descriptions of the applications of 
many X-ray, electron, atom and neutron techniques, will 
be of great value. On the other hand, to those whom one 
suspects of being more numerous, who look upon the 
structure of a single piece of metal as involving a 
heterogeneous conglomeration of allotriomorphic crystals, 
differing in size, orientation and chemical composition, 
containing physical discontinuities varying from Griffith 
cracks to substantial cavities and chemical ones, such as 
metallic and non-metallic inclusions, strangely retained 
gases and residual stresses, this scholarly work may be 
something of a disappointment. 

By the selection of his title, Professor Barrett has let 
himself in for a certain amount of criticism—what he has 
done is analogous to that of an author whose book, entitled 
* The Structure of Houses ™ is, in fact, a learned essay on 
the making of bricks. with or without straw.—P. L. TEED. 


AIRCRAFT INSTRUMENT DESIGN. W. H. Coulthard. 
Pitmans, London 1952. 309 pp. Illustrated. Index. 40s. net. 


The technique of aircraft instrumentation is advancing 
with an ever-increasing tempo in an endeavour to keep 
pace with the aircraft themselves. It may even be suggested 
that increased endeavour is essential to keep in the race. 
Thus, nearly a quarter of a century after the appearance 
of * Stewart's Aircraft Instruments,” a comparable volume 
is more than overdue. Mr. Coulthard is to be congratu- 
lated on his attempt to fill the void and it is no credit to 
British industry that both these books were left to the 
initiative of members of the Instruments Department of 
the R.A.E. 

Mr. Coulthard has presented his work in three main 
parts dealing respectively with the instruments designed 
for use by the Pilot, the Navigator and the Flight Engineer: 
a fourth part, written by Mr. W. O. Broughton, deals with 
Automatic Flight treated as an extension of the classic 
theory of aeroplane stability. 

The pattern followed generally throughout the book is 
to survey in turn each of the important quantities and 
establish the physical bases of its measurement, after which, 
in each case, the author illustrates the principles by 
descriptions of typical modern instruments. 

It is perhaps a disadvantage of this microscopic 
approach that the broad principles underlying the success- 
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ful design of aircraft instruments are never touched upon. 
Surely, in a book of this class entitled * Aircraft Instrument 
Design,” the reader might expect some guidance on display 
such as has been provided in the * Naval Motion Study 
Unit Report No. 47.” More generally still, one might 
expect some classification of displays according to the 
purpose served. Some pilot’s instruments are essentially 
mis-alignment indicators in closed servo loops of which 
the pilot is expected to act as the servo mechanism. For 
too long the instrument designers have been laying a 
macedoine of information in front of the harassed pilot 
and patting themselves on the back every time they added 
another plum. However. during the past decade, the 
designers have been showing increasing awareness of the 
necessity of making their servo loops stable, instead of 
depending upon the skill of the pilot in emulating the 
stage performers who can balance half a dozen unstable 
systems simultaneously. The Sperry “ Zero Reader” is 
a well known example of this modern trend, but the reader 
will search in vain for any mention of this instrument or 
for any mention of the principle on which it is based. 

A!though an interesting analysis of the errors developed 
by vertical gyros during turning flight is included and data 
relating to the most modern gyro-magnetic compass show 
that compass errors of over 20° can be caused by turning 
even at subsonic speeds, the point is not stressed that 
gravity-monitored instruments will have to improve very 
fast if they are to keep pace with future aircraft. In fair- 
ness to Mr. Coulthard it must be said that he may have 
been handicapped by considerations of security. 

No such excuse can be offered for a general looseness 
of expression which is likely to cause either irritation or 
slovenly thinking, according to the attitude of the reader. 
Why should the word “radiation ™ be used to embrace 
the total loss of heat from a heated pitot head and what 
does “centrifugal acceleration” mean? The latter ex- 
pression is used not to describe the behaviour of the 
components of an atom bomb but to account for the 
northerly turning error of a compass. Mention must be 
made of a palpably incorrect formula for the distance 
between two points on a tangent to the earth’s surface 
because innocent students could waste much time making 
false calculations based thereon. 

The indexing can only be described as lamentable. In 
spite of its shortcomings, this book contains much valuable 
material and deserves to be widely read.—F. W. MEREDITH. 


Additions to the Library 


Brennan, D. T1iMeE ENOUGH TO Live. Secker & Warburg. 
1953. 

Brown, Winifred. No Distress SIGNALS. Peter Davies. 
1952. 

British Electrical Development Association. HIGHER 
INDUSTRIAL PRODUCTION WITH ELECTRICITY. 1953. 
Chestnut, H. and R. W. Mayer. SERVOMECHANISMS AND 
REGULATING SYSTEM DesIGN. Vol. I. John Wiley. 

1951. 

Dommasch, D. O. 
Pitman. 1953. 
Hunsaker, J. C. AERONAUTICS AT THE MIpD-CENTURY. 

Yale University Press. 1952. 

Illing, V. C. (Editor). THE SCIENCE OF PETROLEUM. Vol. 
VI. Pt. 1. THE Worvp’s OILFIELDS: THE EASTERN 
HEMISPHERE. O.U.P. 1953. 

Kruschik, J. Die GASTuRBINE. Springer. 1952. 

Lilley, T. et al. PROBLEMS OF ACCELERATING AIRCRAFT 


PRINCIPLES OF AERODYNAMICS. 


PRODUCTION DURING WorRLD War II. Harvard Univer- 
sity. 1947. 

Muchmore, R. B. ESSENTIALS OF MICROWAVES. John 
Wiley. 1952. 

National Physical Laboratory. INDUSTRIAL APPLICATION 
oF AERODYNAMIC TECHNIQUES. (Notes on Applied 
Science No. 2). H.M.S.O. 1952. 

Oswatitsch, K. GASDYNAMIK. Springer. 1952. 

Sechler, E. E. ELAsticity IN ENGINEERING. John Wiley. 
1952. 

Society of British Aircraft Constructors. | STANDARD 
METHOD OF PROPELLER PERFORMANCE ESTIMATION. 

Urry, S. A. SOLUTION OF PROBLEMS IN STRENGTH OF 
MATERIALS. Pitman. 1953. 

Wiener, N. CYBERNETICS. John Wiley. 1948. 

Williamson, G. W. Some PRACTICAL STEPS TOWARDS 
INCREASED PRopuctivity. Belfast Association of 
Engineers. 1952. 
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Reports 


AERODYNAMICS 
BOUNDARY LAYER 


investigations into the effect of continuous suction on laminar 


boundary-layer flow under adverse pressure gradients. B. 
/hwaites. R. & M. 2514. 


The principal problem considered in this paper is that of 
the flow along a surface through which fluid is being 
continuously withdrawn at a constant velocity, the velocity 
of the stream outside the boundary layer having a constant 
and negative gradient. When no suction is applied Howarth 
has solved the flow, and the point of separation is known. 
With continuous suction, the ordinarp methods of using the 
momentum equation, for example Pohlhausen’s method, 
break down. and a new method is used in this paper to find 
the point of separation when suction of constant velocity 
is applied. Only partial success may be recorded in this 
problem, but this account of progress so far is now given 
in the hope that the problem may attract other workers. 
The momentum equation of boundary-layer flow is also used 
to deduce other types of flow in which the velocity of suction 
or the velocity outside the boundary-layer is not constant. 
These examples are inserted by virtue of the simplifications 
of method involved, by which the momentum equation may 
be integrated exactly; one of them is particularly applicable 
to the velocity distribution near the leading-edge of a thin 
aerofoil at high lift-coefficient, in which a small amount of 
suction is sufficient to prevent the stall.—(1.1). 


The asymptotic theory of boundary-layer flow with suction. 

Three parts. E.J. Watson, R. & M. 2619. 
The subject of this report is the steady two-dimensional flow 
of a boundary layer over a permeable surface through which 
the fluid is withdrawn at a known rate of suction. This rate 
of suction is assumed in accordance with the hypotheses of 
the -Foundary layer, to be small compared with the stream 
velocity. Part I deals with the similar solutions of the 
boundary-layer equations, Part II with an arbitrary pressure 
distribution but constant suction velocity. and Part III with 
the general problem, Attempts are made in both Parts I and 
Il to find when separation occurs, but only rough estimates 
can be made as the series do not converge well. In Part II 
the theory is applied to the flow over a porous circular 
cylinder in a uniform stream, and also to the use of suction 
round the nose of an aerofoil to prevent stalling at high 
incidence.—-‘1.1). 


Generalization of boundary-layer momentum-integral equations 
io three-dimensional flows including those of rotating system. 
. Mager. N.A.C.A. Report 1067. 


Boundary-layer equations for application in three-dimensional 
flows are developed. With the use of a fixed velocity profile 
and an empirical friction law an approximate solution in 
closed integral form is obtained for a generalised boundary- 
layer momentum-loss thickness and flow deflection at the 
wall in the turbulent case.—(1.1). 


Survey of transition point measurements at the N.L.L., mainly 
for two-dimensional flow over a N.A.C.A. 0018 profile. H. 
Wijker. N.L.L. Amsterdam, Report A.1269. 


This report contains some remarks on experiments in the 
determination of the transition point. Results following 
Various techniques are compared. Transition points on a 
N.A.C.A. 0018 profile are given for an angle of incidence 
between —6° and +10° and Reynolds numbers between 
0-4 10° and 2-9 x 10°.—(1.1). 


Some experimental investigations on the influence of wall 


houndary layers upon wind tunnel measurements at high sub- 
sonic speeds, E.G. M. Petersohn. F.F.A. Sweden. Report 44. 
Pressure distribution measurements and drag determination 
by means of balance measurements have been carried out 
for a number of models at high subsonic velocity in wind 
tunnels. where the boundary layer of the walls has been 
varied. The phenomenon described should be of a certain 


NOTE:—The figures in parenthesis at the end of each Summary are for office use only. 


importance from the point of view of wind tunnel technique. 
since it is possible to increase the choking velocity for a 
given model by means of thickening the boundary layer.— 
(1x 1. 2). 


Approximate calculation of the influence of wall boundary 
layers upon the blockage interference in a high speed wind 
tunnel. S. B. Berndt, F.F.A. Sweden. Report 45. 


The wall boundary layers in a subsonic wind tunnel work 
to reduce the choking effect of a model. A simple approach 
to a theoretical explanation is used. It involves neglecting 
viscosity, linearisation of the equations of motion and sub- 
stitution of wall layers of constant velocity for the boundary 
layers. Solutions are given for plane and axisymmetric flow. 
Results of numerical application to the case of a parabolic 
arc profile in plane flow confirm experimental results and 
indicate strong boundary layer effects to occur at Mach 
numbers close to one.—(1.1 « 1.2). 


Tests on four aerofoil cascades: Part I. Boundary layer 

characteristics. F.G. Blight and W. Howard. A.R.L. Australia, 

Report E.75. 
Tests on four cascades, described in Part I of this report, 
indicated the presence of a boundary layer separation on 
the blade upper surfaces. The nature of this separation has 
been investigated using oil film, surface pitot and hot wire 
methods, the most consistent and satisfactory results being 
obtained with the oil film.—(1.1 « 1.5). 


Some investigations on thin nose-suction aerofoils, Parts I and 
I]. J. Williams. R. & M. 2693. 


The stalling properties of some thin nose-suction aerofoils 
already tested have been examined, and further theoretical 
investigations have been carried out on thin aerofoils 
specially designed to give high lift with nose-slot suction. 
In Part I, the experimental results from stalling tests on 
thin nose-suction aerofoils are compared and the design 
features of the tested aerofoils are analysed. Part II 
describes a theoretical exploration of possible thin nose-slot 
aerofoils specially designed to have an abrupt fall in velocity 
where suction is to be applied on the upper surface of the 
nose. The theoretical formule required for the calculation 
of the velocity distritution and lift of an aerofoil with a 
sink on its surface, and some detailed notes on the design 
of the Glauert nose-slot aerofoil. are given in the Appendices. 
x 1.10). 


COMPRESSIBLE FLOW 


See also RESEACH 
Experiments on transonic flow around wedges. G. P. Wood. 
N.A.C.A. Technical Note 2829. 


Several aspects of transonic flow around the forward portions 
of wedge profiles were studied. by means of interferometry. 
Measurements were made of the two kinds of flow pattern 
that occur at the leading edge of a wedge at an angle of 
attack. The growth of the supersonic region at a sharp 
convex corner formed by two flat surfaces was observed. 
The pressure-drag coefficients of a wedge of 14°5° semi- 
angle were measured at Mach numbers of 0-768, 0°819, and 
0-854, and were shown to be consistent with those of wedges 
of smaller angle when plotted according to the transonic 
similarity law. Conditions at the bases of the shock waves 
that interacted with boundary layers on the wedge were 
measured. The method of characteristics was used to 
calculate the flow behind an experimentally determined sonic 
line, and the calculated flow field was compared with the 
measured flow field. The accuracy in the location of the 
sonic line necessary to give correctly the pressure distribution 
on the surface behind it was determined.—(1.2). 


Laminar boundary layer on cone in supersonic flow at large 
angle of attack. F.K. Moore. N.A.C.A. Technical Note 2844. 


The laminar boundary-layer ‘flow about a cone at large 
angles of attack to a supersonic stream has been analysed 
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in the plane of symmetry. At the bottom of the cone. 
velocity profiles were obtained showing the expected 
tendency of the boundary layer to become thinner on the 
under side of the cone as the angle of attack is increased. 
At the top of the cone, the analyses failed to yield unique 
solutions, except for small angle of attack. Beyond a certain 
critical angle of attack, boundary-layer flow does not exist 
in the plane of symmetry. thus indicating separation.— 


Shock tube theory and applications. J. Lukasiewicz. N. A. E. 
Canada, Report 15. 


The theory of shock tube flow is presented in a form con- 
venient for aerodynamic design of shock tubes. Tables and 
graphs of the required functions, covering the whole range 
of values of the specific heat ratio, are included. Application 
of shock tubes to investigation of compressible flow is 
examined. In particular, the use of shock tubes as sub- 
and supersonic tunnels, as aeroballistic ranges (including 
hypersonic speeds) and for investigation of propagation of 
continuous waves. shocks (including variable specific heat 


and relaxation time effects) and interactions is considered.— . 


31.2): 


FLuip DyNaMIcs 


Corrections for drag, lift, and moment of an axially sym- 
metrical body placed in a supersonic tunnel having a two- 
dimensional pressure gradient. J. Kolodner, F. Reiche and 
H. F. Ludloff. N.A.C.A. Technical Note 2837. 


Corrections for drag, lift. and moment are derived for an 
axially symmetrical body placed in a supersonic tunnel 
having a two-dimensional pressure gradient.—(1.4). 


Low aspect ratio wings with small thickness at zero lift in 

subsonic and supersonic flow. F. Keune. K.T.H. Sweden. 

Aero T.N. 21. 
The flow considered can be described by two terms: a cross- 
sectional flow orthogonal to the direction of flight and a 
spatial influence. The first term is a two-dimensional flow 
at any distance from the leading edge of the wing, corres- 
ponding to Jones’s theory. The second term depends on the 
variation of the cross-sectional area in chordwise direction, 
analogous to the flow around bodies of revolution.—(1.4). 


INTERNAL FLOW 


Flow surfaces in rotating axial-flow passages. J. D. Stanitz 
and G. O. Ellis. N.A.C.A. Technical Note 2834. 


In order to investigate the deviation of flow surfaces from 
their assumed orientation in the usual type of two- 
dimensional solution, three-dimensional. incompressible. 
non-viscous. absolute irrotational fluid motion is determined 
for flow through rotating axial-flow passages bounded by 
straight blades of finite spacing and infinite axial length 
lying on meridional planes.—(1.5). 


Effect. of changing passage configuration on internal-flow 
characteristics of a 48-inch centrifugal compressor. 1l1—change 
in hub shape. J. Mizisin and D. J. Michel. N.A.C.A. Technical 
Note 2835. 
The passage contour of a 48-inch centrifugal compressor was 
modified by changing the shape of the hub to control the 
deceleration rates along the blade surfaces in order to 
improve the internal efficiency of the impeller. A com- 
parison of internal-flow characteristics at design flow rate 
was made with the original impeller and with a modified- 
blade impeller that had the same area variation in the 
passage. In addition, flow characteristics of the modified- 
hub impeller over a flow range from maximum flow to 
near surge at a corrected tip speed of 700 feet per second 
are presented.—(1.5). 


Etude générale de l'écoulement d'un gaz a travers une tuyére 
quelconque et du passage par la vitesse du son. M. Serruys. 
Publications Scientifiques et Techniques du Ministére de L'Air, 
France, No. 272.—(1.5). 


Die Verteilung der Durchflussmenge in einem ebenen Verzwei- 
gungssystem. H. Reichardt und W. Tollmien. Max-Planck- 
Institut, Géttingen, Mitteilungen Nr. 7.—(1.5). 
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STABILITY AND CONTROL 


The static longitudinal stability and control of an aeroplane as 
affected by the compressibility of the air. J. Buhrman and 
C. M. Kalkman. N.L.L. Amsterdam. Report V.1625. 


A computation method is given for the effect of the com- 
pressibility of the air on the static longitudinal stability and 
control characteristics of an aeroplane. The derived formule 
are based on the Prandtl rule. Their applicability is limited 
to Mach numbers below the critical Mach number since the 
linear theory holds no longer beyond this limit. In order 
to illustrate the magnitude of the influence of the Mach 
numter on static and manceuvre margins and trim changes 
the results of a numerical application of the formule are 
included in the report.—(1.8). 


WINGS AND AEROFOILS 


The scope and accuracy of vortex lattice theory. V. M. Falkner. 

R. & M. 2740. 
The report gives an outline of the development of the 
principles on which potential problems in lifting-plane theory 
are solved by the use of a vortex lattice for the purpose of 
computing downwash. The conditions of convergence 
necessary for an accurate solution are defined, and the main 
purpose of the report is to show that those connected with 
the lattice have been. or can easily be. satisfied.—(1.10). 


Langley full-scale-tunnel investigation of the maximum-lift and 
stalling characteristics of a trapezoidal wing of aspect ratio 4 
with circular-arc airfoil sections. R. H. Lange. N.A.C.A. 
Technical Note 2823. 
Results are given of an investigation at high Reynolds 
numbers and low Mach numbers to determine the maximum- 
lift and stalling characteristics of a trapezoidal wing of 
aspect ratio 4 with 10 per cent-thick, circular-are aerofoil 
sections. The tests included measurements of the lift. the 
drag, and the pitching-moment coefficients of the basic wing 
and of the wing with 0:20-chord droop-nose and trailing- 
edge flaps deflected both alone and in various combinations. 
Scale effects were investigated at Reynolds number ranging 
from 3-27 x 10° to 7°67 x 10°.—(1.10). 


Effects of independent variations of Mach number and Reynolds 
number on the maximum lift coefficients of four N.A.C.A. 
6-series airfoil sections. S. F. Racisz. N.A.C.A. Technical Note 
2824. 
The N.A.C.A. 65-006. 64-009. 64-210, and 642-215 aerofoil 
sections were tested in the Langley low-turbulence pressure 
tunnel to determine the effects of independent variations of 
Mach number and Reynolds number on the maximum-lift 
characteristics for the smooth and rough conditions.— (1.10). 


Effect of a finite trailing-edge thickness on the drag of rect- 
angular and delta wings at supersonic speeds. E. B. Klunker 
and C. Rennemann, Jr. N.A.C.A. Technical Note 2828. 


The effect of a finite trailing-edge thickness on the pressure 
drag of rectangular and delta wings with truncated diamond- 
shaped aerofoil sections with a given thickness ratio is 
studied for supersonic Mach numbers, linearised theory 
being used to evaluate the surface pressures. In order to 
facilitate comparison with wings having sharp trailing edges, 
the position of maximum thickness and base height are 
determined for least pressure drag as functions of a base- 
pressure parameter. Comparison is then made between the 
drag of these wings and similar wings with a sharp trailing 
edge for various aspect ratios and thickness ratios as a 
function of stream Mach number. Calculations of the span- 
wise distribution of drag are included to compare further the 
effect of a base on the drag for different aspect ratios.—-(1.10). 


Span load distributions resulting from constant angle of attack, 
steady rolling velocity, steady pitching velocity, and constant 
vertical acceleration for tapered swepthack wings with stream- 
wise tips. Subsonic leading edges and supersonic trailing edges. 
Margery E. Hannah and K. Margolis. N.A.C.A. Technical 
Note 2831. 
On the basis of the linearised supersonic-flow theory the 
theoretical spanwise distributions of circulation (which are 
proportional to the span load distribution) resulting from 
constant angle of attack, steady rolling velocity, steady 
pitching velocity. and constant vertical acceleration were 
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calculated for a series of thin, swept-back, tapered wings 
with streamwise tips, subsonic leading edges, and supersonic 
trailing edges. The results of the analysis are presented as 
a series of design charts. Some illustrative variations of the 
spanwise distributions of circulation with aspect ratio, taper 
ratio, Mach number, leading-edge sweepback, and axis-of- 
pitch location are also included.—(1.10). 


The flow over two-dimensional aerofoil at incidence in the 

transonic speed range. T. R. Gullstrand. K.T.H. Sweden, 

Aero T.N. 27. 
The transonic integral equation, due to K. Oswatitsch. holds 
only for symmetrical aerofoils without incidence. With the 
help of Green’s theorem, similar integral equations are in 
this report developed for aerofoils at incidence. If there 
are normal shock waves at the aerofoil it cannot be assumed 
that the additional velocity is antisymmetric because the 
shock waves move in different directions on the upper and 
lower surfaces of the aerofoil with change in angle of attack. 
If it is assumed, however, that the additional velocity is 
antisymmetric, except in that part of the flow where the 
shock waves are, approximate equations can be deduced, 
which can be solved in a way similar to that used for the 
other equations, As a result one also gets the extent of 
the shock waves movement with an angle of attack. As 
an example of the method, additional velocity distributions 
are calculated for some Mach numbers for a 6 per cent thick 
symmetrical parabolic are aerofoil. The lift curve slope 
as function of the free stream Mach number is also 
calculated.—(1.10). 


Untersuchungen an einem symmetrischen Tragfliigel mit spaltlos 
angeschlossenem Ruder bei raschen Aenderungen des Ruder- 
ausschlags (ebene Stré6mung). H. Drescher. ©Max-Planck- 
Institut, Géttingen, Mitteilungen Nr. 6.—(1.10). 


The flow over symmetrical aerofoils without incidence at sonic 
speed. T. R. Gullstrand. K.T.H. Sweden. Aero T.N. 24. 


To treat transonic flow it is necessary to take into account 
non-linear terms in the flow equations, which are then very 
difficult to solve. If Green’s theorem is applied, the 
differential equations can be changed to a non-linear integral 
equation. The value of the integrals in the equation can 
be calculated approximately if appropriate assumptions are 
made for the unknown velocity distribution in the region of 
integration. The resulting equation for the velocity distri- 
bution on the aerofoil is then solved by an iteration process. 
In an earlier report by the author the velocity distribution 
on some symmetrical aerofoils at zero lift was calculated in 
the lower transonic range. In the present report the flow 
over some symmetrical aerofoils at zero lift is calculated for 
sonic speed. The solutions are in agreement with available 
experimental results.—(1.10 = 1.2). 


A theoretical discussion of some properties of transonic flow 
over two-dimensional symmetrical aerofoils at zero lift with a 
simple method to estimate the flow properties. T. R. Gull- 
strand. K.T.H. Sweden. Aero T.N. 25. 


With the help of the transonic integral equation it is shown 
that in a region on both sides of sonic velocity the Mach 
number distribution on a symmetrical aerofoil without lift 
is independent of the free stream Mach number. The 
transonic similarity law is discussed and it is shown that in 
predicting the position of the shock wave on the aerofoil the 
law holds only for very thin aerofoils. When there are no 
shocks on the aerofoil the law holds to much larger thickness 
ratios.—(1.10 x 1.2). 


AEROELASTICITY 
Influence of tuned dampers on flexure-aileron flutter. Three 
parts. R. A. Frazer et al. R. & M. 2559. 


In Part I a general theory is developed for the investigation 
of the influence of damping devices on various types on 
flexure-aileron flutter. The numerical applications refer to 
a large transport aircraft, and they are restricted to the case 
of a mass-balanced aileron-carried damper. Part II sup- 
plements Part I and gives results for a partly balanced and 
for a completely balanced aileron-damper system, and Part 
Ilf describes an experimental investigation into the effect 
on flexure-aileron flutter of a tuned damping device attached 
to the aileron.—{2), 
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AIRCRAFT 
See FLIGHT TESTING 
AIRCRAFT DESIGN 
The aerodynamic development of a delta type glider. V. J. 


Smith and E. W. Chambers. 
Note 108. 


This note covers field work and wind tunnel tests on several 
25 and 60 inch span models to provide a stable practical 
glider for use in meteorological research. Various tech- 
niques employed are explained and a series of flight tests 
is described.—(4.1). 


A.R.L. Australia, Aerodynamics 


AIRCRAFT OPERATION 


Ditching tests on a 1/8-size model of the Fokker S-13 crew 
trainer. J. F. Hengeveld. N.L.L. Amsterdam, Report A.1238. 


A dynamically similar model of the Fokker S-13 (1/8 size) 
was subjected to a number of ditching tests at the N.L.L. 
open air ditching station. Conclusions were drawn on the 
effect of the attitude of the plane, on the character of the 
run, and on the influence of such minor variables as flaps, 
propellers, landing gear. etc.—(5.4 x 31.1). 


FLIGHT TESTING 


Harvard 4—flight trials. C.E.P.E., Canada, Report 988. 


This report was prepared to give handling and performance 
data on Harvard 4 aircraft equipped with a propeller spinner, 
an auxiliary fuel tank and a large panel canopy. The 
handling characteristics were fully investigated but no effort 
was made to reduce any of the performance figures to 
standard conditions.—(13.1 = 3.11). 


HYDRODYNAMICS 


Motion of a cylinder under the surface of a heavy fluid. 
Sretensky. N.A.C.A. Technical Memorandum 1335. 


General equations are derived for the approximate solution 
of the flow of a heavy infinite depth about a submerged 
circular cylinder where the centre of the cylinder is sub- 
merged 1:5 diameter or greater. Then by replacing the 
cylinder by a dipole of a certain strength an equation is 
derived for the wave resistance of a circular cylinder moving 
in a fluid of finite depth.—(17.1). 


EON. 


Tank tests on a jet-propelled boat-seaplane fighter (Saunders- 
Roe E6/44). G.L. Fletcher. R. & M. 2718. 


This report describes investigations into porpoising stability. 
water resistance. and seaworthiness of the hull design of 
the E6/44.—(17.2). 


A theoretical and experimental investigation of the effects of 
yaw on pressures, forces, and moments during seaplane landings 
and planing. R. F. Smiley. N.A.C.A. Technical Note 2817. 


A theory for the side force, rolling moment, yawing moment. 
and pressure distribution during the yawed landings and 
planing of seaplanes was developed. For the special case 
of the straight-sided wedge without chine immersion, the 
results of the theoretical analysis are presented in the form 
of generalised curves covering all step landing conditions. 
Experimental impact and planing data are presented for a 
prismatic wedge having an angle of dead rise of 22°5° and 
are shown to be in reasonable agreement with the theoretical 
predictions.—{17.2). 


MATERIALS 
Part VI. 


Sandwich construction and core materials. 
W. J. Pullen, etal. R. & M. 2687. 


A range of struts each consisting of “Balsolite” filler sand- 
wiched between two faces of | in. thick birch plywood was 
tested in order to assess the efficiency of Balsolite 16 as a 
stabiliser in sandwich structures. It is concluded that this 
material compares favourably with other low density 
materials when used as a stabiliser. Modification of the 
material, namely the use of transverse and longitudinal tubes 
alternately. does not appear to be beneficial —({21.5 x 33.2.4). 
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of claw-type yaw probes were determined. Satisfactory ] 
combination survey probes for sensing static and total 
pressure and direction of flow in one or two planes were 
devised.—(31.1 x 1.5). 


MECHANICAL ENGINEERING 


Experimental investigation cf eccentricity ratio, friction, and 
oil flow of short journal bearings. G. B. DuBois and F. W. 
Ocvirk. N.A.C.A. Technical Note 2809. 


An experimental investigation was conducted to obtain STRUCTURES 

performance data on bearings of length-diameter ratios of See also MATERIALS 

1, 1/2, and 1/4 for comparison with theoretical curves. A 

1-375-inch-diameter bearing was tested at speeds up to An analysis of normal accelerations and airspeeds of one type 


6.000 r.p.m. and with unit loads from 0 to 900 Ib./in.. A 
method of plotting oil flow data as a single line is shown. 
Methods are also discussed for approximating a maximum 


of twin-engine transport airplane in commercial operations over 
northern transcontinental route.  R. Steiner. N.A.C.A. 


bearing temperature and evaluating the effect of deflection Technical Note 2833. 
or misalignment on the eccentricity ratio at the ends of the Normal acceleration and air speed data obtained for one > 
bearings.—(23.1). type of twin-engine transport aeroplane in commercial 
operations over a northern trans-continental route are 
analysed to determine the gust and gust-load experiences 


Investigation of 75-millimeter-hore deep-groove ball bearings 
under radial load at high speeds. \—-Oil-flow studies. Z. M. of the aeroplane. The acceleration increments experienced 
Nemeth, E. F. Macks and W. J. Anderson. N.A.C.A. Technical equalled or exceeded the limit-gust-load factor, on the 
Note 2841. average, twice (once positive and once negative) in about 
; 7-5x 10° flight miles. and an effective gust velocity of 30 

At an oil inlet temperature of 100 F. two methods of 


ft./sec. was equalled or exceeded twice in about 1 x 10° 
bearing lubrication, single jet and puddling. were investi- 


A 
gated. In this investigation, 75-millimeter-bore ball bearings 
108 Buckling of low arches or curved beams of small curvature. 
of fom ‘The Y. C. Fung and A. Kaplan. N.A.C.A. Technical Note 2840. 
operating variables on the oil flow through the bearing and A general solution, based on the classical buckling criterion. 
upon the bearing outer-and inner-race temperatures are is given for the problem of buckling of low arches under 
discussed. —(23.1). a lateral loading acting toward the centre of curvature. 
The effects of initial end thrust and axial and lateral elastic 
support are discussed. The buckling load based on the 
energy criterion of Karman and Tsien is also calculated. 
Results for both the classical and the energy criteria are 
compared with experimental results.—-(33.2.1 x 33.2.4). 


POWER PLANTS 


Performance calculations for a doubdle-compound  turbo-jet 
engine of 12:1 design compressor pressure ratio. D. H. 
Mallinson and W. G. E. Lewis. R. & M. 2645. A theoretical and experimental investigation of the influence 
This report describes a theoretical investigation using con- of temperature gradients on the deformation and burst speeds 
ventional component characteristics to discover that division of rotating disks. P. 1. Wilterdink, A. G. Holms and S. S. 
of work between the low and high-pressure compressors of Manson. N.A.C.A. Technical Note 2803. 
equilibrium operation over the normal engine speed range. influence ot 
zi a recently developed method of calculating plastic flow in 


Having decided in favour of a sig mg ratio of 3:1 in the 
low-pressure compressor and 4:1 in the other, a study is 
then _made using more realistic ‘compressor characteristics 


discs by comparing the calculated results with experimental 
observations. Short-time spin tests on parallel-sided, 10-inch- 
diameter discs were conducted under conditions that sub- 


wo determing, the probable performance of such an engine the iss to. range of femperatres from. 70” t 
temperature is 900 deg. C (1,173 deg. K). The equilibrium 1440'F. The agreement between the theoretical and 
running conditions of the engine are investigated with special experimental results was good over ie “fanae OF temper: 
ature conditions investigated.—(33.2.2 x 33.2.4). 
ge reference to the problems introduced by the double- 
i compound type of design.—‘27.1). Lateral instability of rectaneular beams of strain hardening 
material under uniform bending. W. H. Wittrick. A.R.L. 
RESEARCH Australia. Report SM.A98. 
The problem considered in this paper is that of the lateral 
Auxiliary equipment and techniques for adapting the constant- Cross- — under the bending 
semperature anémometer to specific problems in airflow | 
measurements. J. C. Laurence and L. G. Landes. N.A.C.A. of stress-strain curve. Numerical calculations are made for 
Technical Note 2843. two different materials, one representative of an aluminium f 
The constant-temperature hot-wire anemometer amplifier and alloy, the other of annealed mild steel, and the two materials ff | 
accessories have been developed to provide an instrument are compared.—(33.2.2). 
with wide frequency response, good stability. and ease of 
Operation. Auxiliary equipment has been destioned to Evaluation of the theory on the post-buckling behaviour of 
provide heating currents for large wires. to make average- stiffened, flat, rectangular plates subiected to shear and normal 
square computations, and to make double-correlation loads. W.K. G. Floor and T. J. Burgerhout. N.L.L. Amster- 
using this equipment to study periodic phenomena such as : = 
surge, peo stall, and ube ‘ane’ ah centrifugal- and The theory developed in N.L.L. Report S.295 has not been 
axial-flow compressors. The application of the equipment presented in such a form that it could be applied in a 
to the study of non-periodic phenomena such as intensity. convenient way to the majority of practical problems. The 
scale, and spectra of isotropic turbulence is also discussed. numerical evaluation given covered only the special case of 
(31.1). pure shear load acting upon the combination of plate and Is 
stiffeners. Additional numerical calculations were therefore Ff T 
Several combination probes for surveying static and total pe 
and flow direction. Schulze, of the determination of the average stresses and strains or | ey 
and J. R. Erwin, N.A.C.A. Technical Note 2830. the required stiffener cross sections of specified stiffened f mT 
An investigation has teen conducted to provide a basis for plates subjected to specified compressive and shear loads. 
design of combination probes intended to survey static and The case of relatively large compressive loads acting in the 
ee total pressure and direction of flow with special reference plane of the plate in the lateral direction could not fully 
Bie acct to subsonic turbo-machine testing. From the results of be covered because the assumed waveform of the buckles 
this investigation. the factors which determine the sensitivity in the plate proved to be inadequate in this range.—(33.2.4). = 
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UIDED WEAPONS DEVELOPMENT offers good oppor- 

tunities for Senior and Junior Electrical, Electronic, Radio 
and Mechanical Engineers and Draughtsmen, Aerodynamicists, 
Technical Authors, and Computors (female). 

Apply quoting reference G.W. and giving particulars, quali- 
fications and experience, to the Employment Manager. Vickers- 
Armstrongs Limited (Aircraft Section), Weybridge, Surrey. 


DESIGN DRAUGHTSMEN 

Vacancies exist for Senior and Junior Design Draughtsmen 
in the Engineering and Research Departments of the Company. 

Applicants with previous experience in the design of aircraft 
instruments, electrical/electronic equipment and small mechan- 
isms will be preferred, but general technical ability will be the 
first consideration. Permanent progressive posts, with super- 
annuation benefits. 

Apply to Personnel Manager, Smiths Aircraft Instruments 
Ltd.. Bishops Cleeve, near Cheltenham, quoting references, full 
particulars of previous experience, and salary required. 


LECTRICAL DRAUGHTSMEN required for design of 

special aircraft electrical equipment. 

Apply Employment Manager. Vickers-Armstrongs Limited 
(Aircraft Section), Weybridge. Surrey. 


HE ENGLISH ELECTRIC CO. LTD., LUTON, have a 
vacancy for a Senior Aerodynamics Engineer for Guided 
Missile work involving both technical and administrative res- 


ponsibility. Candidates must have a good degree or equivalent 
and be interested in pre-project, performance and assessment 
aspects. Previous experience of supersonics, stability and per- 


formance an advantage. Salary according to experience and 
qualifications. Please write giving full details mentioning 
reference 1078 to The English Electric Company Ltd.. 336/337 
Strand. London. W.C.2 


(Continued) 


DOWTY 
EQUIPMENT LIMITED 


require 
FIRST-CLASS SOUND 
PRACTICAL ENGINEERS 
preferably with some experience of design. development 
and testing of jet engine accessories, particularly flow, 
acceleration and speed governing controls, burners and 
pumps. Exceptional prospects of advancement. 
SENIOR DRAUGHTSMEN 
for development work on Fuel Systems. Experience on 
light hydraulic or pneumatic controls acceptable. 
SENIOR DRAUGHTSMEN 
for Hydraulic and Undercarriage Department. 
Also 
JIG AND TOOL 
DRAUGHTSMEN 
STRESSMAN 


and 
TECHNICAL WRITER 
Important positions with prospects. 

The Company's conditions are exceptionally good; there 
is a realistic pension scheme and the environment is ideal. 
The Company employs a Housing Officer who will render 
every assistance. 
Write, preferably in tabulated form, to: 

PERSONNEL MANAGER 

DOWTY EQUIPMENT LIMITED 
CHELTENHAM 
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SENIOR TECHNICIAN 


The Armaments Divisicn of Sir W. G. Armstrong 
Whitworth Aircraft Ltd.. Coventry, 


Technician to take technical control of Guided Weapon 


require a Senior 


Triais Research and Development in South Australia. 


Applicants should have a Degree in Science or 


Engineering. or equivalent qualifications and consider- 
able experience in Project Development. 
Commensurate salary with excellent prospects. free 


passage for applicant and family. 


Applicants should preferably be between 35 and 45 
years of age. 
Applications to: 
The Chief Engineer, Armaments Division, 
Sir W. G. Armstrong Whitworth Aircraft Ltd., 
Coventry. 


VICKERS-ARMSTRONGS 
LIMITED 
(AIRCRAFT SECTION) 
WEYBRIDGE 


Have vacancies for: 


AIRCRAFT DESIGN DRAUGHTSMEN 
(Senior, Intermediate and Junior) 


RADIO and ELECTRICAL DRAUGHTSMEN 
STRUCTURAL DRAUGHTSMEN 
STRESSMEN 
AERODYNAMICISTS 
SENIOR WIND TUNNEL TECHNICIANS 
WEIGHTS ENGINEERS 
JiG AND TOOL DRAUGHTSMEN 


MECHANICAL TEST ASSISTANTS 
(Electrical and Metallurgical) 


Apply: 
EMPLOYMENT MANAGER 
VICKERS-ARMSTRONGS LIMITED 


(Aircraft Section) 
WEYBRIDGE SURREY 
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UNIVERSITY OF SOUTHAMPTON—Dep: of Aeronautical E 


ESEARCH ASSISTANT required to aid and_ initiate 

researches in the Department of Aeronautical Engineer- 
ing. Applicants. who must have a University degree, will be 
expected to proceed to a higher degree as a result of their work. 
Salary according to qualifications. Applications (3 copies) with 
names of referees and copies of testimonials, should be sent 
to the Secretary and Registrar. The University, Southampton. 


ERODYNAMICIST required by The English Electric Co. 

Ltd., Luton, for work on supersonics with particular 
application to the study of oscillating aerofoil theory. Candi- 
dates must have a good degree, aptitude for mathematics and 
experience in theoretical aerodynamics. Post Graduate quali- 
fications an advantage. Salary according to experience and 
qualifications. Write giving full details. mentioning reference 
457D, to The English Electric Company Ltd., 336/337 Strand, 
London, W.C.2. 


MARKS SECTION 


ACCLES & POLLOCK LTD 


AEROPLANE & MOTOR ALUMINIUM CASTINGS LTD. 
ERDINGTON, BIRMINGHAM. 


ALUMINIUM CAST AUXILIARY DRIVE GEAR BOX 


BIRMINGHAM ALUMIN!IUM CASTING (1903) CO. LTD 


BLACKBURN G GENERAL AIRCRAFT LTD. 


AIRCRAFT MATERIALS LTD. 


STRUCTURAL MATERIALS 
and COMPONENTS 


BOULTON PAUL AIRCRAFT LTD. 


AUTOMOTIVE PRODUCTS CO. LTD. 
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